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Abstract 

Recent progress achieved in the solution of the problem of confinement in various (non-)Abelian 
OA " gauge theories by virtue of a derivation of their string representation is reviewed. The theo- 
^ ' ries under study include QCD within the so-called Method of Field Correlators, QCD-inspired 
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On '. for an analytical description of this phenomenon, is illustrated by an evaluation of field correlators 
^ ■ in these models. 
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1 Introduction 



1.1 General Ideas and Motivations 

Nowadays, there is no doubt that strong interactions of elementary particles are adequately de- 
scribed by Quantum Chromo dynamics (QCD) (see Ref. for recent monographs). Unfortu- 
nately, usual field-theoretical methods are not adequate to this theory itself. That is because in 
the infrared (IR) region, the QCD coupling constant becomes large, which makes the standard 
Feynman diagrammatic technique in this region unapplicable. However, it is the region of the 
strong coupling, which deals with the physically observable colourless objects (hadrons), whereas 
the standard perturbation theory is formulated in terms of coloured (unphysical) objects: quarks, 
gluons, and ghosts. This makes it necessary to develop special techniques, applicable for the 
evaluation of effects beyond the scope of perturbation theory. The latter are usually referred to 
as nonperturbative phenomena. Up to now, those are best of all studied in the framework of the 
approach based on lattice gauge theory 0, which provides us with a natural nonperturbative reg- 
ularization scheme. Various ideas and methods elaborated on in the lattice field theories during 
the two last decades, together with the development of algorithms for numerical calculations and 
progress in the computer technology, have made these theories one of the most powerful tools for 
evaluation of nonperturbative characteristica of QCD (see Ref. for a recent review). However, 
despite obvious progress of this approach, there still remain several problems. Those include e.g. 
the problem of simultaneous reaching the continuum and thermodynamic limits. Indeed, physi- 
cally relevant length scales lie deeply inside the region between the lattice spacing and the size 
of the lattice. However, due to the asymptotic freedom of QCD, in the weak coupling limit, not 
only the lattice spacing, but also the size of the lattice (for a fixed number of sites) becomes 
small, as well as the region between them. However, in order to achieve the thermodynamic limit, 
the size of the lattice should increase. This makes it necessary to construct large lattices, which 
in particular leads to the technical problem of slowdown of simulations on them. As far as the 
problem of reaching the continuum limit alone is concerned, recently some progress in the solution 
of this problem has been achieved by making use of the conception of improved lattice actions 
in Ref. Another problem of the lattice formulation of QCD is the appearance of so-called 
fermion doublers {i.e., additional modes appearing as relevant dynamical degrees of freedom) in 
the definition of the fermionic action on the lattice due to the Nielsen-Ninomiya theorem [0. 
According to this No-Go theorem, would we demand simultaneously hermiticity, locality, and 
chiral symmetry (which will be discussed later on) of the lattice fermionic action, the doublers 
unavoidably appear, which means that all these three physical requirements cannot be achieved 
together. This makes it necessary to introduce fermionic actions which violate one of these prop- 
erties {e.g. Wilson fermions |^, violating chiral symmetry for a finite lattice), checking afterwards 
lattice artefacts associated with a particular choice of the action. Notice however, that recently a 
significant progress in the solution of this problem has been achieved (for a review see and Refs. 
therein). Finally, there remains the important problem of reaching the chiral limit, which becomes 
especially hard if one accounts for dynamical fermions. That was just one of the reasons why the 
main QCD calculations on large lattices have been performed in the quenched approximation, i.e., 
when the creation/ annihilation of dynamical quark pairs is neglected. 

All these problems together with the necessity of getting deeper theoretical insights into non- 
perturbative phenomena require to develop analytical nonperturbative techniques in QCD and 
other theories displaying such phenomena. This is the main motivation for the present work. 
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It is worth realizing that unfortunately up to now a systematic way of analytical investigation of 
nonperturbative phenomena in QCD is still lacking. Instead of that, there exist various approaches, 
enabling one to take them into account phenomenologically for describing hadron interactions. 
Those approaches include e.g. the potential and bag models |10|, |1T|, the large- expansion 
methods [0], the effective Lagrangian approach |13], the QCD sum rule approach [14|, and later 
on its generalization, the so-called Method of Field Correlators [T^, [T^, |TB|. 

The most fundamental problem associated with the IR dynamics of QCD, which is known 
to be one of the most important problems of modern Quantum Field Theory, is the problem 
of explanation and description of confinement (for a review see e.g. [|l9t [18|). In general, by 
confinement one implies the phenomenon of absence in the spectrum of a certain field theory of the 
physical |in) and |out) states of some particles, whose fields are however present in the fundamental 
Lagrangian. It is this phenomenon in QCD, which forces quarks and antiquarks to combine into 
colourless hadronic bound states. The important characteristic of confinement is the existence of 
string-like field configurations leading to a linearly rising quark-antiquark potential as expressed 
by the Wilson's area law (see the next Subsection). Such a field configuration emerging between 
external quarks is usually referred to as QCD string. In the present review, we shall demonstrate 
that the properties of this string can be naturally studied in the framework of the Method of Field 
Correlators and derive by virtue of this method the corresponding string Lagrangian. Notice, that 
the advantage of the Method of Field Correlators is that it deals directly with QCD, and therefore 
enables us to express the coupling constants of this Lagrangian in terms of the fundamental QCD 
quantities, which are the gluonic condensate and the so-called correlation length of the vacuum. 
This approach will also allow us to incorporate quarks and derive a Hamiltonian of the QCD string 
with quarks in the confining QCD vacuum. 

Another fundamental phenomenon of nonperturbative QCD is the spontaneous breaking of 
chiral symmetry, i.e., the U (Nf) x U (A^j)-symmetry of the massless QCD action. Indeed, though 
one could expect this symmetry to be observed on the level of a few MeV, it does not exhibit 
itself in the hadronic spectra. Were this symmetry exact, one would expect parity degeneracy of 
all hadrons, whereas in reality parity partners are generally split by a few hundred MeV. Such a 
spontaneous symmetry breaking has far-reaching consequences. In particular, it implies that there 
exist massless Goldstone bosons, which are identified with pions. The signal for chiral symmetry 
breaking is the appearance of a nonvanishing quark condensate, which plays an important role in 
many nonperturbative approaches [2C, 21]. 

The non-Abelian character of the gauge group SU{3) makes it especially difficult to study 
the problems of confinement and chiral symmetry breaking in the QCD case. To explain the 
mechanisms of these phenomena microscopically, a vast amount of models of the QCD vacuum 
has been proposed (see |jl8[ for a recent review). Those are based either on an ensemble of classical 



field configurations {e.g. instantons |22, 23], see j24] for recent reviews) or on quantum background 



fields ]p5| , Pq] . The most general demand made on all of them was to reproduce two characteristic 
quantities of the QCD vacuum, which are nonzero quark ]^ and gluon Jjl4| condensates, related to 
the chiral symmetry breaking and confinement, respectively. However, at least the semiclassical 
scenario possesses several weak points. First of all, since the topological charge of the QCD 
vacuum as a whole is known to vanish, this vacuum cannot be described by a certain unique 
classical configuration, but should be rather built out of a superposition of various configurations, 
e.g. instantons and antiinstantons. However, such a superposition already does not satisfy classical 
equations of motion and is, moreover, unstable w.r.t. annihilation of the objects with the opposite 
topological charge. Secondly, in order to reproduce the phenomenological gluon condensate JlTl 
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classical configurations should be dense packed (about one configuration per (fm)^), which leads 
to a significant distorsion of the solutions corresponding to these configurations within the original 
superposition ansatz ||2^ [[ And last but not least, a further counterargument against semi classical 
models of the QCD vacuum is that not all of them, once being simulated in the lattice experiments, 
yield the property of confinement (see discussion in Ref. 



Another natural way of investigation of the nonperturbative phenomena in QCD might lie in 
the simplification of the problem under study by considering some solvable theories displaying 
the same type of phenomena. In this way, the problem of chiral symmetry breaking is best of all 
analytically studied in the so-called Nambu-Jona-Lasinio (NJL) type models, which are models 
containing local four-quark interactions [pOj ^ 31, ^ (see Ref. [^] for a recent review). These 



models lead to a gap equation for the dynamical quark mass, signalling spontaneous breaking of 
chiral symmetry. After applying the so-called bosonization procedure (which can be performed 
either by making use of the standard Hubbard- Stratonovich transformation or within the field 
strength approach [Q) as well as a derivative expansion of the resulting quark determinant at 
low energies, this leads to the construction of nonlinear chiral meson Lagrangians |]TB|. The 
advantage of the latter ones is that they summarize QCD low-energy theorems, which is the 
reason why these Lagrangians are intensively used in the modern hadronic physics [^, The 
techniques developed for NJL type models have been in particular applied to the evaluation of 
higher-order derivative terms in meson fields |^T], |3^ , which enabled one to estimate the structure 
constants of the effective chiral Lagrangians introduced in Ref. |3^. Furthermore, in this way in 
Refs. [^, ^ ^ ^ it has been demonstrated that the low-energy properties of light pseudoscalar, 
vector, and axial- vector mesons are well described by effective chiral Lagrangians following from 
the QCD-motivated NJL models. In addition, the path-integral bosonization of an extended NJL 
model including chiral symmetry breaking of light quarks and heavy quark symmetries of heavy 
quarks has been performed (see the second paper of Ref. for a review), which yielded the 
effective Lagrangians of pseudoscalar, vector, and axial-vector D ot B mesons, interacting with 
light TT, p, and ai mesons. 

As far as the theories possessing the property of confinement are concerned, those firstly in- 
clude compact QED and the 3D Georgi-Glashow model |jl9|] and, secondly, the so-called Abelian- 
projected theories |^ |. In the present review, we shall concentrate ourselves on the confining 
properties of the above mentioned non-supersymmetric theories. In all of them, confinement oc- 
curs due to the expected condensation of Abelian magnetic monopoles, after which the vacuum 
structure of these theories becomes similar to that of the dual superconductor (the so-called dual 
Meissner scenario of confinement) ^0|. Such a vacuum then leads to the formation of strings 
(flux tubes) connecting external electric charges, immersed into it. These strings are dual to the 
(magnetic) Abrikosov-Nielsen-Olesen strings The latter ones emerge as classical field config- 



urations in the Abelian Higgs Model, which is the standard relativistic version of the Ginzburg- 
Landau theory of superconductivity. It turns out that the properties of electric strings in the 
dual superconductor are similar to the ones of the realistic strings in QCD, which connect quarks 
with antiquarks and ensure confinement. It is worth noting that this analogy based on the 't 
Hooft-Mandelstam scenario led to several phenomenological dual models of QCD 



42 



(see Ref. 

for a review). Thus, the properties of the QCD string can be naturally studied in the framework 
of the Abelian projection method. Moreover, this approach turns out to provide us with the 



^Recently, some progress in the solution of this problem has been achieved in Ref. |£8|. 

^In what follows, we shall not consider recently discovered supersymmetric theories, also possessing the property 
of confinement [B9| . 
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representations for the partition functions of effective Abelian-projected theories of the SU{2)- 
and SU (3)-gluodynamics in terms of the integrals over string world-sheets. Such an integration, 
which is absent in the approach to the QCD string based on the Method of Field Correlators, 
appears now from the integration over the singular part of the phase of the magnetic Higgs field. 
The reformulation of the integral over the singularities of this field into the integral over string 
world-sheets is possible due to the fact that such singularities just take place at the world-sheets. 
In particular, an interesting string picture emerges in the S'[/(3)-gluodynamics, where after the 
Abelian projection there arise three types of magnetic Higgs fields, leading to three types of strings, 
which (self) interact via the exchanges of two massive dual gauge bosons. An exact procedure of 
the derivation of the string representations for the partition functions of Abelian-projected theo- 
ries in the language of the path-integral, the so-called path-integral duality transformation, will 
be described in details below. In the framework of this approach, we shall also investigate field 
correlators in the Abelian-projected theories and find them to parallel those of QCD, predicted by 
the Method of Field Correlators and measured in the lattice simulations. Topological properties 
of Abelian-projected theories will be also discussed. Furthermore, the effects brought about by 
the summation over the grand canonical ensembles of small vortex loops, built out of the paired 
electric Abrikosov-Nielsen-Olesen strings, in these theories will be studied. In particular in the 
dilute gas approximation, the effective potential of such vortex loops will be derived and employed 
for the evaluation of their correlators. Besides that, we shall study the string representation and 
field correlators of compact QED in 3D and 4D. Notice that due to the absence of the Higgs field 
in this theory, the integration over the string world-sheets is realized there in another way than 
in Abelian-projected theories. Namely, it results from the summation over the branches of the 
multivalued effective monopole potential, which turns out to have the same form as the 3D ver- 
sion of the above discussed potential of vortex loops. Finally, similar forms of the string effective 
actions in QCD, Abelian-projected theories, and compact QED will enable us to elaborate for all 
these theories a unified method of description of the string world-sheet excitations, based on the 
methods of nonlinear sigma models, known from the standard string theory. 

It is worth reahzing, that the string theories, to be derived below, should be treated as effective, 
rather than fundamental ones. The actions of all of them turn out to have the form of an interaction 
between the elements of the string world-sheet, mediated by certain (nonperturbative) gauge field 
propagators. Being expanded in powers of the derivatives w.r.t. world-sheet coordinates, these 
actions yield as a first term of such an expansion the usual Nambu-Goto action. The latter 
one is known to suffer from the problem of conformal anomaly in D 7^ 26 appearing during its 
quantization, which will not be discussed below. It this sense, throughout the present review, we 
shall treat the obtained string theories as effective 4D ones. It is also worth noting that within 
our approach only pure bosonic strings without supersymmetric extensions appear. As far as 
superstrings are concerned, during the last fifteen years, a great progress has been achieved in 



their development (see e.g. for comprehensive monographs). Among the achievements of the 
superstring theory it is worth mentioning such ones as the calculation of the critical dimension 
of the space-time, inclusion of gravity in a common scheme, and, presumably, the absence of 
divergencies for some of these theories. The aim of all the superstring theories is the unification of 
all the four fundamental interactions. In another language, one should eventually be able to derive 
from them both the Standard Model and gravity, whereas all the auxiliary heavy modes should 
become irrelevant. Therefore, the final strategy of superstring theories is a derivation of the known 
field theories out of them. Contrary to this ideology, the aim of the present review is the derivation 
of effective string theories from gauge field theories possessing string-like excitations. As it has 
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been discussed above, such string-like field configurations naturally appear in the confining phases 
of gauge theories. 

Another possible direction of investigation of confinement and chiral symmetry breaking in 
QCD is based on a derivation of self-coupled equations for gauge-invariant vacuum amplitudes 
starting directly from the QCD Lagrangian and seeking for solutions allowing for these proper- 
ties H3, H^. Recently, this approach turned our to be quite useful for the investigation of 



the problem of interrelation between these two phenomena [47]. Once such an interrelation takes 



place, there should exist a relation between quark and gluon condensates as well, which has just 



been established in Ref. We shall briefly demonstrate the method of derivation of such a 

relation later on. 

The organization of the review is as follows. In the next Subsection of the Introduction, we shall 
introduce the main quantitative parameters for the description of confinement in gauge theories 
and quote the criterion of this phenomenon in the sense of Wilson's area law. This criterion will 
then serve as our starting point in a derivation of certain string effective actions in various gauge 
theories. In the last Subsection, we shall consider theoretical foundations of the Method of Field 
Correlators. Section 2 is devoted to a derivation and investigation of the QCD string effective 
action within this method. In Section 3, we investigate the problem of string representation of 
QCD from the point of view of Abelian-projected theories and demonstrate a correspondence 
between the Abelian projection method and the Method of Field Correlators. In Section 4, we 
study the string representation and vacuum correlators of compact QED in 3D and 4D. All these 
investigations eventually bring us to the conclusion that both QCD within the Method of Field 
Correlators, Abelian-projected theories, and compact QED have similar string representations. 
Such an observation then enables us to consider strings in these theories from the same point of 
view and elaborate for them a unified mechanism of description of string excitations. Besides that 
in Section 4, by virtue of the techniques developed for the investigation of the grand canonical 
ensemble of monopoles in compact QED, collective effects in similar ensembles of topological 
defects emerging in the Abelian-projected theories will be studied. This will enable us to improve 
on the calculations of the field strength correlators, performed in Section 3. Finally, we discuss the 
main results summarized in the review as well as possible future developments in the Conclusion 
and Outlook. In five Appendices, some technical details of transformations performed in the main 
text are outlined. 

1.2 Wilson's Criterion of Confinement and the Problem of String Rep- 
resentation of Gauge Theories 

As a most natural characteristic quantity for the description of confinement one usually considers 
the so-called Wilson loop. For example in the case of QCD, this object has the following form 



{W{C)) = j^(^i^P exp ^zgfA,dx}j 



which is nothing else, but an averaged amplitude of the process of creation, propagation, and 
annihilation of a quark-antiquark pair. In Eq. (|1]), stands for the vector-potential of the 
gluonic field |, g is the QCD coupling constant, C is a closed contour, along which the quark- 

■^From now on in the non- Abelian case, = yl^t°,a = 1, . . . — 1, where = (i")*-' is the Hermitean 
generator of the colour group in the fundamental representation, whereas in the Abelian case A^^ is simply a vector 
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antiquark pair propagates, P stands for the path-ordering prescription, which is present only in 
the non-Abehan case, and the average on both sides is performed with the QCD action 0. 

In order to understand why this object really serves as a characteristic quantity in QCD, let 
us consider the case when the contour C is a rectangular one and lies for concreteness in the 
(xi, t)-plane. Let us also denote the size of C along the t-axis as T, and its size along the xi-axis 
as R. Then such a Wilson loop in the case T ^ R is related to the energy of the static {i.e., 
infinitely heavy) quark and antiquark, which are separated from each other by the distance R, by 
the formula 

(Wr^t) ~ e-^°(^>^, T » R. (2) 

In order to get Eq. @, let us fix the axial gauge = 0, so that only the segments of the 
rectangular contour C parallel to the Xi-axis, contribute to {Wrxt)- Denoting 



R 



(3) 



where we have omitted for shortness the dependence of \E'jj on X2 and X3, we are not interested 
in, we get 

{Wn.T) = ^^{^^,mUT)). (4) 
Inserting into Eq. a sum over a complete set of intermediate states |^) {n\ = 1, we get 

n 

(Wr.t) = ^E(%(0) I^)(^|4(^)> = W^mM\n)fe-'^-^, (5) 

where En is the energy of the state \n). At T ^ 00, only the ground state with the lowest energy 
survives in the sum over states standing in Eq. (^, and we finally arrive at Eq. (Q). 

The energy Eq{R) in Eq. (0) includes a /^-independent renormalization of the mass of a heavy 
(anti) quark due to its interaction with the gauge field. To the first order in g'^, up to a colour 
factor, it is the same as in QED |^ and reads 



^E^^ss = C2/-, (6) 
Ana 

where C2 = ^2N^ stands for the Casimir operator of the fundamental representation, and a — > is 
a cutoff parameter {e.g. lattice spacing). The difference E{R) = Eo{R) — AE^^^g therefore defines 
the potential energy of the interaction between a static quark and antiquark. In particular, the 
exponential dependence of the Wilson loop on the area of the minimal surface Smm. [C] encircled 
by the contour C, which we shall denote by iSminJC*]!, 

(H^(C)) -.e-'^l^-'-l^ll (7) 



potential. 

^In what follows, we call the object defined by Eq. (|^) for brevity a "Wilson loop", whereas in the literature it 
is sometimes referred to as a "Wilson loop average" . 

^Throughout the present review, all the investigations will be performed in the Euclidean space-time. 
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(the so-called area law behaviour of the Wilson loop) corresponds to the linearly rising potential 
between a quark and an antiquark, 



E{R) = aR. 



This is the essence of the Wilson's criterion of confinement 

The coefficient a entering Eqs. and is called string tension. This is because the gluonic 
field between a quark and an antiquark is contracted to a tube or a string (the so-called QCD 
string), whose energy is proportional to its length, and a is the energy of such a string per unit 
length. This string plays the central role in the Wilson's picture of confinement, since with the 
increase of the distance R between a quark and an antiquark it stretches and prevents them from 
moving apart to macroscopic distances. 



In order to get an idea of numbers, notice that according to the lattice data |50] the distance 
i?, at which Wilson's criterion of confinement becomes valid, is of the order of 1.0 fm, and the 
string tension is of the order of 0.2 GeV^ (see e.g. |1^). It is worth realizing that the classical QCD 
Lagrangian does not contain a dimensionful parameter of such an order {i.e., of hundreds MeV) 0. 
However, in quantum theory, there always exists a dimensionful cutoff (like the lattice spacing a in 
Eq. (P)), which is related to the QCD coupling constant g through the Gell-Mann-Low equation 



— a 



2 / 1 



da? 



9 pQCB [g 



(9) 



Here /3qcd (g^) stands for the QCD Gell-Mann-Low function, which at the one-loop level reads 



/5qcd {g') 
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167r2' 



(10) 



where Nj is the number of light quarks flavours, whose masses are smaller than 1/a. It is Eq. (|10|), 
which tells us that QCD is an asymptotically free theory, provided that for Nc = S, Nf < 16, 
which indeed holds in the real world. Consequently, the high-energy limit of QCD (the so-called 
perturbative QCD) is similar to the low-energy limit of QED, and the scale parameter following 
from the integration of Eq. (^, 



^QCD = -^exp 



dg^ 
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is measurable in QCD as well as the QED flne-structure constant (=1/137) with the result 
100 MeV < Aqcd < 300 MeV |Q. The phenomenon of the appearance of a dimensionful pa- 
rameter in QCD, which remains flnite in the limit of vanishing cutoff, is usually referred to as di- 
mensional transmutation. All observable dimensionful quantities in QCD {e.g. widths of hadronic 
decays), and, in particular, the string tension are proportional to the corresponding power of 
AqcdQ- Then according to Eqs. ( pUj ) and (|r^), we get 

^E.g. the masses of the hght quarks are of the order of a few MeV. 

^ The aim of all the nonperturbative phenomenological approaches to QCD, such as the Method of Field 
Correlators which will be described in the next Subsection, is to calculate the (7-independent dimensionless ratios 
of these quantities, but not Aqcd itself. 
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which means that all the coefficients in the expansion of the string tension in powers of vanish. 
This conclusion tells us that the QCD string has a pure nonperturbative origin, as well as the 
phenomenon of confinement, which leads to the process of formation of such strings in the vacuum 
itself. 

Throughout the present review, we shall be mostly interested in the models and properties of 
the QCD string and strings in other gauge theories, possessing a confining phase. In the literature, 
this problem is usually referred to as a problem of string representation of gauge theories. 

2 

Contrary to the linearly rising quark- antiquark potential, the Coulomb potential, E{R) oc — ^, 
cannot confine quarks, since in this case the gauge field between them is distributed over the whole 
space. For such a potential, the Wilson loop for the large contour C has the following asymptotic 
behaviour 

{W{C)) e-"°'^^*-^('^\ (13) 



1 



where L{C) = J dsJx'^{s) stands for the length of the contour C, parametrized by the vector- 
^ 

function < s < 1, x^{0) = x^(l). Such a behaviour is called the perimeter law. It is 

dominant at small distances, R < 0.25 fm, where quarks can with a good accuracy be considered 
as separate particles not connected by strings, i.e., in the framework of perturbative QCD. 

It is worth noting that to each order of perturbation theory, it is the perimeter law (piSf ), rather 
than the area law (|^), that holds for the Wilson loops. Because of the ultraviolet divergencies, for 
such a perturbative expansion of Eq. (|I|) in powers of g, one needs a (gauge invariant) regulariza- 
tion. When such a regularization is introduced, the Wilson loop for a smooth contour C {i.e., a 
contour without cusps) takes the form 



{W{C)) = exp 



L{C) 

An a 



{W,,nXC)), (14) 



where (WrenX^)) is finite when expressed via the renormalized charge (^ren.- The exponential factor 
in Eq. (|14|) is due to the renormalization of the mass of a heavy (anti)quark, described by Eq. (^). 
Such a multiplicative renormalization of the smooth Wilson loop has been proved in Refs. |^2|, |5^, 
|5^ . Notice also that if the contour C has a cusp but no self-intersections, then {W{C)) is still 
multiplicatively renormalizable [|55|. Namely, in that case {W{C)) = Z{j) (VFrcn.(C')), where the 
diverging factor Z{j) depends on the cusp angle 7. 

Thus, we conclude that the Wilson loop indeed plays the role of the quantity relevant to the 
description of confinement. Its area law behaviour means that the gauge theory under study {e.g. 
QCD) is in the confining phase, whereas the perimeter law behaviour means that the theory is in 
the Coulomb phase. We also see that due to Eq. (Q), the change of — \n.{W{C)) determines the 
change of the action due to the interaction of an (anti) quark with the gauge field. 



1.3 Method of Field Correlators in QCD: Theoretical Foundations 

The nonlinear character of the QCD action makes it impossible to calculate in this theory the vac- 
uum averages of gauge-invariant quantities in a closed form, e.g. by making use of the path-integral 
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techniques. On the other hand, the standard perturbation theory also becomes unapphcable at 
large distances, where the running coupling constant is around 0.5-1. This makes it necessary to 
develop special approaches to nonperturbative QCD. In the present Subsection, we shall briefly 
describe one of them, the so-called Method of Field Correlators (MFC) in QCD. The idea, which 
lies behind this approach, is that a large class of observables in QCD can be expressed in terms of 
the Wilson loops In this respect, the calculation of some QCD observable consists of two 

steps: the calculation of the Wilson loop for an arbitrary contour and summation of the Wilson 
loops over contours with a certain weight, which is determined by the observable Q For example, 
the Green function of a system consisting of a scalar quark and antiquark with equal masses m, 
which propagates in the QCD vacuum from the initial state {y^y) to the final state (x, x). 



G (x, X] y, y) = (ii (i) (x) $ (x, x) ?/'(x)) (i){y)^ {y, y) ip (y) 



(15) 



where 



'^ix,y) = — P exp 

stands for the parallel transporter factor along the straingt line ^ in the quenched approximation 
reads IBS, 



+00 +00 



G{x,x;y,y) = J ds J dse-"''^'^'^ j Dz j Dzexp l-^ J z^d\-^J fd\ \ {W{C)) . (16) 



In Eq. ([T6|), the dot stands for and the closed contour C consists of the trajectories z^ and z^ 
of a quark and antiquark, z{0) = y, z{s) = x, z{0) = y, z{s) = x, and straight-line pieces, which 
form the initial and final states. One can see that the weight factor in the path-integral standing 
in Eq. (|16D is completely determined by the free theory, and all the dependence on the gauge field 
factors out in the form of the Wilson loop. 

After that, the main idea of the MFC is not to directly evaluate the Wilson loop, but to 
express it via gauge-invariant irreducible correlators (the so-called cumulants) of the gauge field 
strength tensors. Such correlators have been measured in lattice experiments both at large and 
small distances pTl K2[ , which enables one to use them in practical calculations of various physical 



quantities. In order to express the Wilson loop (|l|) via cumulants, one should first make use of 
the non-Abelian Stokes theorem, derived in Refs. []63| , |64| , which yields 



(WiC)) = (tiPexp (^tg I da,,ix{0)F,,ixiO,Xo)^ ^ . (17) 



®It should be mentioned that contrary to the classical theory (where only the electric and magnetic field strengths 
are observable, while the vector potential plays only an auxiliary role for their determination), in the quantum 
theory Wilson loops are observable themselves owing to the Aharonov-Bohm effect . 

^From now on, for quantities including P-ordering along open paths, we adopt the convention that the matrices 
are ordered from the second argument to the first one. 

^"in the following equation, the summation is performed over < v {cf. also Eq. (EO) below). 
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On the R.H.S. of Eq. (]T7|), the integration is performed over an arbitrary surface E bounded by 
the contour C and parametrized by the vector x^{^). Next, ^ = (^^,^^) is a two-dimensional 
coordinate, and 



stands for the infinitesimal surface element, where g{^) is the determinant of the so-called induced 
metric tensor of the surface, defined as 



gabiO = {daX,{mdbX,{0), (19) 



where da = g^s:, a,b = 1, 2. Next, 



V(0 = -^e^\dax,{0){di>x.{0) (20) 



is the extrinsic curvature tensor of the surface, tf^^{$,) = 2, F^^{x,Xo) = $(xo, a;)F^,,(x)$(x, Xq) is 
the covariantly shifted non-Abelian field strength tensor, F^j, = d^A^ — d^A^ — ig [A^, A^], and Xq 
is an arbitrary but fixed reference point, the dependence on which actually drops out. 

The second step in expressing the Wilson loop in terms of the cumulants is to rewrite 
Eq. (jl^) via the cumulant expansion theorem |6^. The aim of this theorem is to express cumu- 



lants, which are usually denoted as {{F^-^^^-^{xi, xq) ■ ■ ■ F^^^^{xn,xo))), through the usual averages 
(xi, Xo) ■ ■ ■ F^„,y„(x„, Xo)). To this end, let us first consider for simphcity the Abelian case, 

where 

{Ff,^u,{Xl,Xo)---Ff,^^^{Xn,Xo)) = (F^ii.i(Xi) • ■ •F^„^„(Xn)) = 



J DA^,r^{A^)F^,,,{x,)---F,^,M, (21) 



with rj (A^) standing for a certain 0(4)-invariant integration measure, and introduce the following 
generating functional for field correlators 

Z [J^,] = ^exp (^ig J d^xJ^,{x)F^,{x)^ ^ . (22) 
The field correlators then obviously have the form 



^JfllUl{Xl) ■ ■ ' ^J^nUn^Xn) 

By virtue of the Taylor expansion. 



(23) 



j=o 



°° 1 /■ /■ ( S' 

Z \J^y\ = 1 + -T / d^^l ■■■ d'^XnJf^.u^iXi) ■ ■ ■ Jf,^u„iXn) JJ ^ 

n=l X'J ^ mui\Xl 



where it has been used that 2[0] = 1, we obtain for the generating functional 



Z \J,.\ = 1 + E ^ / d'x^ ■■■ [ rf^x„J,,,,(xi) ■ ■ ■ J^„.„(xO (F.^^A^i) ■ ■ ■F^„.„(x„)) . (24) 



12 



By definition, an n-th order cumulant is defined as follows 



(25) 



J=0 



In particular, if we specify a current to the form Jfj,u{x) = J da^u{x{C,))6{x — x{^)), where E is a 

s 

certain surface bounded by the contour C, then, by virtue of the Stokes theorem, the partition 
function (^) (where the summation is performed over fi < u) is nothing else, but the Wilson 
loop. The definition (^5[) means that the cumulants of the field strength tensors are the analogues 
of usual irreducible (connected) Green functions, which is quite in line with the interpretation of 
— In {W{C)) as a correction to the free energy due to the interaction of the gauge field with an 
external charged particle (c/. the end of the previous Subsection). Next, following the same steps 
which led from Eq. (ESl) to Eq. (p^), we obtain from Eq. (pSf) 



{WiC))=exp [p^J da,,,,{x,)--- j da^^,^{xu) {{F,,,A^,) ■ ■ ■ F^^,^{x,)))^ . (26) 



Eq. (Pq ) is usually referred to as the cumulant expansion. 

Notice that in Eqs. (p^) and (p6|), it has been assumed that the series on their R.H.S. converge, 
which is true for non-pathological models. To those belongs in particular QCD even at large 
distances, i.e., at large g (see discussion in Ref. |18|). This fact means that the cumulant expansion 
in QCD is indeed a nonperturbative expansion. The important property of the cumulants in such 
models, which distinguishes them from the usual correlators (thus associated with the usual Green 
functions), is that any cumulant decreases vs. the distance between any two points, in which the 
fields in the cumulant are defined. An example of the model where the cumulant expansion 



diverges is the instanton gas |66 



One can now write down the following generating equation for cumulants 

1 + ~~^\ j ^ Xi ■ ■ ■ J d XnJu^ui i^l) ' ' ' Jfini^n i^n) {Ffiiui i^l) ' ' ' Ffj,^u„ {Xn}) 



n=l 



^^p(£^-fr / ^^^i---/ d^^kJ,,^,yAxl) ■ ■ ■ Jf^.u^Xk) {{F^,^^{xi) ■ ■ ■ F^^^^ixk)))] , (27) 



(ig)'' 

\k=l 

which means that varying Eq. (pTf ) several times w.r.t. J^^ and setting then J^^, = 0, one can get 
relations between the cumulants and correlators of various orders. 

In this way, the one- fold variation and setting J^^, = yield {{F^^{x))) = {F^^{x)). Notice 
that due to the 0(4)-invariance of the Euclidean integration measure rj (A^), this average vanishes. 
Next, the two- and three- fold variations and setting J^^, = yield the following relations 

{F^,uM)F^.,u,ix2)) = {{F^,uA^i)F^,.Ax2))) + {F^,,u,M) {F^,u,{x2)) (28) 

and 

{X2)F^.^^^{X3)) = {{Ff,^^^{xi)F^^^^{x2)F^,^^^{xs))) + 
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+ {F^,u^{xi)) {{F^.,u2{x2)F^^^.^{x^))) + {F^^^^{x2)) {F^.^^ixs)) , (29) 

respectively. Eq. (|29|) can be symbolically written as follows 



{Ff,^uAxi)Ff,^u2ix2)F^^^^{x3)) = {{Ff,^^^{xi)F^^^,^{x2)F^^^^{x3))) + 

+ (3) {{F,,,u^{xi)Ff,^^^{x2))) (F^3^3(x3)) + {F^,u^{xi)) {F^^u.Ax2)) {F^^u^ixs)) , (30) 

where the coefficient in curly brackets denotes the number of the terms of the same type, which 
differ from each other only by the order of the arguments (and indices, respectively). Making use 
of this notation, we can next write e.g. the following equation 



{F^,uAxi)F^^^^{x2)F^^^^{x3)F^^^^{x4)) = {{F^^ {X2)F^.^^^{XS)F^^^^{X4))) + 

+ (3) {{F^,.Axi)F^2u,{x2))) {{F^,.,{x3)F^,,,{x^))) + 
+ (4) {{F^,u,{xi)F,^2U2{x2)F^,,,{xs))) {F^,u,{x4)) + 



+ (6) {{F^,,uriXi)F^^^^{x2))) (F^3^3(X3)) (F^,^,(x4)) + (F^i^i(Xi)) {F^^U2{X2)) {F^:,usix3)) {F^^^^{x^)) . 

(31) 

Now one can see that on the R.H.S. of Eqs. (|30|) and ( ^I]) stand with the coefficient equal to unity 
all possible terms, which correspond to various splittings of the set {F^j^^(xi), . . . , F^^t,^{xn)} into 
subsets, so that to every subset corresponds a cumulant. This property holds also for the higher 
correlators, which gives a simple mnemonic rule for calculation of the higher cumulants without 
using Eq. (^). This rule is the essence of the cumulant expansion theorem. 

For the non-Abelian case, the cumulant expansion theorem remains the same except for exclud- 
ing the terms which violate the path-ordering prescription E.g. in the non-Abelian version 
of Eq. (PI), the term 



{{F^,,u^ixi, Xo)Ff,^^.^{x3, Xo))) {F^^U2{x2,Xo)) 

will be absent for a given order of points xi,X2,xs. 

It looks natural to address the question on whether cumulants of various orders are independent 
of each other or there exist any relations between them. Such equations relating cumulants of 
different orders indeed take place. They have been derived by making use of the stochastic 
quantization method [^, ^ in Ref . and further investigated in Refs. [^, |7^ . Alternative 



equations for cumulants following from the non-Abelian Bianchi identities have been proposed in 
Refs. 16^ , Notice that quite recently, there has also been proposed one more approach to a 



derivation of a system of self-coupled equations for cumulants. This has been done in Ref. by 



considering a field of a certain colour in the gluodynamics Lagrangian as the one propagating in 
the background of the fields of all other colours by making use of the background field methods. 
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there have recently been proposed new equations for 
derived by virtue of the stochastic quantization method. Several techniques of 



Finally, besides equations for cumulants 
Wilson loops 



the loop space approach, exploited in a derivation of these equations, have been then applied to 
the solution of the Cauchy problem for the loop equation in 3D turbulence . 

By definition adopted in the theory of random processes, a set of random quantities is called 
Gaussian, provided that all the cumulants of these quantities higher than the quadratic one vanish. 
According to the lattice data |6l|, the stochastic ensemble of fields in QCD can with a good 
accuracy be considered as a Gaussian one, i.e., the bilocal cumulant 



{{F^u{x,xo)Fxp{y,xo))) 



(32) 



is much larger than all the higher cumulants, so that the cumulant expansion converges fastly. 
Due to this property of the QCD vacuum, one can disregard all the cumulants except the bilocal 
one, which leads to the bilocal or Gaussian approximation in the MFC in QCD. However one can 
immediately see that the neglection of all the cumulants higher than the bilocal one in the non- 
Abelian analogue of Eq. (p6D leads to the appearance of the artificial dependences of the R.H.S. 
of this equation on the reference point xq and on the shape of the surface S. 

The first problem can be solved by noting that the bilocal cumulant (|3^) decreases fastly when 



X — y\ ^ Tg, where Tg is the so-called correlation length of the vacuum, Tg ~ 0.13 fm in the 



SU{2)-case and Tg ~ 0.22 fm in the SU {3)-case [|T 



been discussed in the previous Subsection, according to Ref. ^0 



On the other hand, as it has already 
the area law of the Wilson loop 

takes place when its size R is of the order of 1.0 fm. This means that in the confining regime, i.e. 
for the loops of such size, one can in the general case {i.e 
write down the following inequality 



for the dominant amount of cumulants) 



|x — xo| — |y — xol — -R ^ Tg ~ |x — yl . 

According to it, we can with a good accuracy neglect the dependence on the point xq in the bilocal 
cumulant (|3^), i.e., approximate this quantity by the gauge- and translation- invariant cumulant 
as follows 



tr {{F^u{x, Xo)Fxp{y, Xq))) ~ tr {{Fp^{x)<i>{x, y)Fxp{y)^y, x))) . 

As far as the problem of the artificial dependence on the shape of the surface S in the bilocal 
approximation is concerned, it cannot be solved on the basis of simple theoretical arguments, 
but one can choose S to be the surface of the minimal area for a given contour C . As we shall 
see in Subsection 2.1, this enables one to reproduce the area law behaviour of the Wilson loop. 
Since such a surface Smin. is uniquely defined by the contour C , we finally arrive at the following 
expression for the Wilson loop in the bilocal approximation 

{W{C)) = (l^(S^in.)) ^ 




where we have denoted for brevity x = x(^) and x' = x(^'). Due to the colour invariance of 
the Euclidean integration measure in QCD, the bilocal cumulant on the R.H.S. of Eq. ( p3D is 
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proportional to the unity matrix in the fundamental representation, liY^xTVc) ^ind therefore the 
P-ordering is not necessary any more. 

The bilocal cumulant standing in the exponent on the R.H.S. of Eq. (^) can be now 
parametrized by two renormalization-group invariant coefficient functions D and Di as follows 



9' 



((F^,(x)<l>(x, x')Fxp{x')^{x\ x))) = In.^nA {5,xKp - ^^^D ((x - x'f) + 



1 

+ 2 



d , d , 

- x')\6^p - (x - x')p6^x) + Q^ii^ ~ 3;')p'^ma - (x - x')xSpp) 



D, [{x-x'f)^. 



(34) 



The parametrization ( 13^) of the bilocal cumulant is chosen in such a way, that the term containing 
the function Di yields a perimeter type contribution to the Wilson loop (|33|) . Namely, by making 
use of the (ordinary) Stokes theorem, one can prove that |^ 



1 
2 



dcrp^{x) / daxpix' 



d 

{{x - x')x5up - (x - x')p5^x)+ 



dxp 



+^((x - x')p6px - {x- x')xSpp) 



A ((x - x'f) = idxpi dypG ((x - yf) , (35) 



c c 



where 



+ 00 



G(x2) = J d\Di{\). (36) 

In particular, we see that the one-gluon-exchange diagram contribution to the Wilson loop, equal 
to 

exp I -C2-^ i dxp i dy 



c c ^ , 

(which for the contour C without cusps yields the renormalization factor standing on the R.H.S. 
of Eq. (p!^) is contained due to Eq. (|35| ) in the function Di. According to Eq. (|36|), it is equal to 

{^) = C.^.^, (37) 

Eq. ([37| ) is the leading contribution to the function Di. Therefore the effects of renormalization 
of the Wilson loop, discussed at the end of the previous Subsection, are taken into account in the 
MFC by virtue of the function Di. 

Thus, the function Di is nonvanishing already in the order (7^, whereas the function D in this 
order of perturbation theory vanishes. The nonvanishing contributions to the function D emerging 
in the higher orders have been calculated in Ref. |7^. From now on, we shall not be interested 
any more in the perturbative contributions to these functions. The nonperturbative parts of the 



^^From now on, for simplicity, we restore in Eq. ( p^ ) the usual agreement on the summation over all /x and v (A 
and p), not only fi < v. This can always be done by the appropriate normalization of the functions D and Di. 
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functions D and Di have been calculated in the lattice experiments in Refs. , where it has 

been shown that they are related to each other as Di ~ ^D. 



Due to the Lorentz structure standing at the function D in Eq. (|3^ , its contribution to the 
Wilson loop ( |33D cannot be reduced to that of a perimeter type. This function is responsible for 
the area law behaviour of the Wilson loop and gives rise to the QCD string effective action. The 
problems of a derivation of this action from Eq. ( PBD and its further investigation will be the topic 
of the next Section. 



String Representation of QCD in the Framework of the 
Method of Field Correlators 



In the present Section, we shall demonstrate the usefulness of the MFC for the derivation and 
investigation of the local gluodynamics string effective action. Our interpretation of this topic in 

H in Subsection 2.1 



the next three Subsections will mainly follow the original papers |75 



77 



in Subsection 2.2, and in Subsection 2.3. A short review of all these papers can be found in 
Ref. m 



2.1 Gluodynamics String Effective Action from the Wilson Loop Ex- 
pansion 

Let us start with the problem of a derivation of the effective action of the gluodynamics string, 
generated by the strong vacuum background fields in QCD, which ensure confinement and yield a 
dominant contribution to the bilocal cumulant. To this end we mention that, as it has already been 
discussed in the previous Section, the quantity — In {W{C)) is nothing else, but a correction to the 
gluodynamics free energy due to the interaction of a test quark, moving along the contour C, with 
these background fields. This observation enables us to treat this quantity with {W{C)) defined 
by Eq. ( p3D as a nonlocal background-induced gluodynamics string effective action associated with 
the minimal area surface Smin.- In what follows, by the terms "local" and "nonlocal string effective 
action", we shall imply the actions depending on a single string world-sheet coordinate ^ or on 
two such coordinates, respectively. 

Substituting Eq. ( p4D into Eq. (p3D, one can see that the coefficient function D plays the role 
of a propagator of the background field between the points x(^) and which lie on Smin.. 

Though this propagator cannot be completely found analytically, its nonperturbative part can be 
parametrized with a good accuracy as 

D(x) = a.((F;,(0))')e-N/^^ (38) 

where the gluonic condensate ^(f;^(0))^^ is of the order of 0.038 GeV^ This will 

finally enable us to obtain the coupling constants of the few first terms of the resulting local string 
effective action expressed via the correlation length of the vacuum, Tg, and the gluonic condensate. 

In order to derive the desired local effective action, let us proceed with the Taylor expansion 
of the nonlocal action 

^^The colour index "a" should not be confused with the index of the two- vector ^. 



17 



5'efr. (Smin.) = 2 J da^,y{x) J da^v{x)D f- 



X — X' 



7^2 

9 



(39) 



in powers of the derivatives w.r.t. the string world-sheet coordinates ^"'s. Obviously, the nonlocal- 
ity of the initial action will then display itself in the appearance of higher derivatives w.r.t. ^"-^s in 
the final expression for the action. Notice also that in Eq. (PD|), we have explicitly emphasized the 
form of the dependence on Tg, since, as we shall eventually see, the parameter of this expansion 
will be proportional to Tg. 

To perform an expansion, let us first rewrite infinitesimal surface elements on the R.H.S. of 
Eq. (|39D by making use of Eqs. (|T8|), ([T9|) , and (pOf) . Next, we introduce instead of ^'"'s new 
integration variables (^^ 



rp and expand in power series of C^'s the quantities y g{^'), t^^i^'), 

x{^') — x{^), and finally D ^ ^^^2 ^ ■ Such an expansion will automatically be a formal series in 
powers of Tg, e.g. 



(40 -my 



CC'gatiO + TgCC'C{daX,){d,d,X,) + o (t^) 



In what follows, we shall be interested only in the terms in this series not higher than of the fourth 
order in Tg. Then, taking into account that for an arbitrary odd n, J (PCC^ ' ' ' C"D{y) = 0, where 
y = CC^gabiO^ we obtain 

S,s. = 2T^J d'^^)i^2^) J d'CD{y) + T^i^^) -1(9,V(0)(5^V(0) / d'CCC''D{y) + 

+ (^lidadtx,)id,d,x,) + ^idaX,)idMx,)^ J d\ C CXX' D' (y) + 
+ {dax,){d,d,x,)id,x,){d,dfx,) J d\CC'CXWD"iy) 



+2 (S^V^) id,x,){d,d,x,) J d\CC\X'D'{y) + (dad,^)) J d\CC'D{y) 

where ""' stands for the derivative w.r.t. the argument, and we have arranged the terms in powers 
of the derivatives of \J g{C). In order to simplify the integrals over (^'^'s as much as possible, it is very 

useful to fix by reparametrization the conformal gauge for the induced metric, gahiCl = ygiO^ab 
(see e.g. Ref. |T^). After that, by making use of the equality = 2^/1^) '^' becomes possible 
to perform partial integrations over C^'s. It is also useful to replace the ordinary derivatives by 
the covariant ones according to the Gauss- Weingarten formulae 

DaDhXf, = dadbXf, - Tl^dcXf, = Kl^n'^. 

Here, n^'s stand for the unit normal vectors to the string world-sheet, n^^n^^ = 5*-', n^^daX^ = 0, 
i,j = 1, 2, T^f^ is a Christoffel symbol, and i^r*^ is the second fundamental form of the world-sheet. 
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Notice that in the conformal gauge adopted, all the quantities are greatly simplified, e.g. 



r^f, = i (5a 9b + 5lda — Sabd^) ^^y 9i0- particular, the expression for the scalar curvature TZ 
(X^")^ — Kl"'Kl^ of the world-sheet in this gauge takes the form 



7^ 



9i0 



and one can prove the validity of the following equations 



(40) 



9(0 



T''"''idaX^)idMx^) = T'^'"'{DaX^){D,D,DdX^) + 27^ 



T'''>'\DaD,x^){D,DdX^) = ?>{D^DaX^){D'D,x^) - 27^, 



9(0 



9(0 



T"''" [d,.^9{0) {dbX,){d,dax,) = 2 [djn^g{0 

where T"'"^'^ = gabgcd _|_ gacgbd _|_ gadgbc^ Making use of all that, introducing a new integration 
variable z"" = g^^'^C'^, and recovering the metric dependence by virtue of the equality 



(D'^DaX,) [D'D,xJ = g^' (dj,,) {dbt 



(41) 



we arrive after some straightforward algebra at the following expression for the string effective 
action up to the order 



5eff. =CJ j d\^+K I d\^n+^J d^^^g-\dat,,)idbt,,,) + 



\ 



, (42) 



where R is the size of the contour C in the confining regime (c/. Eq. (^). 

In Eq. (^), the coupling constants are completely determined via the zeroth and the first 
moments of the coefficient function D (as well as the omitted higher terms of the expansion are 
determined via the higher moments of the function D) as follows 



(43) 



5 J d'zz'D (z') 



and 



19 




Making use of Eq. (pS]), one can estimate them as 



Due to the lattice data quoted after Eq. (^), we thus get for the S'f/(2)-case the following values 
of the coupling constants a ~ 0.2 GeV^ (c/. the corresponding value before Eq. (^), k ~ 0.003, 
and — = —0.005. The obtained value of the Nambu-Goto term string tension demonstrates the 
agreement of the MFC with the present lattice results, while the two other obtained coupling 
constants are simply small numbers, which confirms the validity of the performed expansion. 

Before discussing various terms in Eq. (^2p , let us comment on the parameter of the performed 
expansion of the initial Eq. (|39|). First of all, as it has already been mentioned, this is an expansion 
in formal power series of Tg, which means that it is valid only when Tg is sufficiently small. 



That was a reason for the author of Ref. to call the limit when Tg is small but a is kept 
fixed as a "string limit of QCD". An n-th term of the expansion has the order of magnitude 

ttg (^(^F^ui^)^ R'^ , which means that the parameter of the expansion is {Tg/Kf. In the 

confining regime, this parameter is of the order of 0.04 (see the lattice data after Eq. (|5^)), i.e., 
is indeed a small number. Therefore in the "string limit", the operators of the lowest orders in 
the derivatives w.r.t. ^°'s dominate in the expansion of the full nonlocal action (^) 0. Contrary, 



in the QCD sum rule limit fT^, the effects brought about by the nonlocality of the functions D 
and Di are disregarded, and these functions are simply replaced by the gluonic condensate. In 
the language of the parametrization (pSj), this means that Tg — > oo, and our expansion diverges. 
This observation clarifies once more the relevance of the string picture of QCD to the description 
of confinement. We also see, that it is the MFC (where the correlation length of the vacuum can 
be considered as a variable parameter), which provides us with such a picture. 

Let us now proceed with the physical discussion of the obtained local string effective action 



The first term on the R.H.S. of this equation is the celebrated Nambu-Goto term with the positive 



string tension a. This term ensures the area law behaviour of the Wilson loop (^, since / d?^y g{^) 
is nothing else, but the area of the minimal surface, bounded by C. The obtained expression for 
the string tension coincides with the one obtained in Ref. [1^. Notice also that though the Nambu- 



Goto term alone is known to suffer from the problem of appearing of a tachyon in its spectrum, 
this problem is absent for the full nonlocal string effective action (|39D. 

The second term is known to be a full derivative in 2D, which can be most easily seen in the 
conformal gauge adopted, i.e., from Eq. (^Ol). This term is a topological invariant proportional to 
the Euler character of the world-sheet 

X = — y d'^^\fg{C}T^ = 2 — 2 X (number of handles) — (number of boundaries). 
In particular, for the surface Smin. under study, x = 1- 



-"^^Examples of the geometric structures, emerging in the expansion in the order Tg, are listed in the Appendix 
to the first paper of Ref. ||75| . 
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The most interesting term in the obtained effective action is the third one, which is usually 
referred to as the rigidity term |]80| , This term is not a full derivative, since it can be rewritten 
(just modulo full derivative terms) as 



an J 



i.e., the integrand does not contain the complete expression for ^/g^^. The rigidity term has been 
introduced into string theory in the above mentioned papers from the general arguments as the 
only possible one, which is invariant under the scale transformation Ax^(^). 

It is worth noting that during the derivation of the effective action (0), we have not used 
explicitly the fact that the surface Smin. = Smin. [C], described by the induced metric gatiO^ is 
the one of the minimal area for a given contour C. However, this surface possesses an important 
property distinguishing it from all the other surfaces bounded by C. Namely, it is defined by the 
(nonlinear) equation D"- D aX ^^.{E,) = together with the boundary condition a;^(^)|^ = x^{s). Due 
to this exceptional property of Smin. [C] and Eq. (^) , we see that at this surface the rigidity term 
vanishes. However, in a complete string picture of QCD including quantum fluctuations above 
confining background, there should appear an average over all world-sheets E's bounded by the 
contour C. It is natural to expect that the dominant contribution to this integral is brought about 
by Smin. [C] 7 i-G., Eq. ( p3D can be represented as 

(H^(Smi„.)) = Jdx^{0{W{^)). (45) 

Here, a;^(^) parametrizes the surface S, and the measure of the average over world-sheets, Dx^{^), 
should just be determined by the above mentioned quantum fluctuations. In writing down Eq. (^3D, 
we have assumed that the statistical weight of the average over world-sheets is given by the MFC- 
inspired Eq. (^) with the replacement Smin. — ^ S. If this assumption really holds, which should 
be justified by further investigations, then Eq. (^2]), with gab being the induced metric describing 
S, yields the first few terms of the local string effective action associated with the statistical weight 
(Vr(S)). Since S here is no more the surface of the minimal area, the rigidity term in this string 
effective action survives. 

As it has been for the first time mentioned in Ref. and then confirmed by related calcu- 
lations in Refs. ||82| and [^, the negative sign of the coupling constant ao of the rigidity term 



is an important property relevant to the stability of the string world-sheet. A simple argument 
in favour of this observation can be obtained by considering the propagator corresponding to the 
action (^2]), which for a certain Lorentz index A reads 

, 1 4 - (46) 

For negative a^, this integral yields 

{xxioxxm -- 



ln(/i|e|) + iro VI«oklel 



27ra 

with /i denoting the IR momentum cutoff and Kq standing for the modified Bessel function, while 
for positive ao an unphysical pole in the propagator occurs. Another arguments concerning the 
necessity of the negative sign of ao have been presented in the above mentioned papers. This 
property (called there "anti-rigidity" or "negative stiffness") has been demonstrated to ensure 
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vanishing of the imaginary part of the frequencies of small fluctuations of the world-sheet, which 
might lead to the instabilities of the latter one. Further, the requirement of the negative sign of 
ao has been employed in Ref . fS^ . There, a new nonlocal string action, which manifests negative 



stiffness, has been proposed as a good candidate for modelling the gluodynamics string. However 
contrary to our calculations, it remained unclear in that paper how the action introduced there can 
be derived either from the Wilson loop expansion or from some quantity relevant to gluodynamics 
[e.g. another vacuum amplitude or the gluodynamics Lagrangian). 

Once being integrated out, the small transversal fluctuations of the world-sheet, mentioned 
above, produce a renormalization of the string tension, which at the two-loop level has been 
calculated in Ref. and reads 



(Jr, 



a 



1 + «o- 



D 



IGtt 



1 + ln 



2{D 
+ «o — 



2){D - 1) / 4A 



2567r2 



In 



where A is an UV momentum cutoff, and D is the dimension of the space-time. It is worth noting 
that the appearance of the rigidity term yields also some modifications in the quark- ant iquark 
potential, string tension behaviour at finite temperature, and thermal deconfinement properties, 
all of which have been surveyed in Ref. 



— ^In , 



^ — and therefore becomes infinitely large due to the /i-dependence. 



An important problem associated with the action (^21) is that the corresponding propaga- 
tor (^) enjoys the limit of small separations, |^| -C —/==■ Indeed, in this limit it goes over into 

(xA(Oa:A(0)) 

This means that normals to the string world-sheet tend to be very short-ranged, i.e., the world- 
sheet becomes extremely crumpled. That is the reason, why this problem is usually referred to as 
the problem of crumpling of the string world- sheet. 

Let us now consider two possibilities of curing this problem. First of them has been put 
forward in Ref. for the case of the effective string theory emerging from compact QED in 
D space-time dimensions (the so-called confining string theory |8^, which will be considered in 
details in Section 4). There, it has been demonstrated that for the case D oo, the correlation 
function of two transversal fluctuations of the string world-sheet in this theory has an oscillatory 
behaviour at large distances. Such a behaviour indicates that the world-sheet is smooth rather 
than crumpled. One might expect that the same mechanism works in all gauge theories, whose 
confining phases admit a representation in the form of some effective string theory with a non-local 
interaction between the world-sheet elements. However, it is not obvious whether this mechanism 
can be extended to the non-Abelian case of gluodynamics, where in the nonlocal string effective 
action instead of the propagator of the massive vector field, appearing in the Abelian case, 
stands the coefficient function D{x). Though according to the lattice data the large 

distance asymptotic behaviour of the latter is indeed similar to the one of the massive vector field 
propagator, there nevertheless remain significant differences. 



Let us therefore turn ourselves to the second possibility, proposed in Ref. [§0| and elaborated 
on in Ref. |7^, which is more applicable to gluodynamics. It is based on the introduction of the 
so-called topological term, which is equal to the algebraic number of self-intersections of the string 
world-sheet, into the string effective action. Then, by adjusting the coupling constant of this 
term, one can eventually arrange the cancellation of contributions to the string partition function 
coming from highly crumpled surfaces, whose intersection numbers differ by one from each other. 

Thus it looks natural to address the problem of a derivation of the topological term from the 
gluodynamics Lagrangian. Such a term has been recently derived in Ref. [B^ for AD compact 
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QED with an additional ^-term. In the dual formulation of the Wilson loop in this theory (which 
is nothing else but the AD confining string theory mentioned above), the latter one occurred to be 
crucial for the formation of the topological string term. However, such a mechanism of generation 
of a topological term is difficult to work out in gluodynamics due to our inability to construct the 
exact dual formulation of the Wilson loop in this theory. Therefore, it looks suggestive to seek for 
some model of the gluodynamics vacuum, which might lead to the appearance of the topological 
term in the string representation of the Wilson loop in this theory. 

In Ref . , this idea has been realized by making use of recent results concerning the evaluation 



of the field strength correlators in the dilute instanton gas model ^0|- latter paper, it has 

been demonstrated that for the case of an instanton gas with broken CP-invariance, the bilocal 
field strength correlator contains a term proportional to the tensor Sp^pXp. This term is absent 
in the case of a CP-symmetric vacuum, since it is proportional to the topological charge of 
the system, V {1%^ — n^, where V is the four- volume of observation, and within the notations 
of Ref. n4 and n4 stand for the densities of instantons and antiinstantons (Ps and J's for 



brevity), respectively. Similarly, the paper |]TB[ also dealt with the approximation of a dilute I — I 
gas with fixed equal sizes p of Ps and J's. In the remaining part of this Subsection, we shall briefiy 
consider the main points of this paper. 



The new structure arising in the bilocal correlator in the / — / gas reads [^| 

A tr ((P^^(x, Xo)Pap(x', Xq))) = 8 K - n^) h ( ^^~^^^ ^ e^^^xp- (47) 

In Eq. (|47|) , the asymptotic behaviour of the function 1^ {z^) at z <^ 1 and 1 has the following 
form 



r?) ^ ^. (4. 



and 



2 



Ir(z')-^Prj\nz\ (49) 



\z 



respectively. 

In what follows, we are going to present the leading term in the derivative expansion of the 
correction to the nonlocal string effective action (pUj), which has the form 



AS-eff. = -lnA(iy(S^i,.)). (50) 

Here, A (Vr(Smin.)) is the corresponding correction to the expression (|33D for the Wilson loop, 
following from Eq. ( ^Tj) in the CP-broken vacuum. 

We shall not be interested in calculating corrections to the Nambu-Goto and rigidity terms 
arising due to additional contribution from the I — I gas to the function D{x). This can be 



easily done by carrying out the corresponding integrals (^) and ( ^4]) of the function D{x) in this 



gas. Notice only that, as it has already been mentioned in Ref. [0, due to the reasons discussed 



in details in Refs. |Q and [T^, a correction to the string tension (^31) , obtained in such a way 
from the I — I gas contribution to the function D{x), should be cancelled by the contributions 
coming from the higher cumulants in this gas. Clearly, this does not necessarily mean that the 



corresponding correction to the rigid string inverse bare coupling constant (^41) vanishes. Indeed 
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one can imagine himself a function D{x), for which / dtD(t) = 0, whereas / dttD(t) ^ 0, where 



t = z^. For example, this is true for the function D{t) defined as 

D{t) = < t < c; D{t) = -^,t> c, 

which can obviously be made continuous by smoothering the jump at t = c with some function, 
odd w.r.t. the line t = c. 

One can now expand the correction (|50|) , emerging from the term (^T]), in powers of the 
derivatives w.r.t. ^"'s in the same manner, as it has been done above for the nonlocal string 
effective action (|39D . Noting that an analogue of the Nambu-Goto term in this expansion vanishes 
since e^y\pt^ytxp = 0, we obtain 



i?2 



Here, 

is the algebraic number of self-intersections of the string world-sheet and 

(3 = 16V ("-4 - fli) J d^zz^Ir (z^) (51) 

is the corresponding coupling constant. 

Note that the averaged separation between the nearest neighbors in the I — I gas is given by 
L = {n^ + ^4)"*. According to phenomenological considerations one obtains for the SU{3)-case, 
p/L ^ 1/3 |91|] (see also Ref. |9^, where the ratio p/L has been obtained from direct lattice 
measurements to be 0.37 — 0.40). The parameter L should then serve as a distance cutoff in the 
integral standing on the R.H.S. of Eq. (plj). Taking this into account, we get from Eqs. (^8]), (^9]), 
and (|5l|) the following approximate value of /3 



/3~ (27rp)^(n4-n4) 




(52) 



Noting that the second term in square brackets on the R.H.S. of Eq. (|52D, emerging due to 
Eq. (^), is much larger than the first one, emerging due to Eq. (|48|), and making use of the value 
p-i ~ 0.6 GeV one obtains /? ~ 57680.4GeV-^ ■ {n^ -n^). 

In conclusion of this Subsection, we have found that in the / — / gas with a nonzero topological 
charge, there appears a topological term in the string representation of the Wilson loop. The 
coupling constant of this term is given by Eq. (|5TD. Together with the Nambu-Goto and rigidity 
terms (see Eq. (^2])), this term forms the effective Lagrangian of the gluodynamics string following 
from the MFC. 



2.2 Incorporation of Perturbative Corrections 

As it has been discussed in the beginning of the previous Subsection, the origin of the nonlo- 
cal string effective action (pOf ), which served as a starting point for the derivation of the local 
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action (^), is essentially nonperturbative. This is because the dominant contribution to the 
coefficient function D{x) is brought about by the strong background fields, which ensure con- 
finement. However, as it has been argued in Ref. |9^, in order to get the exponential growth of 
the multiplicity of states in the spectrum of the open bosonic string, one must account for the 
perturbative gluons, which interact with the string world-sheet. In the present Subsection, we 
shall proceed with studying this interaction by making use of the so-called perturbation theory in 
the nonperturbative QCD vacuum ||2^, |69| . 

To this end, we shall split the total field as = + a^, where i?^ ~ ^ is a strong 
nonperturbative background, and a°'s are perturbative fiuctuations around the latter, a° ~ qB"^- 
Thus, our strategy should be to perform an integration over a^'s in the expression ([T^ ) for the 
Wilson loop, where S is replaced by Smin.. However due to the path-ordering, which remained 
after rewriting the contour integral as a surface one in the version of the non-Abelian Stokes 
theorem proposed in Refs. 
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|, such an integration is difficult to carry out starting with 
Eq. (|17|). That is why, we find it convenient to adopt another version of this theorem, proposed 
in Ref. IQ, where the path-ordering is replaced by the integration over an auxiliary field from 
the SU{N,)/ ^" coset space. In what follows, we shall consider the S'f/(2)-case, where this 

field is a unit three-vector n, which characterizes the instant orientation in colour space, and the 
non-Abelian Stokes theorem takes a remarkably simple form 



{W{C)) 



Dnexp 



iJ 



+ 



abc^a 



(53) 



duA^^ + ge A^AI is a strength tensor of the gauge field, V 
is the covariant derivative, and J = ^,1,1,... is the colour 



Here, T"^^ 



d,5 



ab 



of representation of the SU{2) 



|, 1, |, . . . IS tne colour "spm 
group under consideration, defined via its generators T"'s as T"T" = J(J+ 1)1. The last term in 
the exponent on the R.H.S. of Eq. ( ^3]) is usually referred to as a gauged Wess-Zumino term. 

As we shall see below, such a version of the non-Abelian Stokes theorem will indeed enable 
us to carry out the one-loop integration over perturbative fiuctuations, which will then lead to 
a new {w.r.t. Eq. (^9t)) type of interaction between the string world-sheet elements. Once being 
expanded in powers of the derivatives w.r.t. ^"'s, this interaction will finally yield a correction to 
the rigidity term, keeping the Nambu-Goto term untouched. 

Let us start with the above mentioned splitting of the total gauge field Af^ in Eq. (|53D, which 



owing to the background field formalism 



§, 13 yields 



{W{C))=ArjDB;v{B'^)Dnexpyd'x^--{F;:,y+'- J da.^e^'^n^ {D.n)' {D^nY 

igj 



^min. 



X exp 



da^^n^F;^, X 



X 



' Da;: expj / d'x UofF'^, + la^l {d^^D^,, - 2igF;l 
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^min. 



X 



X 



DO" DO" exp 



(54) 



Here, and D'^ are the background field strength tensor and the corresponding covariant 



derivative, defined identically to J^^^ and with the replacement — > B'^. Secondly, in 
order to avoid double counting of fields during the integration and perform the averages over 
the background fields and quantum fiuctuations separately, we have used the so-called 't Hooft 
identity 



JDB"^viB'\JDalfiB, 



JDB-v{Bi) 

valid for an arbitrary functional /. Here, an integration weight rj (-Bq) should be fixed by the 
demand that all the cumulants and the string tension of the Nambu-Goto term acquire their 
observed values. Notice also that in a derivation of Eq. ( p4|) we have adopted the background 

Feynman gauge, £gaugefix. = -| (-DJ^^aj; 
It should be commented that the term 



-gaugefix. = -I (Dfal) , and denoted F;^ = F'^, {T^f with [T'^f'^ - ^-'^'^ 



—IE 



igj j d(T^,yal{D^nf 

^min. 



(55) 



in Eq. (|54D emerged from the expansion of the Wess-Zumino term as a result of the following 
sequence of transformations 



igJ 



igJ 



igJ / da^uCil (D^n) 



where in the last equality we have used the facts that = 1 and e^'^'^n^n'^ = 0. It is also worth 
mentioning that the terms quadratic in quantum fiuctuations, which emerge from the expansion 
of the field strength tensor JF^^^ and the Wess-Zumino term, cancel each other. 

In what follows, we shall work in the one-loop approximation, and thus disregard in Eq. ( ^4]) the 
terms cubic and quartic in quantum fiuctuations, as well as the ghost term. For simplicity, we shall 
also neglect the interaction of two perturbative gluons with the field strength tensor F^^ (gluon 
spin term). Notice that within the Feynman- Schwinger proper time path- integral representation 
for the perturbative gluon propagator, which will be used immediately below, such a term leads 
to insertions of the colour magnetic moment into the contour of integration (see Ref. [p^). 

Bringing now together the term (^3|) and the term —igJ J da^j^n'^D'^afl from Eq. (|5^ and 

performing Gaussian integration over perturbative fiuctuations, we obtain 



{W{C)) 
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exp 



igj 



da^^^rfF'^^ \ exp 



n Bj; 



(56) 



Here, 



{...)^ = J Dn{...)exp J da^^e'^^^n'^ (D^n)' (D^n) 

(. . jDB;{...)^ (5^) exp -\ J d'x (f;,) 

and the (generally speaking, non-translation-invariant) interaction kernel K^p{x,x') is expressed 
via the perturbative gluon propagator as follows 



K^p{Xj X ) 



+00 



dxpdx'p 



n^{x)n^{x') 



ds I Dze 



P exp I ig J dXzaBa 



be 



(57) 



where z{0) = x' , z{s) = x. In the derivation of Eq. (pBf), we have neglected the interaction of 
the string world-sheet with the background sources of the type D1^F^^{y), where y is an arbitrary 
space-time point outside the world-sheet, which should be finally integrated over. 

In order to derive a correction, emerging due to exchanges by perturbative gluons, to the 
background- induced gluodynamics string effective action (|39|) , let us apply to Eq. ( ^6l) the following 



formula 65 



n=l 



where A and B stand for two commuting operators, and (...) is an arbitrary average. Then, the 
leading correction, we are interested in, corresponds to the complete neglection of correlations 
between the arguments of the first and second exponential factors standing on the R.H.S. of 
Eq. (^). Secondly, it corresponds to putting the n- and -B|JJ- averages of the second exponential 
factor inside it. Taking all this into account and following our definition of the string effective 
action Ses_. as — In (W{C)), we obtain 



S^s, = - In ^ ^exp 



igJ 



da n"- F"- 
2 J "t^^i/"- pu 

^min. 




A5, 



cff. 



(5J 



n B<^. 



where the desired leading correction to the string effective action reads 

ASeff. = ^ / da,,{x) I dap,{x'){{Kpp{x,x'))^)^^. (59) 
V ■ y ■ ^ 

Clearly, in the bilocal approximation, the first term on the R.H.S. of Eq. (^) yields pure back- 
ground part of the string effective action ( ^9]) (c/. Eq. (|53D with replaced by 5^). 
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The obtained correction ( p^D to the pure background string effective action corresponds to 
a new type of interaction between the string world-sheet elements. Namely, instead of the propa- 
gator of the background gluon, represented by the function D{x), in Eq. (p9D stands a propagator 
of the perturbative gluon in the nonperturbative gluodynamics vacuum. Due to the statistical 

weight e P exp ig j dXzaBa of this gluon, it is the region where s is small, which mainly 



contributes to the interaction kernel (|57|). This means that the dominant contribution to the 
obtained correction (|59| ) stems from those points x^(^) and of the world-sheet, which are 

very close to each other. That is in line with the performed derivative expansion of the nonlocal 

<Tg<^R. Taking this into account, we can adopt the 



string effective action (|39D , where 

simplest, local, approximation for the propagator {n"[x]n"'[x' 



i.e., replace it by 



^bc 

T 



j Dnexp 



iJ 

Y 



This expression is a functional of the world-sheet as a whole (z.e., it is independent of x^{^)) and 
therefore can be absorbed into the irrelevant normalization constant M . 

Finally, in order to perform an expansion of the nonlocal correction (|59|) in powers of the 
derivatives w.r.t. ,^"'s and derive from it the first few local terms, it is convenient to pass to the 
integration over the trajectories u^{\) = z^{\) + ^{x' — x)^ — x'^. This enables one to extract 
explicitly the dependence on the points x^ and from the integral over trajectories (necessary 
for the differentiation w.r.t. these points), which yields 



dx^dx'p „ 



dse 



Due X 



X tr P exp 




-|- u] Ba (u + x' + —(x — x') 



with u{0) = u{s) = 0. Then, in the bilocal approximation, the dominant contribution to the 
Nambu-Goto and rigidity terms comes about from taking the derivatives of the free propagation 

factor e 4s and replacing the parallel transporter factor by the one over the closed path. 



which has the form trP exp 



. Finally, substituting the so-obtained expression 



for Kfj,p{x,x') into Eq. (|59| ) and performing an expansion of this nonlocal correction in powers of 
the derivatives w.r.t. ^"'s similarly to the previous Subsection, we arrive at the desired correction 
to the local effective action (|42|) . In this way, it turns out that the Nambu-Goto term string tension 
does not acquire any correction due to perturbative gluonic exchanges, whereas the correction to 
the inverse bare coupling constant of the rigidity term reads 



1 



dss 



Due ( tr P exp 



ig 



dXuaBJu) 



(60) 



^2 



Note that since we have proved in the previous Subsection that the rigidity term vanishes at the 
surface of the minimal area, Smin., the obtained result means that accounting for perturbative 



28 



gluons in the lowest order of perturbation theory does not change the background-induced string 
effective action (|4^ ) associated with this surface. However, the derived correction (|60|) is essential 
for the action ( ^2]) corresponding to the statistical weight (VF(E)) in the integral over world-sheets 
standing on the R.H.S. of Eq. (^S]). 

It is also worth noting that since ~ ^, the parallel transporter factor on the R.H.S. of 
Eq. (pOf ) cannot be expanded in powers of g and should be considered as a whole. The path- 
integral on the R.H.S. of Eq. (|60|) is not simply the perturbative gluon propagator, since the 
integral over the proper time contains an additional power of s, which makes the whole quantity 
dimensionless, as it should be. Notice also that the sign of the obtained correction ( |60|) to the 
inverse bare coupling constant (0) depends on the form of background 5^ entering the Wilson 
loop (which is the most nontrivial content of this correction) on the R.H.S. of Eq. (|60D . However, 
since the derived correction is a pure perturbative effect (due to the factor g'^ present in Eq. (pDD), 
even in the case when it is positive, it cannot change the negative sign of the leading term (^). 



2.3 A Hamiltonian of the Straight-Line QCD String with Spinless 
Quarks 

In the present Subsection, we shall derive a Hamiltonian corresponding to the quark- ant iquark 
Green function ([l5| ) in the confining QCD vacuum. To this end, we shall first write it down in 



the Feynman-Schwinger proper time path-integral representation (p!6|), after which substitute for 
{W{C)) the above obtained expression exp(— 5*05.) with S^s. defined by Eq. (^2]) (and inverse 
bare coupling constant of the rigidity term (Q) modified by Eq. (pO|), i.e., ^ + A^)- It 



is worth noting that for the first time such a Hamiltonian has been derived in Ref. |96], where 



however, only the Nambu-Goto term in the string effective action (|42|) has been accounted for. 
Our aim below will be the derivation of a correction to this result due to the rigidity term, which 
we assume to be nonvanishing. That means that the QCD string in the problem under study is 
excited, since its world-sheet differs from the one of the minimal area. As a byproduct, we shall also 



rederive the leading terms in the Hamiltonian of the QCD string with quarks obtained in Ref. ^6 
However as we shall see, the rigid string theory, being the theory with higher derivatives, leads to 
the interesting and phenomenologically relevant modifications of the Nambu-Goto-induced part 
of the Hamiltonian. We shall also generalize the result of Ref. by considering the case of 
different masses of a quark and ant iquark. 



The two main approximations, under which we shall consider the Green function (IT3), are the 



same as the ones used in Ref. First, we shall neglect quark trajectories with backward motion 



in the proper time, which might lead to creation of additional quark-antiquark pairs. Secondly, 
we shall use the straight-line approximation for the string world-sheet S[C], which, as it has been 
argued in Ref. , corresponds to the valence quark approximation. Such a string may rotate and 
oscillate longitudinally. This approximation is inspired by two limiting cases, / = and / — 00, 
where / is the orbital quantum number of the system. The first case will be investigated below in 
more details, and the correction to the Hamiltonian of the relativistic quark model due to the 
rigidity term in the limit of large masses of a quark and antiquark will be derived. 

Let us now proceed with a derivation of the desired Hamiltonian. To this end, we shall start 
with the expression for the Green function (p!6| ) with {W{C)) defined via Eq. (^) and 
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K = mls + - z^dX, K 



mis 



- / z^dX 
4 



standing for the kinematical factors of a quark and antiquark, respectively. Then, by making use 
of the auxihary field formalism , one can represent it in the following form 



+00 

G{x,x;y,y) = J dT J DzDzD^iD^2Dhabew ~ exp {—a + 2a) J d'^^Vh 



X 



X exp 



-a 



C(q J 



(61) 



where we have integrated over the Lagrange multiplier A*^ (^) = (0 + / (0?/ "-afc — 

and d is the mean value of a{C,). Here t^j^iy = {daWfj) {dbW^), 



K' + K' = ^ I dr 



/ rfr ^ + ^i(r) (1 + z\t)) + ^ + /.2(r) (l + ^^(r 
J /ii(r) ^ ^ u^ir) ^ 



T 



[xo + xo-yo-yo), /ii(r) 



T 



and the no-backtracking time approximation ||9^ 



(62) 



/ii(r) > 0, /i2(r) > 



(63) 



has been used. Similarly to Ref . |Q , we employ in the valence quark sector ( |63D the approximation 
that the string world-sheet S[C] can be parametrized by the straight lines, connecting points 
z^{t) and z^{t) with the same r, i.e., the trajectories of a quark and antiquark are synchronized, 
z^ = (r, z),z^ = (r, z), w^(r, /?) = Pz^ir) + (1 - /5)^^(r), 0<P<1. 

Introducing auxiliary fields [^] i^{t,P) = Ta^,r]{T, P) = -fj^ and making a rescaling — > 

§^z^, one gets from the last exponent on the R.H.S. of Eq. ( |6l]) the following action 



2a" ^A" 



of the string without quarks 



( 2 







aQO^ h 



•■2 2 /•■ \2 •2*2 /••\2 

w r — (wr) + w r — [wr) + 



+2{{ww){rr) — {wr){wr)) + 



a 



I r^r^ — ffr)^ 



— 2r] ({wr)r'^ — {wr){fr) + {wf){rr) — (wr)r' 



(64) 



where a dot stands for and r^(r) = z^(r) — ^^(r) is the relative coordinate. 

Let us now introduce the centre of masses coordinate -R^(t) = ({t)z^{t) + (1 — ({t))z^{t), 
where ({r) = Ci(t) + — C2(''"),0 < ((t) < 1, should be determined from the requirement that 
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decouples from r^. These extremal values of (i and (2 can be obtained from the corresponding 
saddle-point equation. Referring the reader for the details to Appendix A, we shall present here 
only the final result for the path-integral Hamiltonian of the straight-line QCD string with quarks. 
It has the form 



H = H^^^ + —H^^l (65) 



Here 



= - 
2 



{pr ^ + ml) {pr ^ + ml) 2 2 } d3 L 



yUi + /i2 + 0- r / h z/Q H 







(66) 



with 



P = /ii + ^2 - ^^^^'l^' , = /ci/?/?V, pI ^ ^P^, L ^ [r,p] 
pi + /i2 + z/Q ^ 

is the Hamiltonian of the Nambu-Goto string with quarks, which for the case of equal masses of 
a quark and an antiquark has been derived and investigated in Ref. |9^. The Hamiltonian H^^\ 
which is a new one, reads 

=7^w ^^^^ + ;^ ^^^^ + + w ^'^^ 

where p = ^f^^f^, and the coefficients a^'s, A; = 1, . . . , 5, are listed in Appendix A. 

We see that the obtained Hamiltonian (^) contains not only corrections to the orbital momen- 
tum of the system, but also several operators higher than of the second order in the momentum. 
The latter ones emerge as a consequence of the fact that the rigid string theory is a theory with 
higher derivatives. 

The obtained (path-integral) Hamiltonian (|65|)-(|67|) contains auxiliary fields pi, p2 and v. In 
order to construct out of it the operator Hamiltonian, which acts upon the wave functions, one 
should integrate these fields out. This implies the substitution of their extremal values, which can 
be obtained from the corresponding saddle-point equations, into Eqs. (|65|) -(|67|) and performing 
the Weil ordering |98|| . 

In conclusion of this Subsection, let us apply Hamiltonian (|67|) to a derivation of the rigid 
string correction to the Hamiltonian of the so-called relativistic quark model ||9^, i.e., consider 
the case when the orbital momentum is equal to zero. Let us also set for simplicity the mass of a 
quark equal to the mass of an antiquark mi = m2 = m. In order to get H^^^ , one should substitute 
the extremal values of the fields pi, p2, and z/ of the zeroth order in fif^'^^- = ^i^^^- = y/p'^ + 
and t'extr. = cr\r\, into Eq. (0). The limit of large masses of a quark and antiquark means that 
m ^ y/a. In this case, we obtain from Eq. (^) the rigid string Hamiltonian H'--^'' = — ^sJli (p^)^ 
and then from Eqs. (B^) and (pBf) the following expression for the total Hamiltonian H 



if = 2m + a|r| + P^-(— - + ^)(p^)-. (6J 



m 




One can see that the new rigid string-inspired term, quartic in the relative momentum of the 
quark-antiquark pair, may cause a sufficient infiuence to the dynamics of the system. 
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3 String Representation of Abelian-Projected Theories 



3.1 The Method of AbeUan Projections 



The MFC, exploited in the previous Section, enabled us to investigate the string effective action 
associated with a certain {e.g. minimal) world-sheet, but did not yield a prescription of a derivation 
of the full string partition function in the form of an integral over all world-sheets (|45|) . The 
reason for that is that within the MFC one looses the meaning of the path-integral average over 
the gluodynamics vacuum, replacing this average by the phenomenological one. As a consequence, 
it looks difficult to extract singularities corresponding to QCD strings out of the resulting vacuum 
correlation functions. 

In the present Section, we shall proceed with a derivation of such an integral over world- 
sheets in the QCD-inspired effective Abelian-projected theories and find an exact field-theoretical 
analogue of the phenomenological background gluon coefficient function D{x). We shall argue 
that the general features of the gluodynamics vacuum are in line with the ones predicted by the 
so-called dual Meissner picture of confinement, first put forward by 't Hooft and Mandelstam [^ |. 
According to the 't Hooft-Mandelstam scenario, the properties of the QCD string should be similar 
to those of the electric vortex, which emerges between two electrically charged particles immersed 
into a superconducting medium filled with a condensate of Cooper pairs of magnetic monopoles. In 
the case of the usual Abelian Higgs Model (AHM), which is a relativistic version of the Ginzburg- 
Landau theory of superconductivity, such vortices are referred to as Abrikosov-Nielsen-Olesen 
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Up to now, there does not exist any analytical proof of the existence of the condensate of 
Abelian monopoles in QCD, though in lattice QCD there are a lot of numerical data in favour 
of this conjecture (see e.g. Refs. [^9|, |100| , |101| , |102|| ). However in all theories allowing for an 
analytical description of confinement, the latter one is due to monopoles, which either form a gas 
or condense. These examples include compact QED and 3D Georgi-Glashow model ||103| , p!9[] , and 
supersymmetric Yang- Mills theories [^. As far as the origin of magnetic monopoles in QCD is 
concerned, they appear in the so-called Abelian projection method [^, |104|| , which we shall briefly 
describe below. 

The essence of this method is based on a certain partial gauge fixing procedure, which re- 
duces the original gauge group SU{Nc) to the maximal Abelian (or Cartan) subgroup [U{1)]^''~^ , 
i.e., leaves Abelian degrees of freedom unfixed. Then, since the original SU{Nc) group is com- 
pact, the resulting Abelian subgroup is compact as well, which is just the origin of Abelian 
magnetic monopoles in the original non-Abelian theory. To perform such a gauge fixing, one 
chooses a certain composite operator X, transforming by the adjoint representation of SU{Nc), 
X —* X' = UXW, and diagonalizes it, i.e., finds such a gauge (unitary matrix U), that 
X' = diag (Ai, . . . , Aat^). As an example of such an operator may serve e.g. Fi2{x) or, at fi- 
nite temperatures, the so-called Polyakov line |1(J5 



P exp 



dxdA. 



where /? = 1/T is an inverse temperature. Notice that the most interesting numerical results (see 
e.g. [^]) have been obtained for the S'f/(2)-case in the so-called Maximal Abelian gauge, defined 
by the condition of maximization of the functional 
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R[A,] = - J d^x [{Alf + {Aiy 



which reads max R 

u 



A'. 



. Here, A'^ is the transformed vector potential following from the original 
one after the gauge transformation with the matrix U. The condition of the local extremum for 
this functional, (j)^ ± igA'^^ A'^ = 0, where A'^ = A'^ ± iA'^, means that in this gauge we make 
the field A^ as diagonal as possible. 

The matrix U, which diagonalizes the operator X, is obviously defined up to a left multiplica- 
tion by the diagonal SU (Nc) matrix, which belongs to the subgroup [U{l)f^'\ This is just the 
way how the original non-Abelian group reduces to the maximal Abelian subgroup. If we now 
perform a certain Abelian projection, i.e., consider the transformed vector potential 

a; = u(^A^+ '-d,^ U\ (69) 

then it is straightforward to see that under the remained maximal Abelian subgroup, the diagonal 
field components of the matrix-valued vector potential = \AA ,2 = 1,... ,Nc, transform as 

Abelian gauge fields, a^^+^d^ai, whereas the off-diagonal components = {A'^ transform 
as charged matter fields, exp [i (oj — a^)] c^. The latter ones can be further disregarded in 

a derivation of an effective IR Abelian-projected theory owing to the so-called Abelian dominance 
hypothesis ||106| |. As far as the monopoles {w.r.t. the "photon" fields of diagonal elements of the 



vector potential) are concerned, their role is played by the singularities of the matrix U, which 
might occur if some of the eigenvalues Ai, . . . , Aat^ of the operator X coincide. This may happen 
in some 3D point xq, which becomes a world-line of magnetic monopole in 4D. 

Let us start with illustrating the above described ideas at the simplest example of the SU (2)- 
gluodynamics. The initial Yang-Mills action reads 



5'ym 



4* 



-tr / d'xF^, (70) 



with F^^ = Fj^^TK Here, F;^ = d^,Al - d^A'^ + ge'^^A^^A^, = ^,i = l, 2, 3, and r*'s stand for 
Pauli matrices. 

One can perform the gauge transformation (^) so that the gauge-transformed field A'^ obeys 
e.g. the above-mentioned maximal Abelian gauge fixing condition (d^±iga' )A'if = 0, where 



a'^ = A^ 1^. The gauge transformed field strength tensor then reads = U yF^^, + F^Jf^'j 

where the singular contribution has the form -F^"^' = ^ (^[d^,d,y] U^^ U. This contribution comes 
about from the generally singular character of the matrix U of the gauge transformation, mentioned 
above, and describes world-sheets of the Dirac strings ||107|| . Clearly, integration over all possible 



singular gauge transformations results to an integration over F^'J'^ 's. 

Let us next single out the diagonal (neutral) component = of the field A^ by making 
use of the decomposition 



A, = a^T' + A;^T'' = A^ + C^, 



Notice that the maximal Abehan gauge fixing condition can be written as follows = 0, where = 

d^5°'^ ~ ge°'^^a'^. Once being rewritten in this form, this gauge can be easily recognized as the standard background 
gauge with the field a' playing the role of the background. 
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where a = 1,2. Consequently, one has for the field strength tensor 



F^, [A + C] = F^, [A] + {D [A] A C)^^ - ig [C^, C 



(71) 



where 



(O A G)^^ = O^G, - ae^M' and [A] =d^-ig [A^, 
Eq. (fflD can be straightforwardly rewritten as follows 



Here, /^j, = (9 A a)^^ and C^,^ = ge"'^^ A'^^A^^ stand for the contributions of diagonal and off-diagonal 
components of the gluon field to the diagonal part of the field strength tensor, respectively, and 



Iq^ab /\ jg ^YiQ off-diagonal part of the field strength tensor with V 

\ / fJ,U 



Ca 

This yields the following decomposition of the action 
fields, 



ab 



dj'^^-ge 



a63 , 



Ul) taken now on the gauge transformed 



5- 



YM 



4 ' 



smg. 



2 1 



smg. 



a\ 2 



?smg. 



2tr (T^F^jfs-). 



where (^F^ 

Next, since our aim will be the investigation of the confining [i.e., infrared) properties of the 
Abelian-projected SU (2)-gluodynamics (rather than the problems of its renormalization, related 
to the region of asymptotic freedom), we shall disregard the A^-dependent terms 
this approximation, the resulting effective action takes the form 



1061. Within 



Q „ n f sing 
•Jeff. (J-fi^ J fMu 



where we have denoted for brevity /^"^' 
modified Bianchi identities as 



\ J d'x {f,. + fP^-)\ (72) 
(-^^"^') • The monopole current is defined via the 



■M 



with = \e^y\pf\p. Thus the obtained effective theory ( [7^ ) can be regarded as a U{1) gauge 
theory with magnetic monopoles. To proceed with the investigation of the monopole ensemble, 
it is reasonable to cast the partition function under study, Z = j Df^^^^-Dape'^"^-, to the dual 
form This can be done by making use of the first-order formalism, i.e., linearizing the square 
/^^ in Eq. (^2|) by introducing an integration over an auxiliary antisymmetric tensor field Bp^, as 
follows 



2= / Df;^>^-DapDBp,exp 



d X 



smg. 
pu 



(73) 



A"^'" ' 2^ 

^^Notice that the gauge fixing term for the Abehan field is assumed to be included into the integration measure 
Da„. 
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Integration over the a^-field leads to the constraint [B 



I ysing. 



0, whose resolution yields 

^iiv = A B)^^, where B^ is now the "magnetic" potential dual to the "electric" potential 



a^j_. Substituting this representation for B^,, into Eq. ([73|), we obtain 



Z = (^JdB, exp |- J d'x [\f'^^ - })^^ 



(74) 



where from now on F^^, denotes {d A B)^^. In Eq. (|7^) , the integration over /^™^''s has transformed 
to a certain average over the monopole currents, (. . ■)jAi, whose concrete form will be specified 
immediately below. 

Clearly, the next crucial step necessary to proceed with the summation over the ensemble of 
monopoles is to postulate the properties of their ensemble. In a derivation of effective Abelian- 
projected theories, it is usually assumed that those form a condensate of their Cooper pairs of the 
charge 2gjn, where stands for the magnetic coupling constant 1^. This assumption is realized 
firstly by setting for the collective current of magnetic Cooper pairs the expression 



N 
n=l ■' 



where the world-line of the n-th Cooper pair is parametrized by the vector a;^(s), and secondly 
by specifying the measure (. . .) -m to the following form |p.08| , |109| , ^ 



exp {i j d'^xB^fJ 



N=l ■ 



n / — e 

n=l ''q 



u{0)=u{s„) 



X 



X exp 



N 

E 

'=1 



/ ds'i (--ii^is'i) + 2ig.raU^,is'l)B^{u{s'l))] - 4A ^ ds'J d46 [u{s'i) - u{4 

\ / l,k = lQ 



(75) 



Here, the vector m^(s^) parametrizes the same contour as the vector x'^{s). Clearly, the world-line 
action standing in the exponent on the R.H.S. of Eq. (|7|) contains besides the usual free part 
also a short-range interaction term, which after carrying out the path-integral yields the Higgs 
potential of magnetic Cooper pairs. Indeed, Eq. (^) can be rewritten as an integral over the 
effective magnetic Higgs field describing Cooper pairs as follows 



exp [i j dSB^j^y^ = J D^D^* exp | - y d^x 



2 ' ^ ' 



(76) 



where = — 2ig, 



mBfj, is the covariant derivative, and an inessential constant factor has been 



referred to the integration measure []. Finally, substituting Eq. into Eq. (|7|), we arrive at 
the following effective Abelian-projected theory of the S'f/(2)-gluodynamics 

Within this assumption, we restrict ourselves to the sector of the full yet unknown Abelian-projected theory 
where there do not exist antimonopoles. Possible interference between these two sectors requires further investiga- 
tions. 

"'^^A seeming divergency at large proper times produced in Eq. ( |75| ) by the factor e^^'' is actually apparent, 
since the last term in the exponent on the R.H.S. of this equation yields the desired damping. 



35 



Z = J \<^\D\<^\DeDB^expl^- J rf^x [^F^, + ^ + A - r^^)'] | , (77) 

where $(x) = |$(x)| e'^^^-*. Clearly, in Eq. (|77D one can recognize the partition function of the 
dual AHM (DAHM). 

Next, analogous considerations can be applied to the S'f/(3)-gluodynamics. Demanding the 
condensation of magnetic monopoles, one can expect that in this case the resulting Abelian- 
projected theory should be also of the DAHM type, but with the [[/(l)]^ gauge invariance. To 
derive it, let us start with the Cartan decomposition of the action ( [TOl) in the S'f/(3)-case. Then, 
= d^,Ai - d^A'^ + gP^A^Al, and the ^[/(3)-generators = f , z = 1, . . . , 8, obey the 
commutation relations [T^T-'] = if^jkj'k^ where A*'s denote the Gell-Mann matrices. The Cartan 
decomposition of the field A^, which singles out its diagonal part, has the form 



A^ = a^H + (c*;Ea + C^F_„) =A, + C, 

a = 1,2,3, where = (^A'^^,A^^, and the diagonal S'f/(3)-generators, which generate the Cartan 
subalgebra, read H = {Hi,H2) = {T^,T^). We have also introduced the so-called step operators 
E±aS (else called raising operators for positive a's and lowering operators otherwise) by redefining 
the rest (non-diagonal) 5'f/(3)-generators as follows 

Clearly, these operators are non-Hermitean in the sense that {Ea)'^ = E^a- The complete Lie 
algebra of the so-redefined S't/(3)-generators reads 

[H, E±a\ = ±eaE±a, [E±a, E±f,] = T-^^abcE^c, [Ea, E^b] = SabGaH- 

Here, we have introduced the so-called root vectors 

ei = (l,0), e2=(4'-l^)' ^3=(-^,^), 

which thus play the role of the structural constants in the first of the above commutation relations. 
Eq. (|7lD then still holds, and we obtain for various terms on its R.H.S. the following expressions 

F,M] = f...H, f^. = idA a)^, , iD[A] A C)^, = K (P- A C^)^, + (V^ A C")^, , 

-ig [C„ C] = zg |e,H (C'^ A C^,, + ^Satc [q^Cl^E^, - ClClE, 

where = df^ + igeaSL^. Bringing all these expressions together, we arrive at the following Cartan 
decomposition of the field strength tensor 

F^u = (f^!^ + Cfj_u) H + S*^Ea + S^^j^E^a, 

where C^, = igea A C*")^, and S'^, = (D" A ^^ + ^eabcC*J'C:'^. Next, it is worth employing 
the S'f/(3)-version of the maximal Abelian gauge, which for every index a reads 2^^C° = 0. 
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Analogously to the SU (2)-case, this leads to the appearance of the singular contributions to the 
field strength tensor. Accounting for them as well as for the equations trHi^^-ta = tr E±aE±b = 
and tr EaE^i, = ^6ab, we finally obtain the following Cartan decomposition of the SU{3) Yang-Mills 
action 



YM 



'■fiu 



P I f sing. 



where f^s- = 2tr (HF^jfs-) and (F^jf^-)" = 2tr (E^F^jf*?-). Again, disregarding the off-diagonal 
part of the action owing to the Abelian dominance hypothesis, we are left with the 5'f/(3)-analogue 
of the effective action (^), which reads 



5, 



eff. 



- j d'x (f,. + f; 



sing. 



Next, by the dualization of the obtained theory we arrive at the following partition function (c/. 
Eq. @) 



Z = (^j DB^expj- j d' 



X 



(79) 



where F^^, = {d A B)^^ is the field strength tensor of the field B^, which is dual to the field 



andjf = 9^fps-. 



At this point, it is worth mentioning that the root vectors define the lattice at which monopole 
charges are distributed. Taking this into account, it is straightforward to write down the expression 
for the collective current of magnetic Cooper pairs 



n=l a=l 

The average over magnetic currents then has the form 



N 3 . 
29m E E / dxl^^^is)S {X - . 



exp (^i J d'^xB^ 



Ju a=l 



OO 1 

N=l ■ 



N 



+ 00 



n=l 



ds. 



n — e 



M(a)(0) = M('^)(s„) 



X exp 



N 



i=ii ^ ^ 



-4AE jds\jdsl5 [u^''\s[)-u^'^{sl) 





Lk=l 



j D$,D$:exp |-| d^a; E - 2^^7me„B^) $,|' + A - 



X 



(80) 



where the vector u^^\s'^ on the L.H.S. parametrizes the same contour as the vector x^'^'*"(s). 
Finally, it is worth noting that since monopoles are distributed over the root lattice, whose vectors 



are related to each other by the condition = 0, the dual Higgs fields $a's are also not 

a=l 
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completely independent of each other. In Ref. ||104]|, it was argued that the relevant constraint 



for these fields reads Yl, da = 0. Taking this into account we arrive at the following partition 

a=l 

function describing an effective [f/(l)]^ gauge invariant Abelian-projected theory of the SU{3)- 
gluodynamics ||104 



X exp 



d X 



Z 



4 i^" 



De,DB^6 




3 

+ E 

a=l 



where $a 

Notice that in general for the gauge group SU (Nc) the number of independent fields of distinct 
monopoles in the corresponding effective Abelian-projected theory is equal to A^^c — 1. This number 
is just equal to the number of independent closed Abrikosov-Nielsen-Olesen type strings in this 
theory. 

In the present Section, we shall investigate string representations of the models (IT 



and (|TD 

and their extensions due to the introduction of external electrically charged particles, which we 
shall call "quarks" , as well as various correlators in these models. Various topological properties 
of Abelian-projected theories will also be considered. Besides that, we shall briefly discuss the 
phenomenon of chiral symmetry breaking from the point of view of these theories. 



3.2 Nonperturbative Field Correlators and String Representation of 
the Dual Abelian Higgs Model in the London Limit 

3.2.1 String Representation for the Partition Function of Extended DAHM 

Our first aim in the present Subsection will be the derivation of a string representation for the 
partition function of DAHM, extended by introduction of external electrically charged {w.r.t. the 
maximal Abelian U{1) subgroup of the original SU{2) group) particles, which we shall refer to as 
"quarks" . Such an extension can be performed by adding to the action (^) the term i J d'^xa^jj^ 

with j^(x) = e § dx^{s)6{x — x{s)) standing for the conserved electric current of a quark, which 

c 



moves along the closed contour C {cf. the notations to Eq. (|T^)). The electric coupling constant e 
is related to the magnetic one according to the topological quantization condition (which replaces 
the standard Dirac quantization condition in the case of QCD-inspired theories) egm = 4:7m, 
where n is an integer. In what follows, we shall for concreteness restrict ourselves to the case of 
monopoles possessing the minimal charge, i.e., set n = 1. 

Then, performing the dualization (|72D-([74D with the new term included and summing up over 
monopole currents according to Eq. ([75|), we arrive at Eq. ( |77| ) with F^^ replaced by F^^ + F^. 
Here, stands for the field strength tensor generated by quarks according to the equation 

= j^. A solution to this equation reads = — eS^^, where T,^^{x) = J da^,y{x{^))6{x — 
x(^)) is the so-called vorticity tensor current [|1 1 1|| defined on an arbitrary surface T,^ (which is 



just the world-sheet of an open dual Abrikosov-Nielsen-Olesen string), bounded by the contour C. 
Due to the Stokes theorem, the vorticity tensor current is related to the quark current according to 
the equation edvT^^y = j^. In particular, this equation means that in the case when there are no 
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external quarks, the vorticity tensor current is conserved, i.e., due to the conservation of electric 
flux all the strings in that case are closed. Notice that when external quarks are introduced into 
the system, some amount of closed strings might nevertheless survive. 

Let us now proceed with the string representation for the partition function of the theory ([TTD, 
extended by external quarks, in the so-called London limit, A — >^ oo. In this limit, the radial part 
of the magnetic Higgs field can be integrated out, becoming fixed to its v.e.v., |$| — > t], and the 
partition function under study takes the form 



Z = / DB^De''''^ De"'^- exp 



+ F, 



flU 



(82) 



where from now on constant normalization factors will be omitted. In our further interpretation of 
the topic of the string representation for the partition function and field correlators of the extended 



DAHM, we shall mainly follow Refs. ||112|| and ||113|| (for related investigations see Refs. |8^, |109| , 

In Eq. (|82D, we have performed a decomposition of the phase of the magnetic {i.e., dual) Higgs 
field 6 = 6^^^^- + ^''^s-. Here, 6^^^^'{x) describes a certain configuration of electric strings and obeys 
see e. 



the equation 



lie 



sing. 



27rS^j,(x) 



^3) 



This equation is just the covariant formulation of the 4D analogue of the Stokes theorem for the 
field c}p^^"S'(x), written in the local form. Here, S^^, stands for the vorticity tensor current, defined 
at the world-sheet S of a closed electric string, parametrized by the vector x^(^). On the other 
hand, the field 6^'^^'{x) describes simply a single- valued fluctuation around the above-mentioned 
string configuration. Note that as it has been argued in Ref. | |115| |, the integration measure over 
the field 6 factorizes into the product of measures over the fields 6^^^^- and 6'^'^^'. 

Performing the path-integral duality transformation of Eq. (^) along the lines described in 



Ref. ||116|| , we get 



Z = Dxp{^)Dh^^exp< 



d^x 



12t]' 



y 



where T.^u = 4S„„ — and H, 



-'flU 



'Xll'fiv 



+ d,yhxn is the field strength tensor of a massive 



antisymmetric tensor field /i^jy (the so-called Kalb-Ramond field [|117|] ). This antisymmetric spin-1 
tensor field emerged via some constraints from the integration over O'^^^' and describes a massive 
dual gauge boson. As far as the integration over the world-sheets of closed strings, Dx^(^), is 
concerned, it appeared from the integration over 6^^^^' by virtue of Eq. ([5B|), owing to which there 
exists a one-to-one correspondence between 0^^^^- and a;^(^). Physially this correspondence stems 
from the fact that the singularity of the phase of the Higgs field just takes place at the string world- 
sheets (Notice that since in what follows we shall be interested in effective actions rather than 
the integration measures, the Jacobian emerging during the change of the integration variables 



X, 



i(^), which has been evaluated in Ref. [ 115 



will not be discussed below and is assumed 
The details of a derivation of Eq. ( plD are outlined in 



to be included into the measure Dx^(^).) 
Appendix B. Thus, the path-integral duality transformation is just a way of getting a coupling of 
the dual gauge boson, described now by the field /i^;^, to a string world-sheet, rather than to a 
world-line (as it takes place in the usual case of the Wilson loop). 
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Finally, the Gaussian integration over the field h^^p in Eq. (^) (see Appendix C) leads to the 



following expression for the partition function (|8^ 



exp 



e2 



j dx^j dy^D^^{x -y) 
c c 

X j L)x^(Oexp -{nrif j d'^x j dSf:^^{x)Dl^\x - y)f:^^{y) 



^5) 



Here, D^^(x) = ^J^^^^ Ki{m\x\) is the propagator of the dual vector boson, whose mass m is equal 
to 2gmf], and -ft'j's, i = 0,1, 2, 3, henceforth stand for the modified Bessel functions. 

Note that since quarks and antiquarks were from the very beginning considered as classical 
particles, the external contour C explicitly enters the final result. Would one consider them on 
the quantum level, Eq. (^) must be supplied by a certain prescription of the summation over 
the contours. This can be done by linearizing the current x current interaction standing in the 
first exponential factor on the R.H.S. of Eq. (pSf) by integration over an auxiliary massive electric 



vector field interacting with the current and further averaging over j^'s by virtue of some 
formula d la Eq. (|TSp. 

Clearly, the above mentioned exponential factor is the standard result, which can be obtained 
without accounting for the Abrikosov-Nielsen-Olesen type electric strings. Contrary to that, the 
integral over string world-sheets on the R.H.S. of this equation stems just from the contribution 
of strings to the partition function and is therefore the essence of the desired string representation 
(c/. Eq. (^)). The respective string effective action describes both the interaction of the closed 
world-sheets S's with the open world-sheets E^'s and self-interactions of these objects. 

In particular, comparing the obtained nonlocal string effective action, standing in the second 
exponential factor on the R.H.S. of Eq. (|85|) , with Eq. (|39|) we obtain by virtue of Eqs. (^31) 
and (0) the following values of the string tension of the Nambu-Goto term and the inverse bare 
coupling constant of the rigidity term for the surface E: 

a = 2717]^ Ko{c) ~ 2717]^ In - (86) 

c 

and ^ = — Y^2-- Obviously, these two quantities, corresponding to the world-sheet S^, are in the 
factor 16 larger. Here, c stands for a characteristic small dimensionless parameter. In the London 
limit, this parameter is of the order of the ratio of m to the mass of magnetic Cooper pair (which 
plays the role of the UV momentum cutoff somehow analogous to the inverse lattice spacing 1/a 
{cf. Eq. (^)), i.e., c ~ Moreover, in a derivation of Eq. ( PB| ) we have assumed that not only 
- ^ 1, but also In^ ^ 1, i.e., the obtained expression for the string tension is valid with the 
logarithmic accuracy. Note that the logarithmic divergency of the string tension in the 3D version 
of AHM, the Ginzburg-Landau model, is a well known result, which can be obtained directly 
from the definition of this quantity as a free energy per unit length of the Abrikosov vortex (see 
e.g. ||119|] ). The physical origin of this result is that in the centre of the vortex, the condensate 
is destroyed, and dual gauge bosons remain massless. Possible phenomenological consequences of 



this effect were recently discussed in [120 



Clearly, both the string tension and the inverse bare coupling constant of the rigidity term 
are nonanalytic in Qm, which means that these quantities are essentially nonperturbative simi- 
larly to the QCD case {of. Eq. ([121) ). Notice also that the finite temperature behaviour of the 
magnetic Higgs field v.e.v. ||119[| , ri{T) oc y^l — ^ with Tc standing for the critical temperature. 
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obviously governs the corresponding behaviour of the string tension ( P^ ) and, in particular, the 
deconfinement phase transition. 

Finally, it is worth noting that since contrary to the previous Section, the obtained nonlocal 
string effective action is no more associated with the world-sheet of the minimal area, Smin. [C]; 
the rigidity term generally does not vanish. 



3.2.2 String Representation of the Bilocal Field Strength Correlator 

The partition function of extended DAHM, investigated in the previous Subsection, can obviously 
be applied to the derivation of the string representation for the bilocal correlator of the field 
strength tensors. Indeed, owing to the Stokes theorem, this partition function can be written as 



Cm J 



M ' 



where (. . .)„ 



psmg. 



with Scs. 



and (. . .)jM given by Eqs. (^) and (|7S]), respectively. Applying to this expression the cumulant 
expansion, we have in the bilocal approximation: 



Z 



exp 



^7) 



Following the MFC, let us parametrize the bilocal cumulant by the two Lorentz structures similarly 
to Eq. (H): 



+ 



vAx'). 



Owing to Eqs. ( pS] ) and p^), Eq. ( pH] ) eventually yields 



Z^expl-Jd'xJ d'y ^Sj,(a:)Sj,(y)P {{x - yf) + (x)jj(y) | d\V,{X) 



{x-yY 

On the other hand, this expression should coincide with Eq. (pS]) divided by Z 
i.e., it reads 

Z = exp {- / d'x J dSDt\x - y) [(4vrr^)^Sj,(x)Sj,(y) + (x)jj(y)^ 
8(7rr/)2| d'x J d'yD^^\x-y)J:^,ix)J:M 



X 



89) 



X ( exp 



where 



(...) 



/ Dx.iO (. . .) exp \-i7rrjf J d'x J d'yE,,{x)D(^\x - y)T.^,{y) 



J Dx^iO exp -{nr^y J d^x J d^yT.,,{x)D^A\x - y)T.,M 
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This is just the standard normahzation condition, encoded in the measure Da;^(^). 



(90) 
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In the confining regime of tlie Wilson loop describing external quarks, we are interested witli, tlie 
area of the surface is much larger than the typical area of the surface E, swept out by closed 
electric strings. Owing to this, the average over world-sheets standing on the R.H.S. of Eq. (|90D 
can be disregarded w.r.t. the first exponential factor. Then, the comparison of the latter one with 
Eq. (|89[) straightforwardly yields for the function V the following expression 



m? Ki{m\x\) 
Atc'^ \x\ 



(91) 



whereas for the function Vi we get the equation / d\T)i(\) = 4D^^(x), which leads to: 

Ki{m\x\) m 



Vi {x 



m 



+ — I KQ{m\x\) + K2{m\x\) 



(92) 



We see that in the limit |x| ^ ^, the asymptotic behaviours of the coefficient functions (|9TD 
and (P^) are given by 



V 



m 



-m\x\ 



4^/2712 (^rn \x\)^ 



(93) 



and 



m 



-m\x\ 



(94) 



For bookkeeping purposes, let us also list here the asymptotic behaviours of these functions in the 
opposite case, \x\ <^ —. Those read 



V 



m 



Att'^x'^ 



(95) 



and 



2 1 |4 ■ 



(96) 



One can now see that according to the lattice data [^, |62|, the asymptotic behaviours (|93|) 
and (|9^) are very similar to the large distance ones of the nonperturbative parts of the functions 
D and Di [cf. Eq. (^)). In particular, both functions decrease exponentially, and the function V 
is much larger than the function T>i due to the preexponential power-like behaviour. We also see 
that the dual gauge boson mass m corresponds to the inverse correlation length of the vacuum 
T~^. In particular, in the string limit of QCD, when Tg —>■ while the value of the string tension is 

kept fixed, m corresponds to \f^^- Clearly, the found similarity in the large-distance asymptotic 
behaviours of the bifocal cumulant of the field strength tensors in DAHM and the gauge-invariant 
cumulant in QCD supports the original conjecture by 't Hooft and Mandelstam concerning the 
dual Meissner nature of confinement. 

Moreover, the short distance asymptotic behaviours ( P^D and (^) also parallel the results 
obtained within the MFC of QCD in the lowest order of perturbation theory. Namely, at such 
distances the function Di to the lowest order also behaves Eq. (|37|) ) and is much 

larger than the function D to the same order. However, it should be realized that the effects of 
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asymptotic freedom (which is the most important UV feature of QCD, distinguishing it from all 
the Abelian gauge theories) obviously cannot be obtained from DAHM. 



Hence we see that within the approximation |S| ^ 



relevant to the confinement of external 



quarks, the bilocal approximation to the MFC is an exact result in the London limit of DAHM, 
i.e., all the cumulants of the orders higher than the second one vanish. Higher cumulants are 
naturally calculable upon performing the average over world-sheets on the R.H.S. of Eq. (|90D 
(In particular, this average yields also the important modification of the bilocal cumulant, which 
will be discussed below, in Subsection 4.3.). Clearly, higher cumulants appeared in this way are 
dependent on the properties of the ensemble of closed strings, one averages over. In particular, 
this means that the relations between cumulants of various orders do depend on these properties 
as well. However, it turns out that such relations may be established even without performing 
the average over world-sheets. This occurs to be possible in the London limit of DAHM with an 
additional term, which describes an interaction of an axion with two dual gauge bosons [|121|| . The 
reader is referred to that paper for a detailed discussion of the relation between the bilocal and 
threelocal cumulants in such an extended model. 



3.2.3 Dynamical Chiral Symmetry Breaking within the Method of Field Correlators 

The QCD Lagrangian with Nj massless flavours is known to be invariant under the global symme- 
try transformations, which are the U{Nf) x U{Nf) independent rotations of left- and right-handed 
quark fields. This symmetry is referred to as chiral symmetry. The above mentioned rotations 
of the two-component Weyl spinors are equivalent to the independent vector and axial U{Nf) 
rotations of the full four-component Dirac spinors, under which the QCD Lagrangian remains 
invariant as well. At the same time, the axial transformations mix states with different P-parities. 
Therefore, we conclude that if the chiral symmetry remains unbroken, one would observe parity 
degeneracy of all the states, whose other quantum numbers are the same. The observed splittings 
between such states occur, however, to be too large to be explained by the small bare or current 
quark masses. Namely, this splitting is of the order of hundreds MeV, whereas the current masses 
of light u- and d-quarks are of the order of a few MeV 0. This observation tells us that the chiral 
symmetry of the QCD Lagrangian is broken down spontaneously. This phenomenon of the spon- 
taneous chiral symmetry breaking (SCSB) naturally leads to the appearance of light pseudoscalar 
Goldstone bosons, whose role is played by pions, which are indeed the lightest of all the hadrons. 
Besides confinement, the explanation of SCSB is known to be the other most fundamental problem 
of the modern theory of strong interactions. 

The order parameter of SCSB is the so-called quark condensate (else called chiral condensate) 
~ —(250 MeV)^. This is nothing else, but the quark Green function taken at the origin or a 
closed quark loop in the momentum representation. Clearly, if the quarks from the beginning are 
treated as massless {i.e., only the kinetic term is present in their propagators), the phenomenon of 
SCSB leads to the appearance of a non-zero dynamical quark mass, which in general is momentum- 
dependent. At zero momentum, the value of the dynamical quark mass is of the order of 350-400 
MeV, which is just the value of the so-called constituent quark mass. 

Let us now proceed to the quantitative description of SCSB. Integrating the quarks out of the 
QCD Lagrangian obviously yields 

^^The current mass of the s-quark, which is around 150 MeV, is still smaller than the typical splitting values at 
least in a factor of three. 
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/=i " 

where f)^^'^'^- = j^D^^'^'^' with D^^^'^' = — igA^ standing for the covariant derivative in the 
fundamental representation. As it follows from its definition, the chiral condensate for a given 
flavour / then has the form 

where V stands for the four-volume of observation. In what follows, let us for simplicity restrict 
ourselves to the case of one flavour only. Next, let be an eigenfunction of the Dirac operator, 
corresponding to a nonvanishing eigenvalue zD^^'^'^'n = A„\E'„. Then, since 75 anticommutes 
with iD^^^'^', the function \E'„' = j^'^n is also an eigenfunction of the Dirac operator, corresponding 
to the eigenvalue A„' = — A„. Owing to this observation, the determinant can be rewritten as 
follows 



det (it) 



fund. 



+ im 



exp 



1 
2 



exp 



m 



+00 



dXu{X) In (A^ + 



m 



where i/(A) = X)^(A — A„) stands for the so-called spectral density of the Dirac operator. By 

n 

virtue of Eq. (p7|), we thus obtain for the chiral condensate 

+00 



1 



dX{u{X)) 



m 



A^ + 



(9J 



m—>0 



where (z^(A)) denotes the spectral density averaged over the full QCD partition function, includ- 
ing the weight given by the determinant itself. The latter can be disregarded in the quenched 
approximation (justified at large- Ai'c), in which one neglects the backward influence of quarks to 
the dynamics. 

For a flnite-volume system, the R.H.S. of Eq. (p8|) obviously vanishes. However, for the case 
when the volume increases, the spectrum becomes continuous, and we should use the formula 



m 

lim -— 



sign (m)7r5(A). 



Taking this into account, one flnally gets the celebrated Banks-Casher relation ||122 | 

1 



V 



sign (m)7r {iy{0)) ■ 



(99) 



Thus, we see that the chiral condensate is proportional to the averaged spectral density of the 
Dirac operator at the origin. Notice that the sign function in Eq. (^) means that at small m, the 
QCD partition function depends on m nonanalytically, which is typical for the situation when the 
symmetry is spontaneously broken. 
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Up to now, there exist several microscopic models of SCSB in QCD |T^. The most elaborated 
and therefore popular of them is the one due to instantons p2| It is based on the observa- 
tion that in the background field of one (anti)instanton, the Dirac operator has an exact zero 
mode ||124|| . The estimate of (z^(0)) for the case of a gas of J's and J's has been performed in 
Ref. |22| directly by calculating the overlap of the zero modes and has the form (z^(0)) ~ 0. 
This result then has been rederived in Ref. ||125|| by solving a closed equation for the averaged 
quark propagator as an expansion in powers of the so-called packing fraction parameter 
where N stands for the total number of /'s and J's. By virtue of Eq. (pUj), the chiral condensate 
reads (ijip'^ 

Notice that this result can be also derived if we first average over the / — / gas and only after 
that calculate the chiral condensate by making use of the so-obtained effective theory In that 
case, quark interactions are due to the scattering of two or more (anti) quarks over the same I or 
/, i.e., there arise four- (or more) fermion interaction terms. The range of such interaction, which 



is usually referred to as 't Hooft interaction |p.24||, is obviously of the order of p. In the case of 



two flavours, this interaction is the four-fermion one and yields an effective theory similar to the 
well known Nambu-Jona-Lasinio model [^, 

There have been also derived other types of nonlocal NJL models in an approximate way from 
QCD by making use of path-integral techniques in bilocal fields. In particular, exact bosonization 
of 2D QCD has been performed in Ref. |p.26|| , and several attempts to bosonization of 4D QCD have 
been proposed ||127| , The general strategy of these approaches includes several steps. Firstly, 
one casts the effective four-fermion interaction (which can be either local or nonlocal with a certain 
interaction kernel) into the form of a Yukawa interaction by introducing a set of collective bosonic 
fields, which is just the essence of the bosonization procedure. Secondly, one can integrate over 
the fermions and derive an effective action in terms of these collective fields. Next, for large- A'^c the 
saddle-point of this effective action yields the so-called Schwinger-Dyson (or gap) equation, which 
determines the dynamical mass, responsible for SCSB. After that, an expansion of the resulting 
meson action in small field fluctuations around this stationary solution yields the so-called Bethe- 
Salpeter equation, which determines the spectrum of meson excitations. This finally leads to a 
description of the chiral sector of QCD in terms of Effective Chiral Hadron Lagrangians pi], 153, 
|37| containing higher order derivative terms with fixed structure constants. The form of these 
Lagrangians agrees with that of the phenomenological Hadron Lagrangians postulated at the 
end of the Sixties on the basis of pure group-theoretical arguments when considering nonlinear 
realizations of chiral symmetry |TB[. However, contrary to those Lagrangians, the new ones are 



not obtained by symmetry principles alone, but rather derived from an underlying microscopic 
quark (diquark) picture, which enables one to estimate masses and coupling constants of composite 
hadrons. 

However, it is worth noting that though the QCD-motivated NJL type models model well the 
SCSB phenomenon in QCD and describe with a good accuracy the hadron spectrum and coupling 
constants, they do not reproduce the confinement property. Recently, an attempt of a derivation 
of an effective Lagrangian for a light quark propagating in the confining QCD vacuum, which 



could account simultaneously for both phenomena, has been done |47|. Within this approach, one 



Recently in Ref. |123| , there has been proposed a model based on the so-called dyonic gas. The advantage 
of this model w.r.t. the instanton one is that, as it has been argued there, it yields not only SCSB, but also the 
confinement property. 

Here, we adopt the same notations as the ones at the end of Subsection 2.1. However, now p and L stand not 
for fixed, but for the averaged size and separation in the I — I ensemble, respectively. 
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starts with an expression for the Green function of a system consisting of an infinitely heavy quark 
and hght antiquark. The propagation of a hght (anti)quark alone in the QCD vacuum then leads 
to SCSB, whereas the presence of a heavy quark enables one to take into account the QCD string, 
joining both objects, and thus to incorporate confinement. Averaging over gluonic fields by making 
use of the cumulant expansion in the bilocal approximation, it is then straightforward to derive 
an effective NJL-type Lagrangian for the light quark with a certain nonlocal kernel. However, 
the form of this kernel is now determined via the coefficient function D{x) {cf. Eq. (0)), which, 
as it has been demonstrated above, yields the string picture {cf. Eqs. (|^), (^), (PD), i.e., is 
responsible for confinement. Such an interpolation between confinement and SCSB within this 
approach leads to a relation between chiral and gluonic condensates, which reads 

{^^) cc ~D{0)Tg. (100) 
Following Ref. let us briefiy demonstrate how this relation can be obtained in the Abelian- 



projected theories. To this end, we for simplicity consider the case of Abelian-projected SU{2)- 
gluodynamics, studied above. (Our analysis can be straightforwardly extended to the Abelian- 
projected S't/(3)-gluodynamics, which will be investigated in the next Subsection.) To proceed 
with, let us assume for a while that a monopole has an effective infinitesimal size b (which is 
known to become finite for the 't Hooft-Polyakov monopole). Then, every monopole with the 



world-line of the length T produces T/b quasizero modes [|128|| . The total number of quasizero 
modes produced by all the monopoles from the 3D volume V^^^ in the interval of modes AA thus 
has the form (z^(A)) AA = jV^^^n^, where is the 3D density of monopoles. Since V = TV^^\ 
one gets an estimate 

MO)) r^3 
-^^AXb- 

Since the quasizero modes of all the monopoles are mixed due to the interaction, the denominator 
of the R.H.S. of this relation is of the order of unity. On the other hand, the 3D density of 



monopoles can be estimated via the following equation 



which follows from Eq. (|^) due to equations of motion. Indeed, owing to Eq. (|102|) , one has 



J d'x {j^ (x, X4) (0, X4)) ~ V{0)Tg. (103) 

Here, as it has been argued after Eq. (^), the correlation length of the vacuum, Tg, for the Abelian- 
projected S'f/(2)-gluodynamics is equal to the inverse mass of the dual gauge boson, m~^. Notice 
also that in a derivation of the first estimate in Eq. ( |103| ), one uses the following observation |^ . 
The correlator (j/i(a;)j^(0)) estimates the probability of finding a monopole at the 4D point x, if 
there is one at the origin. Integrating over d^x, one finds the probability of having a monopole at 
the origin, while another one is anywhere. Fixing 0:4 means that the probability refers to a given 
moment. It has been also assumed that one magnetic monopole yields one quasizero fermion mode 
per unit length of its world-line. This is true for an isolated monopole, and this result has been 
extrapolated to the QCD vacuum as a whole. 

Finally, substituting Eq. ( p.03|) into Eq. ( |101j ) and making use of Eq. (l99|), we arrive at Eq. (|100[) . 
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3.3 Nonperturbative Field Correlators and String Representation of 
the 5't/(3)-QCD within the Abelian Projection Method 

In the present Subsection, we shall extend the results of Subsection 3.2 to the case of the SU{3)- 



QCD. Our interpretation of this subject will mainly follow Ref. |p.29| , |130|| . We shall start with 
the string representation for the partition function (^T|) in the London limit. 



3.3.1 String Representation for the Partition Function of the Abelian-Projected 
SU (3)-Gluodynamics 

Similarly to the S'f/(2)-case, in the London limit, the monopole fields become infinitely heavy, 
and their radial parts can be integrated out. After that, we are left with the following partition 
function 



Z 



X exp< / d'^x 



H' O, (1 



4 /^^ 



„2 3 

y E id A 



2g. 



\a=l 



a=l 



(104) 



Similarly to DAHM, in the model (|8TD, there exist string-like singularities of the Abrikosov- 
Nielsen-Olesen type. That is why, in Eq. ( |104|) we have again decomposed the total phases of the 
monopole fields into singular and regular parts, 9a = O^'^^' + G'^a^', and imposed the constraint of 
vanishing of the sum of regular parts by introducing the integration over the Lagrange multiplier 
k{x). Analogously to the DAHM, in the model ( |104D , 6'^™^ 's describe a given electric string 
configuration and are related to the world-sheets S^'s of strings of three types via the equations 



smg. 



x) = 27rS)J^(a;) = 27r y d(j^i,{xa{0)S{x - Xa{0), 



(105) 



where Xa = a;^(^) is a four- vector, which parametrizes the world-sheet Sq. 

The path-integral duality transformation of the partition function ( [L04| ) parallels that of Sub- 
section 3.2.1. The only seeming problem brought about by the additional integration over the 
Lagrange multiplier occurs to be trivial due to the explicit form of the root vectors. Indeed, let 
us first cast Eq. ( |104| ) into the following form 



Z 




X 



^ {c;y + ^c; {d,ea - 2^?„e.B,) + E C^- 



a=l 



(106) 



and carry out the integration over In this way, one needs to solve the equation 5^C^ = k, 

which should hold for an arbitrary index a. The solution to this equation reads 



k{y) 
{x - yf 
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where h^p stands for the Kalb-Ramond field of the a-th type. Next, making use of the constraint 

3 

J2 ^a™^' = 0, replacing then the integrals over 6'^™S''s with the integrals over a;"(^)'s by virtue of 



a=l 



Eq. ( |105| ), and discarding again for simplicity the Jacobians ||115|| emerging during such changes 
of the integration variables, we arrive at the following representation for the partition function 



X 



X 



X 



Here, -ff^j.^ = ^tJ.KyX + ^x^liu + ^y^Xn stands for the field strength tensor of the Kalb-Ramond field 

3 

. Clearly, due to the explicit form of the root vectors, the sum eaB^ vanishes, and the 



integration over the Lagrange multiplier thus yields an inessential determinant factor. Notice also 

3 



that due to Eq. (105), the constraint OT^^' = resulted into a constraint for the world-sheets 

a=l 

3 

of strings of three types Yl, ^t.u ~ 0- This means that actually only the world-sheets of two types 

a=l 

are independent of each other, whereas the third one is unambiguously fixed by the demand that 
the above constraint holds. 

Straightforward integrations over the dual gauge field as well as over the Kalb-Ramond 
fields then yield the following desired string representation for the partition function 



2 = J Dxi.ios 

\a=l / 



X exp 
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2 J da^^iXaiO) J da^^ixaiO) 



Ki {niB \Xa{C) - Xa{i')\) 



\Xa{i) ~ Xa{i')\ 



(107) 



Here, = y/GgmV is the mass of the fields B^, which they acquire due to the Higgs mechanism. 
Finally, it is possible to integrate out one of the three world-sheets, for concreteness a;^(^). This 
yields the expression for the partition function in terms of the integral over two independent string 
world-sheets 



Z = j DxliODxliOx 



X exp < -Qm-n 
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+ / do- 



S2 



(105 



According to Eq. ( 108 ), in the language of the effective string theory, the partition function (|104|) 
has the form of two independent string world-sheets, which (self-) interact by the exchanges of 
massive dual gauge bosons. Notice also that as it follows from Eq. (|108|) , the energy density 
corresponding to the obtained effective nonlocal string Lagrangian increases not only with the 
distance between two points lying on the same world-sheet, but also with the distance between two 
different world-sheets, which means that also the ensemble of strings as a whole displays confining 
properties. This observation will be elaborated on in the next Section, where the collective effects 
in the string ensemble will be studied. 



3.3.2 String Representation of Abelian-Projected S'?7(3)-QCD and the Bilocal Field 
Strength Correlator 

An external quark of a certain colour c = R, B,G (red, blue, green, respectively) can be introduced 
into Abelian-projected 5'f/(3)-gluodynamics by adding to the initial action (|78D the interaction 
term iQ^^^ J d^xa^j^. Here, jp{x) = § dx^{s)5{x — x{s)), and the vectors of colour charges read 

The electric coupling constant g here is again related to the magnetic one, Qm, via the topological 
quantization condition gg-m = 4:7rk with k standing for an integer, henceforth set to unity (the 
generalization to an arbitrary k is straightforward). Then, applying the Stokes theorem and 
the cumulant expansion in the bilocal approximation, we get for the partition function of the 
Abelian-projected 5'[/(3)-QCD with an external quark of the colour c the following expression: 






exp 



{c)if){c)j 



d'x d'yE,^{x)E,,{y)((f;,{x)fl{y) 



(109) 



Here, i,j = 1,2 are the [?7(l)]^-indices, referring to the Cartan generators H, and S^,^ is the 
vorticity tensor current defined at a certain world-sheet S, bounded by the contour C. The average 
on the R.H.S. of Eq. ([lOD is defined as follows: (. . .)^^ j„ = (/ Da^ exp (-^eff. [a^, f^^'] ) (. . ■)).^^ 

with S'eflf. and (...):« given by Eqs. (^) and (|80|), respectively. Upon the parametrization of the 
bilocal cumulant as 
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we can write for Eq. ( p.09| ) the following expression 
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where it has been used that for every c, (Q*-'^-')^ — 3 

On the other hand, the partition function can be calculated exactly. Indeed, the dualization 
of the action ( [78D with the term iQ'^'^^ / (i'^xa^j^ added, leads to Eq. ( [79D with F^,^ replaced by 
F^j^+F^fJ. Here, F|fj stands for the field strength tensor of a test quark of the colour c, which obeys 
the equation d^^'^^^ = Q^''''j„ and thus can be written as F^^^^^ = — Q'^'^^S^j,. Then, the summation 
over magnetic currents in the sense of Eq. ( pOD yields Eq. ( pTD with the same extension of F^^. In 
the London limit under study, the path-integral duality transformation of this action and further 
integration over the Kalb-Ramond fields (see Ref. | JL3(]| | for details) yield 



z, = J Dx;{os {t.^^^u 
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(112) 



Here, = S^^, — 2s^^T.^^ with the following numbers s^^'s: s^^'s: s 
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bi — 63 
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53*^^ = 1, which obey the relation Q*^'^^ = ^eaS^^\ Taking 

into account that for every c, (si'^-*) = 2, we eventually arrive at the following expression for the 
partition function (c/. Eq. (pOD): 
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with the average over world-sheets defined as 



(...) 



!Dxi{05[Y.^l 



exp 



{nl^Y J d^x J d^y^,{x)D(^l {x - y)S^,(y) 



a=l 



JDx-iOSiE S« )exp -{n7^yjd^xJd^yJ:i,{x)Di^l{x~y)^,iy) 



Comparing now Eq. ( |111| ) with Eq. ( |113| ), we see that in the approximation ^ valid in 
the confining regime for an external quark, the fuctions D and Di are given by Eqs. (|9T|) and (|92D 
with the replacement m ^ itib- Besides that, it is obvious that the bilocal cumulant (|110|) is 
nonvanishing only for the gluonic field strength tensors of the same kind, i.e., for i = j = loii = 
j = 2. Hence, for these diagonal cumulants (whose large-distance asymptotic behaviours match 
those of the S'?7(3)-gluodynamics as well as Eqs. (|93|) and (0) correspond to the S'f/(2)-case), 
the vacuum of the Abelian-projected S'f/(3)-QCD in the London limit does exhibit a nontrivial 
correlation length T„ = —. 

Finally, it is also worth noting that there has also been derived a string representation of 
the Abelian-projected S'f/(3)-QCD, extended by the introduction of the 6-term. The reader is 
referred to Ref. [|13CI|| for details of this investigation. 
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3.4 Representation of Abelian-Projected Theories in Terms of the 
Monopole Currents 

In the present Subsection, we shall derive a representation for the partition functions of the 
Abelian-projected theories directly in terms of the monopole currents. In our interpretation we 
shall follow Ref. imi. 



Let us again start our analysis with the 5'f/(2)-case, i.e., DAHM and for simplicity consider the 
unextended case, where intuitively the resulting monopole effective action should contain besides 
the free part, quadratic in the monopole currents, also the interaction of these currents with the 
closed string world-sheets E's. As our starting point will serve Eqs. (B.2) and (B.3) of Appendix B, 
where, however, we shall not perform any hypergauge fixing for the field /i^,^. Then, the partition 
function has the form 



12r]^ 

(114) 

Notice that according to the equation of motion for the field A^, the absence of external electric 
currents is expressed by the equation dfj,J-'^u = 0, where J-'fj,iy = dfj^A^, — d^Afj, + gmhf_iu- Regarding 
JF^jy as a full electromagnetic field strength tensor, one can write for it the corresponding Bianchi 
identity modified by the monopoles, d^^f^^, = gmd^h^^,. This identity means that the monopole 
current can be written in terms of the Kalb-Ramond field /i^^^ as 

ifi = gmdjiu^i, (115) 

which manifests its conservation. (Notice that this current is related to the field from Eq. (B.l) 
of Appendix B as = —gmCfj,. This means that the (5-function in the last equality on the R.H.S. 
of this equation just imposes once more the conservation of the current j^.) 

It is also instructive to write down the equation of motion for the Kalb-Ramond field in 
terms of the introduced full electromagnetic field strength tensor. This equation has the form 
•^i/A = ^d^Hf^i^x + 2^S;^A- By virtue of conservation of the vorticity tensor current for the closed 
string world-sheets, d^T^^y = 0, this equation again yields the condition of absence of external 
electric currents, d^T^u = 0. 

Let us now turn ourselves to a derivation of the monopole current representation for the 
partition function of DAHM. To this end, we shall first resolve the equation ^e^u\pdi,hxp = —j^ 
w.r.t. hfj_y, which yields 



Next, we get the following expressions for various terms on the R.H.S. of Eq. (|1 14|) 

Bringing all this together and performing in Eq. (|114| ) again the hypergauge transformation h^^ 
hfj,u + dfj^Xu — dyXp with the gauge function = — -^A^, which eliminates the field A^, we finally 
arrive at the desired monopole current representation, which has the form 
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5m 



;ii6) 



The first term in the exponent on the R.H.S. of Eq. ( |116D has the form of the Biot-Savart energy 
of the electric field generated by monopole currents ||109|| , the second term corresponds to the 
(gauged) kinetic energy of Cooper pairs, and the term 



47r2 



X 



y-x\ 



(117) 



describes the interaction of the string world-sheet with the monopole current j^. Obviously, this 
interaction can be rewritten in the form S'int. = / d'^xjp^H^^^' , where H^^^' is the four-dimensional 
analogue of the magnetic induction, produced by the electric string according to the equation 



Notice that if we approximate the current by the classical expression, i.e., set 



:ii8) 



j^{x) = 2gm f dy^{T)5{x -yir)), 



(119) 



the interaction term (|117|) takes the form S'int. = 2(7mL(S,r), where L(S,r) is simply the Gauss 
linking number of the world-sheet S with the trajectory of the monopole Cooper pair, F 0. 

It is straightforward to extend the above analysis to the case of the Abelian-projected SU{3)- 
gluodynamics, where the string representation for the partition function ( |104[ ) has the form 

Z = J Dx'^iOS [e Da;Dh% exp j- | d'x (h;,,)' - 



Ka=l 



(120) 



with = e^a^. Analogously to the argumentation following after Eq. (|114|) , Eq. (|120|) means 
that the arising monopole currents can be expressed in terms of three Kalb-Ramond fields as 

ja 

a=l ^ 



Qmxl ^dvhl^ 0. Finally, rewriting Eq. ( |120| ) via these currents and resolving the constraint 
by integrating over one of the world-sheets (for concreteness, again a;^(0)) W6 obtain 



Z = j DxliODxliODh%exp< 



27r2 



d'xd'yj^i 



x] 



[X 



y) 



m 



B 



^^Topological interactions of this kind are sometimes interpreted as a 4D analogue of the Aharonov-Bohm effect. 
In particular, this interaction, albeit for the current of an external electrically charged particle with the world-sheet 
of magnetic string, emerges in the string representation for the Wilson loop of this particle in AHM |115[. 
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Clearly, these currents are related to the current , which enters Eq. ([79|), as — e^j 
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+ 27ri 




:i2ii 



We see that the last four terms on the R.H.S. of Eq. (|121|) describe an interference between 
various possibihties of the interaction between string world-sheets and monopole currents in the 
Abelian-projected S'?7(3)-gluodynamics to occur. 

For illustrations, let us establish a correspondence of the above results to the 3D ones. Namely, 
let us derive a 3D analogue of Eq. (|116|) , i.e., find a representation in terms of the monopole 
currents of the dual Ginzburg-Landau model. There, Eq. ( p3D is replaced by [ p.09|| 



27r5^(x) 



;i22) 
)) being 



Here, on the R.H.S. stands the so-called vortex density with 5^(x) = / (ix^(r)(5(x — x(r 

L 

the 5-function defined w.r.t. the electric vortex line L, parametrized by the vector x(r). This 
line is closed in the case under study, i.e., in the absence of external quarks, which means that 
d^S^ = 0. Performing by virtue of Eq. (|122|) the path-integral duality transformation of the 
partition function of 3D DAHM in the London limit. 



we get for it the following representation 



(123) 



Z = J D^{T)DhD^expi^- J d' 



X 



1 



2 {df,hu 



duh. 



2TTihf,6f, 



(124) 



Notice that the Kalb-Ramond field has now reduced to a massive one-form field h with the mass 
m = 2gmi], as well as the A^-field has reduced to a scalar ip. Analogously to the 4D case, the field 
E = gmV2h + Vv9 can be regarded as a full electric field, defined via the full dual electromagnetic 
field strength tensor as = \e^y\J-',j\. The absence of external quarks is now expressed by the 
equation VE = 0, following from the equation of motion for the field (p. Correspondingly, the 
monopole currents are defined as jv = d^J^^i, = gmV^^fn/xd^hx and are manifestly conserved. 
Notice also that the condition of closeness of the vortex lines, d^6^ = 0, unambiguously exhibits 
itself as a condition of absence of external quarks, VE = 0, by virtue of the equation of motion 



for the field h, which can be written in the form E,. 



1^2 



Next, after performing the gauge transformation h ^ h + V7 with the gauge function 7 
-^V^ the field ip drops out. Expressing via as follows 



_d_ 

dx,. 



d'y 



jxiy) 



■ ■ ■ Ay/2TTgjn ' Ox^ J |x - y| 

and substituting this expression into the R.H.S. of Eq. (|124|) , we finally arrive at the desired 
representation for the partition function of 3D DAHM in terms of the monopole currents 
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Z = j Dx(r)Dhexp|- ^ j d?xd^yj^{i^)-^ 



1 /■ ,S -2 V^^^ 



>S'int.(-^! Jm) 



(125) 



The interaction term of the electric vortex line with the monopole current now takes the form 



e^^x d xd ?/j^(x) 



(y 



'-5 



This interaction term can be again rewritten as ^int. = / rf^xjH™''', where the magnetic induction, 
generated by the electric vortex hne, obeys the equation e^uxd^Hl"'^- = 6^. Note again that if we 
set for the monopole current the classical expression ( |119| ), the interaction term takes the form 
^int. = 'igm.L{L, r) with r) standing for the Gauss linking number of the contours L and F. 

In conclusion of this Section notice that our results demonstrate the usefulness of the Abelian 
projection method and path-integral duality transformation for the solution of the problem of 
string representation of non-Abelian gauge theories and provide us with some new insights con- 
cerning the vacuum structure of these theories. They give a new field-theoretical status to the 
MFC as well as to the 't Hooft-Mandelstam scenario of confinement. 



4 Ensembles of Topological Defects in the Abelian- 
Projected Theories and String Representation of Com- 
pact QED 

In the previous Section, we have demonstrated the relevance of the massive Kalb-Ramond field 
coupled to the string world-sheet to the construction of the string representation of Abelian- 
projected SU (2)-gluodynamics. We have also argued that the large-distance asymptotic behaviour 
of the propagator of this field matches the one of the bilocal field strength correlator in MFC of 
QCD. All this tells us that the Kalb-Ramond field coupled to the string world-sheet is indeed 
quite adequate for modelling the QCD string effective action. 

In the present Section, we shall find the dual formulation of one more model allowing for the 
analytical description of confinement, which is compact QED |p.03| , in three and four space-time 



dimensions. There, confinement is also due to monopoles, which however form not a condensate 
of their Cooper pairs (as it is argued to take place in the Abelian-projected theories), but a dilute 
gas. The resulting dual action also turns out to be some kind of a Kalb-Ramond field action, albeit 
quite nonlinear one. As we shall see, this nonlinearity will eventually realize the independence of 
the Wilson loop, describing an external electrically charged test particle, of the form of a certain 
surface bounded by the contour of this Wilson loop. Namely, it will be demonstrated that this 
independence is achieved by the summation over the branches of a certain (highly nonlinear) 
multivalued effective potential of the monopole densities in 3D or monopole currents in 4D, both 
of which are unambiguously related to the Kalb-Ramond field via the modified Bianchi identities. 
This is the essence of the string representation of compact QED, which is thus alternative to the 
DAHM one. Note that the inapplicability of the method of a derivation of the string representation 
of DAHM to compact QED is clearly not surprising due to the absence of magnetic Higgs field in 
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the latter case (no condensation of monopole Cooper pairs). Correspondingly, there is no way of 
getting the integration over string world-sheets from the integration over the singular part of the 
phase of this field. 

However in the low-energy limit, the real branch of the above mentioned effective monopole 
potential goes over into a simple quadratic functional, and the resulting action of the Kalb-Ramond 
field becomes analogous to the one of Eq. (^^. This will enable us to derive the string tension 
and the bilocal field strength correlator in 3D- and 4D compact QED in this limit. Next, we 
shall apply the partition function of the type ( [5^ ) to the description of the string world-sheet 
excitations in DAHM, 4D compact QED, and QCD and derive the effective action quadratic in 
these excitations. In our interpretation of the above mentioned topics, which will be covered in 
the next two Subsections, we shall mainly follow the papers |]88| , |132|| and [|131| . 



In the third Subsection, we shall demonstrate that the developed techniques of investigation 
of the grand canonical ensembles of monopoles are quite appropriate to the description of similar 
ensembles of topological defects in the Abelian-projected theories. This is a natural extension 
of the results of the previous Section, where these defects have been treated as individual {i.e., 
noninteracting) ones. Here, contrary to that, collective effects in the dilute gases of such defects 
will be studied. In our interpretation of this topic, we shall follow Refs. ||133|, |134 . 



4.1 Vacuum Correlators and String Representation of Compact QED 
4.1.1 3D-case 

In the present Subsection, we shall derive a representation of 3D compact QED in terms of 
the monopole densities and then employ it for the construction of the string representation of 
the Wilson loop in this theory. Besides that, vacuum correlators in the low-energy limit of 3D 
compact QED will be also investigated. 

The most important feature of 3D compact QED, which distinguishes it from the noncompact 
case, is the existence of magnetic monopoles. Their general configuration is the Coulomb gas with 



the action [103 



Smon. = 9'^^qaqb{^ ^) (z^, Zfe) + 5*0 9a) (126) 

a<b a 

where A is the 3D Laplace operator, and 5*0 is the action of a single monopole, 5*0 = Here, 
similarly to Ref. [|103|1 , we have adopted the standard Dirac notations, where eg = 2™, restricting 
ourselves to the monopoles of the minimal charge, i.e., setting n = 1. Then, the partition function 



of the grand canonical ensemble of monopoles corresponding to the action (|126|) reads 



c 



N N 

=1+E E TT^n 

N=lqa=±l • j=l 



d^Zi exp 



71 

2^ 



FPgas(y) 



(127) 



N 



where Pgas(x) = J2 (x — z^) is the monopole density, corresponding to the gas configuration. 

a=l 

Here, a single monopole weight C oc exp (— Sq) has the dimension of (mass)^ and is usually referred 
to as fugacity. Notice also that we have restricted ourselves to the values = ±1, since at large 
values of the magnetic coupling constant g, monopoles with Ig] > 1 turn out to be unstable and 
tend to dissociate into those with |g| = 1. 
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Next, Coulomb interaction can be made local, albeit nonlinear one, by introduction an auxiliary 
scalar field ||103 | 



Zxx 



Dxexp 



2 (^mX)^ - 2Ccos(^x) 



The magnetic mass m = g\/2C, of the field x, following from the quadratic term in the expansion of 
the cosine on the R.H.S. of Eq. 



[28|) , is due to the Debye screening in the monopole gas. The next, 
quartic, term of the expansion determines the coupling constant of the diagrammatic expansion 
for this gas, which is therefore exponentially small and proportional to (^''exp (—const. (7^). 

Before proceeding with the investigation of 3D compact QED, let us present the derivation of 
its partition function ( |127| ) (as well as the related Eq. ( |128| )) from the partition function of the 
usual 3D QED extended by the external monopoles. The statistical weight of such a model {i.e., 
its partition function for the given monopole charges ga's and the number of monopoles A^, where 
the summation over the grand canonical ensemble is not yet performed) is the same as the one 
of the Ginzburg-Landau theory extended by external monopoles, but without the Higgs field (c/. 
the previous Section) and reads 



DA. 



(129) 



Here, F^^ = d^A^ — d^A^, and F^^ is the monopole field strength tensor, responsible for the 
violation of the Bianchi identities. 



gas- 



(130) 



Next, the term ^F^^ 



in the Lagrangian on the R.H.S. of Eq. (|129|) can be linearized by introducing 
the integration over an auxiliary vector field as follows 



exp 



f d^xF^ 



DB^ exp 



(fx [ + Y^s^,xB^F,x 



After that, the integration over the A^-field leads to a constraint {iee^y\B^ + F^'^) = 0, whose 
resolution yields B^ = j^e^y\Ff;^ + d^x- Substituting this representation for the field B^ into the 

remained action / d^x ^B'l + [P^tf] and making use of Eq. ( |130| ), we finally arrive at the 



following expression for the partition function ( |129|) 



Z [pgas] = j Dxexpl- J d^x 



lid.xr 



WXPg^s 



(131) 



Clearly, the summation over the grand canonical ensemble of monopoles leads now to Eq. 
On the other hand, integration over the field x Eq. ( |131|) yields the argument of the exponent 
standing on the R.H.S. of Eq. ( [127] ). This conclusion completes the derivation of Eqs. (|127|) 
and ( |128| ) from the partition function 



Let us now cast the partition function (|127|) into the form of an integral over the monopole 
densities. This can be done by introducing into Eq. (|127|) a unity of the form 



J DpS (p(x) - Pgas(x)) = J DpDXexp 



N 



La=l 



J d^xXp 



(132) 
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Then, performing the summation over the monopole ensemble in the same way as it has been 
done in a derivation of the representation (|128|) , we get 



DpD\ exp < / (TxiTypiidi 

I 2e^ J |x ^ y| 



p(y) + / d^x (2C cos A - i\p) \ . (133) 



Finally, integrating over the Lagrange multiplier A by resolving the corresponding saddle-point 
equation. 



sm A = -, 

we arrive at the following expression for the partition function 



Dp exp 



Here, 



V^real[p] = J d^X lp\n 



2C 



2C 



1 + 



(134) 



(135) 



(136) 



is the parabolic-type effective monopole potential, whose asymptotic behaviour at p ^ C reads 



Koal[p] 



(137) 



Owing to Eq. ( |132| ), the obtained representation ( p.35| ) is natural to be called as a representation 
for the partition function in terms of the monopole densities. 

Note that the reason for the notation "Keai[p]" is that during the integration over the field A in 
Eq. ( |133| ) we have chosen only the real branch of the solution to the saddle-point equation ( [134| ) and 



disregarded the complex ones. Owing to this, Eq. (|135|) describes actually not the full expression 
for the partition function, but only its real part. The corresponding full expression will be discussed 
below. 

The obtained representation for the partition function can be straightforwardly applied to 
the calculation of the coefficient function 1)™°°- (x^), related to the bilocal correlator of the field 
strength tensors as follows 



+ 



x,,6x 



V'f (x^ 



(138) 



where the average over the monopole densities is defined by the partition function ( |135| ), whereas 
the y4^-average is defined as 



JDAp{...)exp [-^Jd'xF^ 



JDApexp [-^Jd^xF^ 
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In Eq. ( |138| ), J-'fj,^ = F^^ + stands for the full electromagnetic field strength tensor, which 
includes also the monopole part 



This monopole part yields the R.H.S. of the Bianchi identities modified by the monopoles, 

df,n^ = 2np, (139) 

where TC^ = \e^^x^ux stands for the full magnetic induction. Eqs. (|138|) and (|139|) then lead to 
the following equation for the function p^o"^- 

AI?---(x2) =-4vr2(p(x)p(0))^, (140) 

which is a 3D analogue of the 4D equation ( |102|) . The correlator standing on the R.H.S. of 
Eq. ( [L40| ) can be found in the limit of small monopole densities (low-energy limit), p <^ (. By 
making use of Eqs. ( |135| ) and ( [L37| ), we obtain 



/■ p-m|x| 

(p(x)p(O))^ = -f 

^ ZTT |x| 

Then, demanding that X>'^°'^' (x^ oo) — * 0, we get by the maximum principle for the harmonic 
functions the desired expression for the function X)"^"^' in the low-energy limit 

, . p-m|x| 
pmon. ^^2^j ^ 27rC^^. (141) 

We see that in the model under study, the correlation length of the vacuum Tg, i.e., the distance 
at which the function decreases, corresponds to the inverse Debye mass of the field "^"^ 

(c/. the case of Abelian-projected theories, where the role of Tg was played by the inverse mass of 
the dual gauge boson(s), generated by the Higgs mechanism). The coefficient function "D^^" (x^) 
will be derived later on. 

Let us now proceed to the problem of string representation of 3D compact QED. To this end, 
we shall consider an expression for the Wilson loop in this theory and try to find a mechanism 
realizing the independence of this object of the form of a certain string world-sheet S, bounded by 
the contour C. By virtue of the Stokes theorem, the Wilson loop can be rewritten in the following 
form 

{W{C)) = (exp r-J da^^J^^A \ = /exp li J dcx^nA\ 

= {W{C))^^ (exp J d'xpi^U^)^ )^ , (142) 
where the free photon contribution reads 



{WiC))^^ = Lxpiif A^dxA) =exp . (143) 
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In Eq. (|ig), da, 



\e^,u\dayx, and r7(x) 



stands for the solid angle under 



which the surface S shows up to an observer at the point x. Notice that due to the Gauss law, 
in the case when S is a closed surface surrounding the point x, '/^(x) = 47r, which is the standard 
result for the total solid angle in 3D. 



Equation ( 142 ) seems to contain some discrepancy, since its L.H.S. depends only on the con- 
tour C, whereas the R.H.S. depends on an arbitrary surface S, bounded by C. However, this 
actually occurs to be not a discrepancy, but a key point in the construction of the desired string 
representation. The resolution of the apparent paradox lies in the observation that during the 
derivation of the effective monopole potential ( |136| ), we have accounted only for the one, namely 
real, branch of the solution to the saddle-point equation ( p.34| ). Actually, however, one should sum 
up over all the complex-valued branches of the integrand of the effective potential (|136D at every 
space point x. This requires to replace Keai[p] by the total expression 



+00 . 
V[p]= jd 



X < p 



In 




+ 27rm 



-2C 




(144) 



which, in particular, being substituted for l^eai[p] into Eq. ( |135| ) yields the full partition function, 
mentioned above. As far as the Wilson loop, given by the expression 



{W{C)) = {W{C)). X 



Dpexp 



TC 

2? 



FP(y) + v[p] 



d^xp{'x.)ri(x.) 



(145) 



is concerned, such a summation over the branches of V[p] thus restores the independence of 
this object of the choice of the world-sheet. (Notice that from now on we omit an inessential 
normalization factor, implying everywhere the normalization {W{0)) = 1.) 

It is worth noting that the obtained string representation (|145| ) has been for the first time 
derived in another, more indirect, way in Ref. [p8|. It is therefore instructive to establish a 



correspondence between our approach and the one of that paper. 

The main idea of Ref. |^ was to calculate the Wilson loop starting with the direct definition 
of this average in the sense of the partition function ( |127| ) of the monopole gas. The corresponding 
expression has the form 



1+ 



-N N 



+ Y1 JJ^Ii d^z^exp -— / d^xd^ypg^i^) _ Pgas(y) + - / d^xpg^{x)r]{x) 

N=lqa=±l ■ i=l'' IX VI /J 



J Dxexpi^- J d 
Dip exp i — I d^x 



eW 1 



(146) 
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where (p = gx + 2- 

Next, one can prove the following equality 



exp 



in 



dx^ f dy^ 



c 



c 



d X d^if - -8^7] 



j Dh^^ exp - J d^x {i(pef,^xd^,Kx + g^h^^^ - 2mh^^'i:^^ 



(147) 



which makes it possible to represent the contribution of the kinetic term on the R.H.S. of Eq. (|146|) 
and the free photon contribution ( |143| ) to the Wilson loop as an integral over the Kalb-Ramond 
field. The only nontrivial point necessary to prove the equality (|147|) is an expression for the 
derivative of the solid angle. One has 



dxvi^) 



_d_ 

dyx 



da xiy) 



d \ d 



dyu J dyr, |x - y 



1 r 1 

r + / daxiy)A- -. 

-y J x-y 



(148) 



Applying to the first integral on the R.H.S. of Eq. 

. d 



da,,^ daxTT- 



Stokes theorem in the operator form. 



£^J.Xudyu, one finally obtains 
d r , 1 



dxvi^) 



ex 



dx. 



dyv 



An J daxiy)6{^-y). 



c 



exp^-- 



Making use of this result and carrying out the Gaussian integral over the field hfj_^, one can 
demonstrate that both sides of Eq. ( |147|) are equal to 

^ j d^x {d^Lp f + ^J da^d^(p + j rfo-^(x) j rfo-^(y)5(x - y) 

S EE 

thus proving the validity of this equation. 

Substituting now Eq. ( |147| ) into Eq. ( |146| ), it is straightforward to carry out the integral over 
the field if. Since this field has no more kinetic term, such an integration can be performed in 
the saddle-point approximation. The corresponding saddle-point equation has the same form as 
Eq. ( p.34| ) with the replacement p £fj.u\d^hux- As a result, we obtain the following expression 
for the full Wilson loop 

{W{C)) = {W{C)). {W{C))^^^^ = 



J Dh^^exp <- J d^x {g^h?^^ + V [e^,yxd^,Kxi) + 27ri j da^^h^^ 



(149) 



where the world-sheet independence of the R.H.S. is again provided by the summation over the 
branches of the multivalued action, which is now the action of the Kalb-Ramond field. 

Comparing now Eqs. ( |145|) and ( 149 ), we see that the Kalb-Ramond field is indeed related 
to the monopole density via the equation e^yxd^h^x = P- Thus, a conclusion following from the 
representation of the full Wilson loop in terms of the integral over the Kalb-Ramond field is 
that this field is simply related to the sum of the photon and monopole field strength tensors as 
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h^u = '^^i^u- In the formal language, such a decomposition of the Kalb-Ramond field is just the 
essence of the Hodge decomposition theorem. 

Let us now consider the low-energy limit of Eq. (|149|) and again restrict ourselves to the real 
branch of the effective potential, i.e., replace V [e^^xdfih^x] by Vj-ca.i[e^^xd^h^x]. This yields the 
following expression for the Wilson loop 



^ H%x + a'^hl^^ - 27ri/i^^S^^ 



(150) 



.6C 

Notice, that the mass of the Kalb-Ramond field resulting from this equation is equal to the Debye 
mass m of the field x from Eq. ( |128| ). 

One can now see that Eq. ( [L50| ) is quite similar to the 3D version of Eq. ( |114| ) (with the A^- 
field gauged away) we had in the DAHM case. However, the important difference from DAHM 
is that restricting ourselves to the real branch of the potential, we have violated the surface 
independence of the R.H.S. of Eq. ( |150| ). This problem is similar to the one we met within the 
MFC in Subsection 1.3, where in the expression for the Wilson loop, written via the non-Abelian 
Stokes theorem and cumulant expansion, all the cumulants higher than the bilocal one were 
disregarded. There, the surface independence was restored by replacing E by the surface of the 
minimal area, Smin. = Smin. [C], bounded by the contour C. Let us follow this recipe, after which 
the quantity 

S^,,. = - In (W^(C)),__^^|^^^_ (151) 

can be considered as a low-energy string effective action of 3D compact QED. 

The integration over the Kalb-Ramond field in Eq. ( |150[ ) is now almost the same as the one 
of Appendix C and yields 



(W^(C))low-energyL^^_ = ^Xp -\ J daxA^) J da,,{y) {J^xA^)J^Uy))A,„ 

\ ^min. ^min. 

which is consistent with the result following directly from the cumulant expansion of Eq. (|142|) . 
Here, the bilocal correlator is defined by Eq. ( |138| ) with the function 25™°"^- given by Eq. ( |141| ) and 
pfuii = j)P^°^- _|_ j)mon.^ where the photon and monopole contributions read 



pphot. |^^2^ 

and 



e2 



27r|xP 



(a;^) = -^L+^)e-™N, (152) 



respectively. Since the approximation p <^ C, in which Eq. ( |141| ) has been derived, is just the 
low-energy limit, in which Eq. ( |150| ) follows from Eq. (|149|) , coincidence of the function T>'^°'^', 
following from the propagator of the Kalb-Ramond field, with the one of Eq. ( 141 ) confirms the 
consistency of our calculations. 
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Notice that by performing an expansion of the nonlocal string effective action ( |151D in powers 
of the derivatives w.r.t. the world-sheet coordinates ^°'s, one gets the following string tension of 
the Nambu-Goto term 



a = 7r'^. (153) 
9 

Similarly to the expression ( ]86| ) for the string tension in the Abelian-projected SU{2)- 
gluodynamics, which is nonanalytic in the magnetic coupling constant, Eq. (|153|) is nonanalytic 
in the electric coupling constant. This observation manifests the nonperturbative nature of string 
representation of these theories. 

Note that the inverse bare coupling constant of the rigidity term corresponding to the string 
effective action (|151| ) can be also calculated by virtue of the results of Section 2, and the result 



reads = — ^ ^ 3 • This expression is also nonanalytic in e, and its negative sign supports 



the stability of strings. However, as it has been demonstrated in Section 2, at the surface of the 
minimal area, Smin. [C] , at which string effective action (|151|) is defined, the rigidity term vanishes. 
We see that the long- and short distance asymptotic behaviours of the functions ( |141| ) and ( |152| ) 



display the properties analogous to those of the corresponding functions in QCD within MFC pjl] , 
|62| . Namely, at large distances both of the functions ( |141| ) and ( p.52| ) decrease exponentially 



with the correlation length m~ , and at such distances X>^°'^' -C T'™°'^- due to the preexponential 
factor. In the same time, in the opposite case |x| <^ m~^, the function p™""^- is much larger 
than the function T>"^°'^\ which also parallels the MFC results. Notice, however, that the short- 
distance similarity takes place only to the lowest order of perturbation theory in QCD, where its 
specific non-Abelian properties are not important. Note also that the above described asymptotic 
behaviours of the functions T>^°^' and p^o^^- clearly match those of the corresponding functions, 
which parametrize the bilocal correlator of the dual field strength tensors in the Abelian-projected 
SU{2)- and S'[/(3)-gluodynamics (c/. the previous Section). 

4.1.2 4D-case 

In the present Subsection, we shall generalize the above considerations to the case of 4D compact 
QED. Note that in this case, contrary to what happens in 3D, monopoles are no more point- 
like particles at rest, but are rather represented by the closed loops. The disorder effects in the 
ensemble of such loops, leading to the confinement of an external electrically charged test particle, 
become strong enough only in the so-called strong coupling limit, i.e., at sufficiently large values 
of the electric coupling constant. This situation differs from the case of 3D compact QED, where 
confinement takes place at arbitrary values of the coupling constant (see Ref. [0 for an extended 
discussion). 

The action we start with is again of the form of AHM without the Higgs field, extended by 
external monopoles (c/. Eq. 



'S'gas 



\jd'x{F,^ + F^:;)\ (154) 



Here, the monopole field strength tensor describes the violation of the Bianchi identities by the 
collective current of monopoles, which form a dilute gas 
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N 



a=l 



(155) 



In Eq. (|155|) , we have parametrized the trajectory of the a-th monopole by the vector a;^(r) = 

1 

y'fi + where = J drx'^^ij) denotes the position of the trajectory, whereas the vector z^j^ij) 

corresponds to its shape, both of which should be eventually averaged over. The procedure of a 
derivation of the statistical weight entering the partition function of the grand canonical ensemble 
of monopoles is analogous to the corresponding 3D one, described in the previous Subsection. 
Firstly, one can linearize the term ^ / d^xF"^^ in the action ( |154|) by introducing an integration 
over an auxiliary antisymmetric tensor field B^^ as follows 



DA^DB^y exp 



d X 



B,„. + -Bfj^yF^j, 



.2 



1 1 / n2- 

^ 2 t^" fJ.i' ^ ^ \ tJ-i^ J 



Next, the constraint 



dfj. {b^iv 



resolved by setting B 



+ F^u, where F^, 



emerging after the integration over the field can be 
dyB^. This finally yields 



d^B, 



Z 



,-gas 
J fj, 



DBn exp 



d X 



_ P2 _ • n -gas 



This expression should now be averaged over the monopole gas in the following sense 



O 



gas 



N 

n 

i=l 



dYF>z'fx 



ga=±l 



(156) 



(157) 



Here, /i [z"^] is a certain rotation- and translation invariant measure of integration over the shapes 
of monopole loops, whose concrete form will not be specified here (For example, one can choose 
it in the form of the properly normalized measure of an ensemble of random loops, representing 
trajectories of scalar particles, i.e., 



Dz'fi 



O 



dsj 

Si 



Du{s[)exp -- / ii^is'^ds'i 0[u{s[)] 



u{0)=u{si) 

where the vector u^{s'j) parametrizes the same contour as the vector z^^{t).). 

One can now write down the contribution of monopoles to the partition function of their 
grand canonical ensemble. Owing to Eq. (|157|) it reads 



1+E 

N=l 



(20 



N 



J d'^y j Dzulz] cos (^g dz^B^{x] 



N 



(158) 



Here, C, oc e~^° again stands for the fugacity (Boltzmann factor of a single monopole) of dimension 
(mass)^ with the action of a single monopole given by 5*0 = ^^if^. Note that the constant in this 
formula is proportional to the length of the monopole loop. 
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In order to evaluate the path-integral over z^'s in Eq. ( |158D , let us employ the dilute gas ap- 
proximation, which requires that typical distances between the monopole loops are much larger 
than their sizes. This means that generally \y"'\ ^ \z°'\^ where \y\ = \fl/^i- Let us denote charac- 
teristic distances \y\ by L, characteristic sizes of monopole loops = /c/r-\/i^) by a, and perform 

V / 

the Taylor expansion of B^{x) up to the first order m. a/L (which is the first one yielding a 
nonvanishing contribution to the integral § dz^B^{x) on the R.H.S. of Eq. ( |158|) ), 



B^{x) = B,{y) + z,d,B^{y) + O (^(^-j j • (159) 
Then, the substitution of this expansion into Eq. ( |158| ) yields 

J Dz^[z]cos {g j dz^B^{x)^ ~ j Dz^[z] cos {gP^^[z]d„Bf,{y)) = 

oo (_l)n I- 

= ll-7A^9'''d,.B,M---du,^B,,M Dzf,[z]P,,,,[z]---P,,^,Jz], (160) 

n=0 l^^-J- 

where P^u[z] = § dz^Zy stands for the tensor area associated with the contour parametrized by 
Due to the rotation- and translation invariance of the measure /^[z], the average of the 
product of the tensor areas can be written in the form 



/[a ) r- 
Dz^{z\P^^y^[z\ ■ ■ ■ PfMinUinM = (^2n — 1)!! ■ ■ ■ lA'2n-ii'2n-i,At2n!^2n + pcrmutatious 

(161) 

Here, l^u,Xp = | {S^ixSup — ^pp^ux) {cf. Appendix C), and the normalization factor (2n — 1)!! is 
explicitly extracted out since the sum in square brackets on the R.H.S. of Eq. ( |161| ) contains 
(2n — 1)!! terms. Substituting now Eq. ( |161| ) into Eq. ( |160D and contracting the indices, it is 
worth noting the following. Due to the Hodge decomposition theorem, the field -B^ can always be 
represented as d^ip + efj,iy\pduhxp, where (p and hxp stand for a scalar and an antisymmetric tensors, 
respectively. Owing to the conservation of the current j^^^, the contribution of the first of them 
to the partition function ( |156| ) vanishes. The remaining part of the 5^-field automatically obeys 



the condition = 0. Taking this into account, estimating [d^Byf' as and denoting 

(<^ a) by A~^, where A acts as a natural UV momentum cutoff, we finally obtain 

\Bp{y)\\ 



Dzfi[z] cos (^g (j) dz^B^{x)^ 



cos 



A 



This yields the desired expression for the partition function of the grand canonical ensemble of 
monopoles 



Zmon. = / DB^ exp <^ - / d^x 



^ , (162) 

The Debye mass of the -B^-field following from the expansion of the cosine in the action on the 




R.H.S. of Eq. (|T62]) reads m = 

■^■^One can check that for the plane contour, = —Puf_i = —S, fi < v, where 5* is the area inside the contour. 
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Let us now turn ourselves to a derivation of the representation for the partition function ( |162D 
in terms of the monopole currents similar to the representation of its 3D analogue (|128|) in terms 
of the monopole densities, given by Eqs. (|135|) and ( p.44|) . Such a representation will enable us to 



obtain the bilocal correlator of the field strength tensors and to construct the string representation 
for the Wilson loop of an external electrically charged test particle in the 4D-case under study. To 
proceed with note that, analogously to the 3D-case, the integration over the field in Eq. ( |1 561 ) 
leads to the Coulomb interaction between the monopole currents. 









J fj. 


= exp 





1 



{x — x') 



Averaging this expression over the grand canonical ensemble of monopoles we obtain 

oo 



J Dj^DX^,exp 



X — X' 



/\2 



'. j d'^xXf.j^ + 2C J d' 



X cos 



A 



(163) 

where the term fixing the Fock-Schwinger gauge for the Lagrange multiplier is assumed to be 
included into the integration measure. Notice that Dj^ here is the standard integration measure 
over the vector field, which is of the same form as -DA^. 

Clearly in order to integrate now over A^, one should solve the saddle-point equation 



lA, 



sm 




iA 
'2C 



(164) 



This can be done by noting that its L.H.S. is a vector in the direction A^, which means that it 
can be equal to the R.H.S. only provided that the direction of the vector A^ coincides with the 
direction of the vector j^. Therefore, it is reasonable to seek for a solution to Eq. ( |164| ) in the form 



A, 



IA 



^1 li^l 



riiy . Then, Eq. (|164|) reduces to the scalar equation sin 



solution of the latter one yields the desired current representation 



2C 



1 



d'^xd'^x' 



1 



Straightforward 



(165) 



where the effective complex-valued potential of monopole currents reads 



+ 00 

E ld'x{A\j, 



In 



A 



1 + 



'A. 
2C 









+ 27rm 


-2C^ 





(166) 



Analogously to the 3D-case, by virtue of the representation ( [165| ), ( |166| ) it is now possible to 
construct the string representation for the Wilson loop of an external electrically charged test 
particle. Once being written via the Stokes theorem, the Wilson loop reads 
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{WiC)) = /expr-lda,,J^,A\ ={WiC))^/expr-lda,,h,A\ , (167) 




where JF^,^ = eF^j, + h^^ is the full field strength tensor, extended by the Kalb-Ramond field h^^, 
which obeys the equation d^h^y = j^. Next, on the R.H.S. of Eq. ( |167| ), the A^-average is defined 
as 



(...) 



JDA,{...)exp{-yd'xF^ 



whereas the j^-average is defined by means of the partition function ( |165| )-( [T66| ). Among the 
two multipliers on the R.H.S. of Eq. ( p.67| ), the first one is the standard "perimeter" (Gaussian) 
contribution to the Wilson loop, brought about by the free photons, i.e., 



{WiC)). = exp 



c c / 



87r2 

As far as the second one is concerned, expressing h^^ via j^, it is possible to rewrite it directly as 

expi^-i J d'^xj^rii^^ , (168) 



where 



l_ _d_ 



^^(^) = ^^M^VT^ / daxpixi^))— — ^^7^ (169) 



s 



stands for the 4D solid angle, under which the surface S is seen by an observer at the point x (If S is 
a closed surface surrounding the point x than by virtue of the Gauss law, da^u dS^d^ — dS^dfj,, 
one can check that for the conserved current j^, / d'^xjfj,!]^ = J dS^j^, as it should be. Here, 
dS^ stands for the oriented element of the hypersurface bounded by S.). Again, an apparent S- 
dependence of Eq. ( |168| ) actually drops out due to the summation over branches of the multivalued 
potential (|166|) . This is the essence of the string representation of the Wilson loop in the 4D 
monopole gas. 

Let us now consider the low-energy limit of the Wilson loop, i.e., the limit A ^ (, and 
investigate the (stable) minimum of the real branch of the potential. This corresponds to extracting 
the term with n = from the whole sum in Eq. ( |166|) . Then, since we have restricted ourselves to 
the only one branch of the potential, the S-independence of the Wilson loop is spoiled. In order 
to restore it, let us choose E to be the surface of the minimal area S = Smin. [C]. Then the Wilson 
loop takes the form 

(W^(^))low-ener«v = (C)) . X 



X 



(170) 
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where now rj^ is defined by Eq. ( |169D with the replacement S Smin.- Recalhng the expression 



for jp via h^u, Eq. ( [17CI| ) can be written as follows 



(W^(^))low-energy = {C)) J Dh,^ CXp - / d'x (^^^^A + ^/^J.) + ^ / ^^M^ V 



(171) 

Note that the mass of the Kalb-Ramond field following from the quadratic part of the action 
standing in the exponent on the R.H.S. of Eq. (|171|) is equal to the Debye mass m of the field 
following from Eq. (|162|) . 

Integrating over the Kalb-Ramond field we obtain 

(W^(C))iow-energy =exp| J d(rxu{x) J da pp{y) {Txu{x)T pp{y)) ^^ -^ 

\ ^min. ^min. 

Here, owing to the Appendix C, the bilocal correlator, once being parametrized as 

{^xu{x):Fppm^^^^^ = (6xp6,p - 6xp6,p)v^"-- {x') + 



1 

+ 2 



XpSxp Xp5xp 



with P^^" = x>f^"^' + T'™°'^-, yields for the coefficient functions T>™°'^' and the expressions 

given by Eqs. ( [^ ) and (0), respectively, with the Higgs mass of the dual gauge boson replaced 
by the Debye one. Consequently, similarly to the case of the Abelian-projected theories, the large- 
(and short) distance asymptotic behaviours of the functions D'^°'^' and J)™^'^' match those of the 
corresponding functions in QCD within the MFC (see discussion after Eq. (|92|)). We also see that 
in the case of 4D compact QED under study, the role of the correlation length of the vacuum Tg 
is played by the inverse Debye mass of the -B^-field, . As far as the function pP^°*- (^x"^^ is 

2 ____ 

concerned, it turns out to be equal to g^||^|4 (c/. Eq. (|37D). All this demonstrates the relevance of 
concepts of the MFC to the description of confinement in compact QED. 

4.2 A Method of Description of the String World-Sheet Excitations 

The results of the previous Sections tell us that strings in QCD within the MFC and the approach 
based on the method of Abelian projections, as well as in compact QED, can be with a good 
accuracy described by the same action of the massive Kalb-Ramond field interacting with the 
string world-sheet. This makes it reasonable to develop a unified mechanism of description of the 
world-sheet excitations in all these theories. This can be done on the basis of the background-field 



method proposed in Ref. |135| for the nonlinear sigma models. 



Let us for concreteness work with the action standing in the exponent on the R.H.S. of Eq. (|84D 
for the case when there are no external quarks, i.e., the string world-sheets are closed. This action 
has the form 

S = U'x {-^^H%x + glK, - ^^/^M^S^,) , (172) 
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(c/. Eq. ( |114D where the A^-field has been absorbed by fixing the gauge of the field hf^^). In 
order to develop the background-field method, one should define a geodesies passing through the 
background world-sheet and the excited one x^{^) = y^{C,) + where Zfj,{^) stands for 

the world-sheet fiuctuation. Such a geodesies has the form p^{^,s) = + sZfj_{^), where s 

denotes the arc-length parameter, < s < 1. The expansion of the string effective action (|172|) 
in powers of quantum fluctuations z^{^)^s can be performed by virtue of the arc-dependent term 
describing the interaction of the Kalb-Ramond fleld with the string, which reads 

S[pit s)] = -m I rf'e V[p(e, s)y' {dap.it s)) (9,p.(e, s)) . 
Then the term containing n quantum fluctuations has the form 

1 d"" 



n\ ds'' 



s=0 



In particular, we obtain 



SW = / daMm^ulvm, 



(173) 



(174) 



and 



-zvr / d'^z,{Oe'''{day,iO)[idbZxmHu,x[y{0] + ^^a{OidM^^ (175) 



where, owing to the closeness of S, during the derivation of Eqs. ( |174|) and ( |175|) several full 
derivative terms were omitted. 

Notice that as it was discussed in Ref. ||135|| , the terms ( |173| )-( pT5| ) are necessary and sufficient 
to determine all one- loop UV divergencies in the theory ( |172| ). That is why in what follows we 
shall restrict ourselves to the derivation of the effective action, quadratic in quantum fluctuations 

In order to get such an action, we shall flrst carry out the integral 



J Dh^y exp 
It turns out to be equal to 



— d X 



exp 



dyaiO' 



dyaiOdyfsi^'Y 



;i76) 



Here, we have again taken into account that the world-sheets under study are closed, so that the 
boundary terms vanish. 
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Secondly, one should substitute the saddle-point value of the integral 



J Dhf^i, exp 



- / d^x 



HfluX + dm^fiu ] ^^^^ 



into Eq. ( |175| ). This saddle-point value reads 



(177) 



where m = 2(7^^ is again the mass of the Kalb-Ramond field. Upon the substitution of Eq. (|177|) 
into Eq. (|175|) and accounting for Eq. (|176|) , we finally get the following action quadratic in 
quantum fluctuations 



S, 



quadr. 



{nr]f 



{y - y')a{y - y')p 
{y - y'f 



3V, {{y - y'y) + 



m 



57r^ \y — y 



- {3Ki (m \y - y'\) + K3 (m \y - y'\)) 



d_ 

'dya 




{y'){y - y')v + ^ ^i, \, \ v) + {X ^i, u \, ^ u)\Vi [{y - y'f^ 



where y = y{^), y' = y{^'), and the function Vi is defined by Eq. (|92D . It is remarkable that though 
the interactions between the points lying on the background world-sheet are completely described 
via the function T>, the dynamics of the world-sheet fluctuations is governed by the function 
Vi, which in the case of open world-sheets is responsible for the perimeter-type interactions. 
This phenomenon can be interpreted as an interpolation between the world-sheet and world-line 
dynamics, which is absent on the background level. 



4.3 Abelian-Projected Theories as Ensembles of Vortex Loops and 
Modified Expressions for the Field Strength Correlators 

In our investigations of Abelian-projected theories performed in the previous Section, the interac- 
tion of topological defects (Abrikosov type electric vortices in 3D and Nielsen-Olesen type strings 
in 4D) has not been taken into account, i.e., these objects have been treated as individual ones. 
Owing to that, string representations of the partition functions of Abelian-projected theories as 
well as field strength correlators in these theories were insensitive to the properties of the ensemble 
of strings as a whole and depended actually on a certain fixed string configuration only {e.g. a 
single string). On the other hand, it is known [|119|| that in the case of zero temperature under 
study Abrikosov vortices in the Ginzburg-Landau theory form bound states, consisting of a vortex 
and an antivortex, which are usually referred to as vortex dipoles. Such vortex dipoles are short 
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living (virtual) objects, whose typical sizes are much smaller than the typical distances between 
them. This means that similarly to monopoles in compact QED, vortex dipoles form a dilute 
gas. In the 2D-case, the summation over the grand canonical ensemble of such dipoles in the 
dilute gas approximation has been performed in Ref. [|136|| , and the result has the form of the 2D 



sine-Gordon theory of the scalar field similar to Eq. (|128|) , albeit with an additional mass term of 
the field x- The aim of the present Subsection is to perform analogous summations in the 3D- and 
4D-cases of the Abelian-projected SU{2)- and SU (3)-gluodynamics. Besides the representations 
in terms of the effective sine-Gordon theories, the related direct representations in terms of the 
integrals over the vortex dipoles (small vortex loops built out of two Nielsen-Olesen type electric 
strings in the 4D-case) similar to Eqs. (|135|) , ( |144| ) (( p.65|) , (|166|) in 4D) will be derived, which will 
enable us to evaluate correlation functions of vortex dipoles (loops). Such correlation functions 
will further be applied to the more precise determination of the correlators of field strength tensors 
in Abelian-projected theories. In particular, it will be shown that the effect of Debye screening in 
the gas of vortex dipoles (loops) leads to the modification of the correlation length of the vacuum, 
at which these correlators decrease. Namely, it changes from the inverse mass of the dual gauge 
boson, which it acquires due to the Higgs mechanism, to its inverse full mass, which accounts also 
for the above mentioned Debye screening effect. 

Let us start our analysis with the dual Ginzburg-Landau theory (3D Abelian-projected SU{2)- 
gluodynamics) in the London limit, whose partition function is given by Eq. ( |123| ). As it has 
been demonstrated in the previous Section, the path-integral duality transformation casts this 
partition function into the form (|124|) . After gauging away the field (p by the gauge transformation 



h ^ h — ^"^^Vv? the partition function takes the form 



Z = J Dx(r)Dh exp I - y d 



where we have denoted for brevity Hf^,^ = d^hy — duh^j,. Note that according to the Hodge 
decomposition theorem, the field /i^ can always be represented as d^cp + e^u\duip\- Due to the 
closeness of vortices, the field decouples from 5^. Moreover, the 0- field decouples from the 
field ip^ as well and yields only an inessential determinant factor, which is not of our interest. 
Therefore, this field can be disregarded from the very beginning, which means that the h-field 
obeys the equation Vh = 0. 

To proceed from individual vortices to the grand canonical ensemble of vortex dipoles, one 
should replace 5^ in Eq. (|178| ) by the following expression describing the gas of N vortex dipoles 



N 

W = E / dzl{T)5 (x - x'^(r)) . (179) 

a=l 

Here, n^'s stand for winding numbers, and we have decomposed the vector x°(r) as x"(r) = 

1 

y° -|- z'^(r), where y*^ = / (irx°(r) denotes the position of the a-th vortex dipole. In what follows, 



we shall restrict ourselves to the vortices possessing the minimal winding numbers, = ±1. 



That is because the energy of a single vortex is known to be a quadratic function of the flux ||119|| , 

and therefore the existence of two vortices of a unit flux is more energetically favorable than the 

existence of one vortex of the double flux. Besides that, as it has been mentioned above, we shall 

work in the dilute gas approximation, according to which characteristic distances |y|, which we 

1 

shall denote by L, are much larger than the characteristic sizes of vortex dipoles, / rfrv^, which 
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we shall denote by a. 

Within these two approximations, by substituting Eq. (|179|) into Eq. ( |178| ) one can proceed 
with the summation over the grand canonical ensemble of vortex dipoles. This procedure can be 
performed with quite general form of the measure of integration over the shapes of vortex dipoles, 
which should obey only the requirements of rotation- and translation invariance. Clearly, such a 
summation essentially parallels the summation over the grand canonical ensemble of monopoles in 
4D compact QED, described in Subsection 4.1. As a result, we arrive at the following expression 
for the grand canonical partition function 



Z = I Dhexp 



_}_ZJ2 ,n 2y^2 




;i80) 



Here, ( o: e ° again denotes the fugacity of a vortex dipole, which has the dimension (mass) 



with So standing for the dipole's action Q, and A 



being the UV momentum cutoff. Thus, 



the summation over the grand canonical ensemble of vortex dipoles, built out of electric Abrikosov 
vortices, with the most general form of the measure of integration over their shapes yields in the 
dilute gas approximation the effective sine-Gordon theory ( |180| ). In particular, this way of treating 
the gas of vortex dipoles leads to increasing of the mass of the field h. Namely, expanding the 
cosine in Eq. ( [18CI| ) we get the square of the full mass, = -|- mj^ = Q'^f]'^, where m = 2gmf] 
is the usual mass of this field (equal to the mass of the dual gauge boson), and mo = is 
the additional contribution coming from the Debye screening. We have also introduced the full 
magnetic charge Q - 



I 2C 



Our next aim is to derive the representation for the partition function ( |180D directly in the 
form of an integral over the vortex dipoles. This can be done by making use of the following 
equality 



exp 



(Tx 



Dj exp <^ - 



d'^xd^yjf.ix)- rj^(y) + 27ri / df'xh^j^ 



in whose derivation it has been used that Vh = 0. Substituting it into Eq. 
straightforwardly resolve the resulting saddle-point equation for the field h, |^ sin [j^ ^ 
which yields the desired vortex representation 



one can 



tA 



Z = Dj exp 



TTT] 



-m|x-y| 

|x-y 



-j,{y) + V[2nj^ 



1811 



Here, the complex- valued potential of the vortex dipoles is given by Eq. ( |166| ) with the replacement 
d^x — i> d^x. 



The obtained representation (|181|) can be now applied to the calculation of correlators of vortex 
dipoles. Indeed, it is possible to demonstrate that if we introduce into Eq. ( |178| ) (with 6^ replaced 
by S^^^) a unity of the form 

^^This action can be estimated as era, where a is the energy of the Abrikosov vortex per its unit length, i.e., the 
3D analogue of the string tension. With the logarithmic accuracy one has |119| a o: if In^ (c/. Eq. (H)), where 
A is the magnetic Higgs field coupling constant. 
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:i82) 



and integrate out all the fields except j, the result will coincide with Eq. ( |181| ). This is the 
reason why the correlators of j's are nothing else, but the correlators of the vortex dipoles. Such 
correlators can be calculated in the low-energy limit, i.e., when A|j| -C C- Moreover, we shall 
perform an additional approximation by restricting ourselves to the real branch of the effective 
potential of vortex dipoles, which corresponds to extracting from the whole sum standing on the 
R.H.S. of Eq. ( |166| ) the term with n = 0. Within these approximations, we arrive at the following 
expression for the generating functional of the correlators of j's 



J D} exp {- / d'xd^yj,{^r-^j^j,{y) + J d^x (-2C + ^ + J,j,)] } 



/Djexp{- 



^ / d^xd^yj,i^r-^;^Uy) + J dH (-2C + ^ 



)]} 



= exp 



J d^x J rfVM(x)/C(x-y)J^(y) 



;i83) 



where /C(x) = ^S^id'^ - "^^)^^- Next, since due to Eq. ([T8|) = (9^5s^^ = 0, the Hodge 
decomposition theorem requires that should have the form = e^y\dyip\. Then owing to the 
same theorem, the coupling / d^xJ^j^ is nonvanishing only provided that can be represented 
as e^uxdylx- This coupling then takes the form / d^xI^T^^^y, where T^^^-x) = d^d^ — S^^d^'^. On 
the other hand, once being substituted into the R.H.S. of Eq. (|183|) , such a representation for 
yields 



Z[3] = exp 



d^x J dV/.(x)/.(y)T^.(x)/C(x-y) 



Thus, varying Z[3] twice w.r.t. I, we get 

T^^(x)TAp(y) ((^^(x)(^p(y)) = -2T^a(x)/C(x - y). 

Due to the rotation- and translation invariance of space, it is further natural to write down for 
the correlator {{py{x.){pp{y)) the following ansatz: 5i,p/(x — y). This yields 



|x-y| 



;i84) 



The desired correlator of j's reads 



{3M)jy{y)) = ^f^aisSupad^dy ((/?^(x)yp^(y)) . 
Taking into account the above obtained results, we eventually have 



(j^(x)j.(0))=T^.(x)/(x) 



m 



D 



fc2 



m 



d'y- 



-M\y\ 



y\\y\ 



where in a derivation of the last equality the definition of /C and Eq. ( p.84| ) have been employed. 
The remained integral can most easily be calculated by dividing the integration region over |y| 
into two parts, [0, |x|], [|x|,+oo), and expanding j^^z^ in Legendre polynomials P„'s on both of 
them. Then, the integration over the azimuthal angle singles out from the whole series only the 



72 



zeroth term, / Pn{cos9)dcos9 = 26no- After that, the integration over |y| at both intervals is 



straightforward and yields 



-M\y\ 



y\\y\ 



M2| 



(1- 



-Af|x| 



Note that an alternative method of calculation of this integral (which is the 3D analogue of the 
respective 4D calculation, to be discussed below) can be found in Ref. | |133[| . Finally, we arrive at 
the following expression for the bilocal correlator of vortex dipoles: 



0-^(x)j,(0)) 



1 ( V 



2 I 2 ■ 

m + rrij^e 



-Af|xp 



vr \A7ir]M J 

Clearly, when m/j — 0, this correlator vanishes, which means that it is the Debye screening in the 
ensemble of vortex dipoles, which is responsible for their correlations. 

Let us now proceed to the 4D-case of small electric vortex loops in the DAHM. This model is of 
more interest for us, since it is the 4D-case, where we have calculated field strength correlators in 
the previous Section, and thus the final aim of further investigations will be to find modifications 
of these correlators due to the interactions in the ensemble of vortex loops. To study the grand 
canonical ensemble of such loops, it is necessary to replace S^j^ in Eq. (|172|) by the following 
expression 



X] 



N 
a=l 



da,,{x%0)S{x-x%0), 



;i85) 



where below in this Section we shall set ^ G [0,1] x [0,1]. After that, the summation over the 
grand canonical ensemble of vortex loops is similar to the summation over the ensemble of vortex 
dipoles in the dual Ginzburg-Landau theory. Referring the reader for the details to Appendix D, 
we shall present here the result of this procedure, which has the form 



Z 



Dh 



exp 



d^x 



Ut]' 



tt2 I „2 l2 



2( COS 



\h 



liv I 



A2 



;i86) 



Here \h, 



and the fugacity C (Boltzmann factor of a single vortex loop) has now the di- 



mension (mass)^. We have also introduced a new UV momentum cutoff A equal to \J with L and 
a denoting the characteristic distances between vortex loops and their typical sizes, respectively. 
The square of the full mass of the field /i^i, following from Eq. ( |186| ) reads = + m|) = Q^rf. 



Here, 



2(yfm?7 is the usual Higgs contribution, mn 



-^^^ is the Debye contribution, and 



Q = 2\J + is the full magnetic charge. 

The representation for the partition function (|186|) in terms of the vortex loops can be obtained 
by virtue of the following equality 



exp 



d X 



127]' 



= J DSf^^expl^- {nT]Y J d'^x J d'^yS^^{x)D^^}{x - y)S^^{y) + m J 



d'^Th ^ 



(187) 
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It can be derived by taking into account that d^h^i, = 0. Indeed, owing to the Hodge decomposition 
theorem, the Kalb-Ramond field can always be represented as follows: h^i, = d^ip^ — d^ipfj, + 
£tiu\pd\ipp. Clearly, the field y?^ decouples not only from E^^^ (due to the conservation of the latter 
one), but also from ip^. The yj^-field thus yields only an inessential determinant factor, which is 
not of our interest. Therefore this field can be disregarded, which proves the above statement. 

Substituting Eq. ( |187D into Eq. ( |186| ), we can integrate the field h^i, out, which yields the 
desired representation for the partition function ( p.86| ). 



with the effective potential V given by Eq. 

Similarly to the 3D case, correlation functions of S'^j^'s, calculated by virtue of the partition 
function ( |188| ), are nothing else, but the correlation functions of vortex loops in the gas. This can 
be seen by mentioning that if we insert into the partition function 2 = J Dx^{^)Dh^ye~^ with 
the action 5* given by Eq. (|172|) and S^j, replaced by S^^^ the following unity 



DS^yDLu exp 



v^gas 



;i89) 



and integrate out all the fields except S^u, the result will coincide with Eq. ( p.88| ). Such correlation 
functions of the vortex loops can be most easily calculated in the low-energy limit, 15^,^1 -C (, 
by restricting oneself to the vicinity of the minimum of the real branch of the potential (|166D , 
where it reduces just to a quadratic functional of S^u. To calculate these correlation functions, let 
us derive an expression for the respective generating functional. It reads 



flUj 



J DS,, exp {- [(vrr^)^ / d'x J d'yS,,{x)D(^\x - y)S,,{y) + J d'x {^Sf^, + J,,S,,)] } 



/ DS,, exp {- [{tttj)^ j d^x J d^yS,,{x)D^l^ {x - y)S,,{y) + ^ / J } 



exp 



where 



d'^x j d^yJf,u{x)G{x - y)Jf,u{y) 



(190) 



C 



{d'-m')D^}{x) 



(191) 



Next, since d^T^f^ 



7r2A4 

0, Eq. (|189|) requires that d^Sfj^u 
theorem leads to the following representation for 5*^,^: 
same theorem, the coupling / d'^xJ^^S^u will be nonvanishing only provided that J^y = e^y\pd\Ip. 
This coupling then reads 2 / d'^xIpT^^ip^. On the other hand, substituting the above representation 
for into the R.H.S. of Eq. ( |19CI|) , we have 



as well. Then, the Hodge decomposition 
Sfj.u = £pu\pd\(pp. Owing to this fact and the 



exp 



2 d^x d'^yIp{x)Iy{y)Tp^{x)Q{x - y) 



Thus, varying Z[J^y\ twice w.r.t. J^^ and setting then J^^ = 0, we get 
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T^^{x)Txp{:y) {ipu{x)ipp{y)) = ~Tpx{x)g{x - y). 

Again, due to the rotation- and translation invariance of space-time, it is natural to seek for 
{(p^{x)(pp{y)) in the form of the following ansatz: S^pg{x — y). This yields the equation d'^g = Q, 
whose solution reads 



9{x) - 



^2A4-(^^'-^') / d'yDlt\x-y)D'i!{y), 
where Dq^\x) = D^^{x) at m = 0, i.e., it is just ^^1^2 ■ The last integral can obviously be rewritten 



as 



J d'zDi'\z)Dtf{z-x). (192) 

As we will see below, it will be necessary to know the more general expression, namely that for 
the integral 



d'zD^^\z)D^^\z-x). (193) 
Its calculation is outlined in Appendix E, and the result reads 

^{Dl^\x)-Dfhx)). (194) 



Note that as it obviously follows from Eq. ( p,93| ), since this result depends on x only as it 
should be symmetric w.r.t. interchange m M. One can see that this really holds: Eq. (|194|) is 
invariant under this interchange, since during it mj^ changes its sign. 

Setting now in Eq. (|194|) m = 0, we get ^ (^Dq'\x) — D^^(x)j, which yields for the desired 
integral (|192|) the same result with the substitution mo ^ M. Thus, the final expression for the 



function g reads 



g[x 



id' - m') {off (x) - Di^\x)) . (195) 



The desired correlator of Sp,^'s has the form 



(7rMA2^2 
.as the 

{Sp^{x)Sxpiy)) = Sp^acpexp^adldy {ipp{x)ip^{y)) 

and therefore 

{Spy{x)Sxp{Q)) = -epyap£\p^pdld^g{x) = 

= {hJpp - S^pSpx) G{x) + {6pxdpd^ + d^pdpdx - Sppdxd^ - dxudpdp) g{x), (196) 

where it has been used that d'^g{x) = Q{x). 

This result can immediately be applied to the calculation of the string contribution to the 
bifocal cumulant (PS]). Indeed, applying to the average on the R.H.S. of Eq. ( PO]) the cumulant 
expansion in the bilocal approximation, one gets: 



26 



Note that this result can also be obtained directly by making use of the method presented in Appendix E, 



which was done in Ref. |130| 
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Z ^ expj — y" d^x J d^yD^^i^x — y) 



(4vrr^)'Sf,(x)Sj,(y) + ij^f(a:)jJ(y) 



+ 



Comparing this expression with Eq. (pT]), we see that owing to Eq. ( |196| ), the additional string 
contribution to the cumulant (|88|) has the form 

A ((^.(x)/ap(2/))),^,^.m = {AT^g^r,^)' f d^zd\D^^\x - z)Dt\y -u)x 



X {{5^xSup - Sf,p5ux) G{z -u)+ [S^pdld^ + Suxd^d^ - S^xd^dl - S^pd'^d^}^ g{z - . 
Comparing this equation further with Eq. (|88D and taking into account that 

{x - y),V, {{x - yf) = -^-d;G {{x - y)') , 

where the function G is defined by Eq. ( |36D with the replacement Di — > Vi, we arrive at the 
following system of equations, which determine string contributions to the dunctions T> and G: 



AV ((x - yf) = (4W) Y d'zd\Dl^\x - z)D^^\y - u)Q{z - u), 

AG({x-yf) = {Sng^v'f J d'zd%D^^\x - z)Dt\y - u)g{z - u). 
Inserting now Eq. ( p.91|) into Eq. (|197| ), we get 



(197) 



(19^ 



AP Ux - y)- 



By virtue of the result of Appendix E, we have 



j d%Dt\y-u)D^S{^-u). 



AV [x^) 



m 



47r^ 



M 



1 VI TTl 

—Ki{M\x\) - —Ki{m\x\) 



x\ 



Adding this result to Eq. (|9TD, we finally obtain for the finction V the following full result: 



P^-" x^ 



m'^M Ki{M\x\) 

47^2 



(199) 



Analogously, inserting Eq. ( |195|) into Eq. ( [198| ), we have 



AG({x-yf)=C 



^ Am 

or further by virtue of Appendix E, 



,2\ 2 



d\Dl^\y-u) 



Di'\x - u) - D^hx - u] 
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Together with Eq. (P^D, this yields the following full result for the function Vi. 



full 



m 



D 



+ 



m 



IT 



2M2|x|4 27r2Ma;2 



™)+|(,,„(M|.|) + A,(MW)) 



(200) 



It is worth noting that the functions AD and ADi contained the terms exactly equal to 
Eqs. (|91|) and (p2|), respectively, but with the opposite sign, which just canceled out in the full 
functions ( |199| ) and (|200|) . We also see that, as it should be, the functions ( |199D and (|200|) go over 
into Eqs. (0) and (|92D, respectively, when niD 0, i.e., when one neglects the effect of screening 
in the ensemble of vortex loops. An obvious important consequence of the obtained Eqs. (|199|) 
and ( |200[ ) is that the correlation length of the vacuum, Tg, becomes modified from 1/m (according 
to Eqs. (|9TD and (0)) to 1/M. (It is worth emphasizing once more that this effect is just due to 
the Debye screening of magnetic charge of the dual vector boson in the ensemble of electrically 
charged vortex loops, which makes this particle more heavy, namely enlarges its mass from m to 
M.) Indeed, it is straightforward to see that at |x| ^ 



V 



full 



{mMf 



-M\x\ 



AV2n2 {M\x\y- 



and 



full 



m 



D 



7r2M2 



+ 



{mMf 



-M\x\ 



X\ 



2V27ri {M\x\)^ 

It is also remarkable that the leading term of the large-distance asymptotics of the function 
D^"^^ is a pure power-like one, rather than that of the function T>i, given by Eq. (p^). Another 
nontrivial result is that the screening does not change the short- distance asymptotic behaviours of 
the functions (0) and (p^), i.e., the short- distance asymptotics of the functions ( |199| ) and (|200| ) 
are given by Eqs. (|95D and (|9^), respectively. 

Let us now turn ourselves to the case of the Abelian-projected S'?7(3)-gluodynamics. In the case 
when Abrikosov-Nielsen-Olesen type electric strings in this theory are considered as noninteracting 
objects, the expression for the partition function is given by Eq. ( |108| ), which can be shortly written 



as 



Z = j Dxl{i)Dxl{i)x 



X exp < -Qrari 



d'xdS T}(x)T}(y) + T}(x)T?(y) + T?(x)T?(y) 



Ki{niB\x - y\) 



\x-y\ 



(201) 

In order to proceed from the individual strings to the grand canonical ensemble of interacting 
vortex loops, one should replace S^^(x), where from now on a = 1, 2, in Eq. (|201|) by 



N 

^;rix) = E < / d'y.u (xtio) 5 (x - 4(0) 
k=i •' 



Here, n^'s stand for winding numbers, which we shall again set to be equal ±1. Performing such 
a replacement, one can see the crucial difference of the grand canonical ensemble of small vortex 
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loops in the model under study from that in the Abelian-projected SU (2)-gluodynamics, studied 
above. Namely, the system has now the form of two interacting gases consisting of the vortex 
loops of two kinds, while in the S'?7(2)-case the gas was built out of vortex loops of the only one 
kind. 

Analogously to that case, we shall treat such a grand canonical ensemble of vortex loops in 
the dilute gas approximation. According to it, characteristic sizes of loops are much smaller 
than characteristic distances between them, which in particular means that the vortex loops are 
short living objects. Then the summation over this grand canonical ensemble can be most easily 
performed by inserting the unity 

1 = / Ds;j {s;^ - s^r) (202) 

into the R.H.S. of Eq. ( |201|) (with S^^^ replaced by S^^*^*^) and representing the 5-functions as the 
integrals over Lagrange multipliers. Then, the contribution of vortex loops of each kind to the 
grand canonical ensemble takes the following form 



Z 



y a gas 




X exp <; -gniv\l 2 / d'^xd^ Sl^{x)Sl^{y) + Sl^{x)Sl^{y) + Sl^{x)Sl^{y) 



Ki{mB\x - y\) 
\x - y\ 



i J d'x\% {S% - 



gas 



(203) 



After that, the desired summation is straightforward, since it technically parallels the one of 
Abelian-projected 5'f/(2)-gluodynamics. We have 



oo / N . „ 



X 



X 



nl=±l 



N 

Y.nl da,^ (4(0) (x^(0) 



k=l 



X {the same term with the replacement of indices 1^2} 



exp I 2C j d^y 



cos 



A2 



+ COS 



A2 



(204) 



Here, the world-sheet coordinate of the k-th vortex loop of the a-th type ^ x1{C,) has been decom- 
posed as xl{^) = y% + zl{^), where the vector y% = S d^^x%{^) describes the position of the vortex 
loop, whereas the vector z'^{S,) describes its shape. Next, on the L.H.S. of Eq. ( p04| ), /i [z^] again 
stands for a certain rotation- and translation invariant measure of integration over the shapes of 
vortex loops, and C, oc e"'^'' denotes the fugacity (Boltzmann factor of a single vortex loop 0) of 

^^For brevity, we omit the Lorentz index. 

^*It is natural to assume that the vortex loops of different kinds have the same fugacity, since different 6'^'"S 's 
enter the initial partition function (104) in the same way. 
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dimension (mass)^ with Sq being the action of a single loop. In Eq. ( p^04| ) , we have also introduced 

the UV momentum cutoff A = -y/^ where a is a typical size of the vortex loop, and L 

is a typical distance between loops, so that in the dilute gas approximation under study a <ti L. 

3), we have denoted 



Finally in Eq. 



Next, it is possible to integrate out the Lagrange multipliers by solving the saddle-point equa- 
tion following from Eqs. ( P03| ) and (|204|) , 



sm 




' 2C 



After that, we arrive at the following representation for the partition function of the grand canon- 
ical ensemble 



^grand — ^ ^iv ^-^P 



9rari\\- I dVy \sU^)Sl^iy) + SUx)SUy) + SU^)Sl{y) 



X 



Ki{mB\x - y\) 
\x - y\ 



+ V 



+ V 



(205) 



which owing to Eq. (202) is natural to be referred to as the representation in terms of the vortex 



is given by Eq. (|166|) . 



loops. In Eq. (|205|) , the effective potential of vortex loops V 

It is further instructive to illustrate the difference of such a partition function of two interacting 
gases of vortex loops from the case of Abehan- projected SU (2)-gluodynamics by studying a related 
representation in terms of a certain effective sine-Gordon theory. This can be done by introducing 
the new integration variables S^^^, = ^ (^Sj^^ + S"^^ j and S"^^ = | (^Sj^^ — Sf^Jj , which diagonalize 
the quadratic form in square brackets on the R.H.S. of Eq. (|20^). Then Eqs. (|203|) and (|204|) yield 



Z, 



grand 



+2C J d^x 



+ COS 




(206) 



A^j,. The partition function of 



where we have denoted h^^ = ^ (a^^ + A^^) and = A^^ 
the desired sine-Gordon theory can be obtained from Eq. ( p06| ) by making use of the following 
equality 




(207) 



(c/. the R.H.S. with the quadratic part of the action of the Kalb-Ramond fields on the R.H.S. 
of Eq. ( |1 201 ) with the fields gauged away). It can easily be proved by noting that due to the 
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Hodge decomposition theorem and the equation d^S^^^, = (which follows from Eq. ( |202D and 
conservation of S?;P^), duh'f,,, = 0. Substituting further Eq. (|207|) into Eq. ( P06| ) and performing 



the rescaling h^^^, we arrive at the following representation for the partition function of the 

grand canonical ensemble of vortex loops in terms of the local sine-Gordon theory, equivalent to 
the nonlocal theory 



^grand = J Dk^^^SXpl- J 



d X 



-2C 



cos 



WV2 



+ cos 



2 3 2 
+ 2^rn 



(208) 



As we now see, an essential property of the obtained sine-Gordon theory, which distinguishes it 
from an analogous theory describing the grand canonical ensemble of vortex loops in the Abelian- 
projected S'L'^(2)-gluodynamics, is the presence of two interacting Kalb-Ramond fields, while in 
the S'?7(2)-case there was only one self-interacting field. Notice that upon the expansion of the 
cosines on the R.H.S. of Eq. ( p08|) , it is straightforward to see that only the interaction terms of 

the type (hl^y^ {^%^ survive. In another words, despite of the mixing of the Kalb-Ramond 
fields in the arguments of the cosines, no terms linear in any of these fields appear in the action. In 
particular, the full masses of both Kalb-Ramond fields. Mi and M2, can be read off from Eq. (|208D 



by expanding the cosines up to the quadratic terms. The result reads = m\ + = Qa^^, 



3C, ma 



^a/C are the Debye masses, and we have introduced the magnetic 
Q2 



where mi 

charges Q, = Jogl + l^, Q, = ^6gl + ^. 

Equation (|205|) can now be used for the evaluation of correlators of vortex loops, which due 
to Eq. ( |202| ), are nothing else but the correlators of S^^^s. Those are calculable in the low-energy 

limit, S'l^j, -C C; by considering the real branch of the potential (|166|) , i.e., extracting from 
the sum over branches of this potential the term with n = 0. This branch has a simple parabolic 
form, and in the vicinity of its minimum (corresponding to the low-energy limit) the generating 



functional for correlators of S'^j.'s reads 



Z 



ja 



DS;,exp 



9mV 




Ki{mB\x - y\ 
\x-y\ 



d X 



s}JA + sLj- 



where J^^ is a source of S^^,, and = J^^ ± J^^,. Such two Gaussian integrals can be calculated 
by virtue of the following equality 



DSa^ exp 



Ki{mB\x - y\) 



\x - y\ 



/A^ 

S^,u{y) + I d^xi T^Sl^ + J^yS, 



2C 



= exp 

and the result reads 



M2C 
87r2A4 



d^xd%Jf,i,{x)Jf,u{y) {dl 



Ki{M2\x-y\ 
\x - y\ 
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7 



exp 



rf^x j dS [j^u{x)J^^{y)Gi{x -y) + ^^^.(a;) J^^(i/)^2(a; - y)] | 



where Qa{x) = 2^ (9^ — m\) D^^}^{x). Owing to the conservation of and the Hodge decom- 
position theorem, we again have 5"^^, = Sfj^^xpdxfp and J^^, = e^^xpdxlp, which yields 



Z 



ja 



exp< —2 / d'^x / d'^y 



i;{xmy)T,,{x) {G^ix -y) + G^ix - y)) + 



+2ll{x)ll{y)Tp,ix) (Giix -y)- G^ix - y)) 



On the other hand, the couphng / d^xJ^^S^^ takes the form 2 / d'^xI^T^^ipl, and varying Z 
twice w.r.t. /^'s we arrive at the following system of equations: 



ja 



Tp^{x)Txp{y) (ip^{x)ip (y)) = Tp^{x)Txp{y) {ip^{x)(p Jy)) = -T^xix) {Giix - y) + G2{,x - y)) , 



Tp^{x)Txp{y) [^l{x)(pl{y)) = -Tpx{x) {Gi{x - y) - G2{x - y)) . 
Adopting for the correlators of <y9^'s the following ansdtze, 

{vl{x)^l{o)) = {vl{x)v%o)) = {v>lix)v%y)) = s.pUx), 

we get: 

c 



2A4 



ml 



Ml - - .... 

This result makes the choice of notations ^^f±{x)" quite natural. Finally, the correlators of vortex 
loops read 

SU^)Sl{0)) = (sl^{x)Sl{0)) = {6,xS,p - 6M iSiix) + G2{x)) + 



+ {Spxdpd^ + 6^pdpdx - 6^,pdxd^ - S^xdf^dp) f+{x) 



SU^)Slp{0)) = {6.xS,p - 6M iOiix) - G2{x)) + 



+ {5p\dpdy + 5ypdpdx - 5ppdxdy - S^xd^^dp) f-{x). 

This result can now immediately be applied to the calculation of the modification of the 
correlator ( |110| ) due to the screening in the ensemble of vortex loops. Indeed, applying to the 
average over world-sheets, standing on the R.H.S. of Eq. ( [L13| ) the cumulant expansion in the 
bilocal approximation, we have due to Eq. ( |111| ): 
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A {{fU^)fliy))) = -247rW5i^^^^ / d'zd\DiZ{x - z)Dl^^{y - u) {S%{z)S\,{u)) . 

Taking further into account the equahties 

{slA^)S%{y)) = {sl^{x)Sl{y)) = {sl^{x)Sl{y)) , 
and the facts that for every c, s'f^s'f^ = 2, s^i^s''2^ + s^f^s^^^ + s'^^s^^^ = —1, we can write 

This leads to: 

AD = 487rW/ d^zd^D^^Kx - z)D^J,liy - u)g2{z - u), 

AG ((x - yf) = gevrW/ d^zd^uD^^Hx - z)D^i\{y - u) {f^{z -u)-f^{z- u)) , 

where G is given by Eq. ( ^61) with the replacement Di ^ Di. Carrying now the integrals analo- 
gously to how it was done above in the 5'[/(2)-case, we get 

Ab {x')=ml {d^SM - D(^l{x)) , 




Together with the old expressions for the functions D and Di (given by Eqs. ( pi]) and (p2D, 
respectively, with m ruB), which did not account for the screening effect in the ensemble of 
vortex loops 0, we finally obtain that D^"^^ and D^^^^ are given by Eqs. ( |199|) and ( pOOp , respectively, 
with the replacements m ^ rriB, itld m2, and M —>■ M2. Therefore the whole discussion, 
following after Eq. ( |20CI|) , remains the same modulo these replacements. In particular, when m2 
vanishes, i.e., one disregards the effect of screening, the old expressions for the functions D and 
Di are recovered. 

5 Conclusion and Outlook 

The main problem addressed in the present review was an attempt of an analytical description 
of confinement in QCD and other gauge theories. As a guiding principle for our investigations 
served the so-called Wilson's picture of confinement, according to which this phenomenon can be 
described in terms of some effective theory of strings, joining coloured objects to each other and 
preventing them from moving apart to macroscopic distances. In this review, we have proceeded 
with a derivation of such string theories corresponding to various gauge ones, including QCD, 
i.e., with the solution of the problem of string representation of gauge theories. We have started 

It is remarkable that these expressions again become exactly canceled by the corresponding terms in AD and 
ADi. 
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our analysis with the nonlocal string effective action, arising within the so-called Method of Field 
Correlators of QCD, where the interaction between the string world-sheet elements is mediated by 
the phenomenological background gluon propagator. By performing the derivative expansion of 
this action, we have derived the first few terms of the string Lagrangian. The first two nontrivial of 
them turned out to be the Nambu-Goto and rigidity terms with the coupling constants expressed 
completely via the gluonic condensate and correlation length of the QCD vacuum. The signs of 
these constants ensure the stabihty of strings in the so-obtained effective string theory. After 
that, we have investigated the problem of crumpling of the string world-sheet by a derivation 
of the topological string term in the instanton gas model of the gluodynamics vacuum. Next, 
by making use of perturbation theory in the nonpcrturbative QCD vacuum, we have calculated 
perturbative corrections to the obtained string effective action. Those led to a new form of the 
nonlocal string effective action with the propagator between the elements of the world-sheet being 
the one of a perturbative gluon in the confining background. By the derivative expansion of this 
action, we got a correction to the rigidity term coupling constant, whereas the string tension of 
the Nambu-Goto term occurs to get no corrections due to perturbative gluonic exchanges. Finally, 
we have derived the Hamiltonian of QCD string with spinless quarks at the ends, associated with 
the obtained string effective action including the rigidity term. In the particular case of vanishing 
orbital momentum of the system, this Hamiltonian reduces to that of the so-called relativistic 
quark model, albeit with some modifications due to the rigidity term, which might have some 
influence to the dynamics of the QCD string with quarks. All these topics have been elaborated 
on in Section 2, and form the essence of the string representation of QCD within the Method of 
Field Correlators. 

In Section 3, we have addressed the problem of string representation of Abelian-projected 
theories. In this way, we have started with the string representation for the partition function 
of the simplest model of this kind, namely the Abelian-projected S'f/(2)-QCD, which is argued 
to be the dual Abclian Higgs Model with external electrically charged particles. The advantage 
of this approach to the string representation of QCD w.r.t. the one based on the Method of 
Field Correlators is a possibility to get an integration over string world-sheets, resulting from the 
integration over the singular part of the phase of the Higgs field. After the string representation for 
the partition function in the London limit, we have derived such a representation for the bifocal 
cumulant of the field strength tensors. iFiom this, in the confining regime for a test external 
quark, we got the expressions for the two functions, which parametrize the bifocal cumulant 
within the Method of Field Correlators. The obtained results demonstrate that the large-distance 
asymptotic behaviour of the bifocal field strength cumulant matches the one of the corresponding 
gauge- invariant cumulant in QCD, predicted by the Method of Field Correlators and measured in 
the lattice experiments. In particular, the role of the correlation length of the vacuum, introduced 
in the Method of Field Correlators as a phenomenological input, turned out to be played in the 
Abelian-projected SU (2)-gluodynamics by the inverse mass of the dual gauge boson. These results 
support the method of Abelian projection on the one hand and give a new field-theoretical status 
to the Method of Field Correlators on the other hand. After that, we have discussed the other 
fundamental nonpcrturbative phenomenon of QCD, the chiral symmetry breaking, from the point 
of view of the Abelian-projected theories and the Method of Field Correlators. In particular, in 
the SU (2)-case under study, we have quoted the derivation of the relation between the chiral and 
gluonic condensates, following from these two approaches. 

Next, by making use of the Abelian projection method, we have addressed the problem of 
string representation of the SU (3)-gluodynamics. Namely, we have casted the related dual model. 
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containing three types of the dual Higgs fields, into the string form. Consequently, the latter 
one turned out to contain three types of strings, among which, however, only two were actually 
independent. As a result, we have found that both the ensemble of strings as a whole and 
individual strings display confining properties in the sense that these two distinct types of strings 
(self)intcract via the exchanges of the massive dual gauge bosons. We have further also derived 
string representation for the Abelian-projected 5't/(3)-QCD with quarks, from which there was 
again obtained the related string representation for the bilocal cumulants of field strength tensors. 
In the confining regime for a test quark, it coincides for diagonal cumulants with the one of the 
SU{2)-case modulo the modification of mass of the dual vector boson. This means that only w.r.t. 
those cumulants, which are built out of gluonic fields referring to the same generator of the Cartan 
subalgebra the vacuum of the model under study exhibits a nonvanishing correlation length. The 
latter one again turned out to be equal to the inverse mass of the dual gauge bosons, which they 
acquire by means of the Higgs mechanism. 

In conclusion of this topic, we have derived another useful representation for the partition 
functions of the Abelian-projected theories in the form of the integral over monopole currents. 
Besides the part quadratic in these currents, which represents the Biot-Savart interaction and 
the kinetic energy of monopole Cooper pairs, the obtained Lagrangians contain also the terms 
describing the long-range interaction of magnetic currents with the string world-sheets. If one 
localizes the (quantum) monopole currents along the classical trajectories, these terms take the 
form of the Gauss linking numbers between the currents and world-sheets. This means that the 
latter objects may be viewed as solenoids, which scatter the former ones (c/. the Aharonov-Bohm 
effect). 

In Section 4, we have studied another model, allowing for an analytical description of con- 
finement, which is the compact QED. In this way, we have successively investigated the 3D- and 
4D-cases and in both of them derived the string representations for the Wilson loop of an external 
electrically charged test particle. The essence of these representations is a certain mechanism 
realizing the independence of the Wilson loop of the shape of an arbitrary surface bounded by 
its contour. In this way, it has been argued that this mechanism is based on the summation over 
the branches of a certain multivalued effective potential, which in the 3D-case is the potential 
of the monopole densities, whereas in the 4D-case it is the potential of the monopole currents. 
In the former case, we have established a correspondence of this approach to another recently 
found one, the so-called confining string theory. After that, in both cases, we have calculated the 
bilocal cumulants of the field strength tensors in the low-energy limit of the model under study. 
Those also turned out to be in line with the general concepts of the Method of Field Correlators 
and therefore match the corresponding results known from the lattice measurements in QCD and 
found analytically for the effective Abelian-projected theories in Section 3. Note only that con- 
trary to the Abelian-projected theories, in compact QED monopoles are not condensed (in the 
sense that the dual Higgs field does not exist), but form a gas, and consequently the dual gauge 
field acquires a nonvanishing mass due to the Debye screening in such a gas, rather than due to 
the Higgs mechanism. Correspondingly, the correlation length of the vacuum in compact QED is 
equal to the inverse Debye mass of the dual gauge boson. 

On the basis of the found correspondence between string representations of various gauge 
theories and the above discussed similarities in the large-distance asymptotic behaviours of the 
bilocal cumulants in QCD, Abelian-projected theories, and (low-energy limit of) compact QED, 
we have elaborated on a unified method of description of the string world-sheet excitations in 
the respective string theories. This method, which employed the techniques of nonlinear sigma 
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models, enabled us to derive the effective action, quadratic in the world-sheet ffuctuations. 

Finally, we have applied the methods of summation over the grand canonical ensemble of 
monopoles in compact QED to the summation over the grand canonical ensemble of topological 
defects (vortex dipoles built out of electric Abrikosov vortices in 3D and vortex loops built out of 
electric Nielsen-Olesen strings in 4D) in Abelian-projected theories. In this way, we have evaluated 
correlation functions of these defects in the low-energy limit. This enabled us to improve on the 
calculation of the bifocal field strength correlators in Abelian-projected theories by a derivation 
of the corresponding string contributions. In this way, it has been demonstrated that the true 
correlation length of the vacuum in these theories is given not simply by the inverse mass of the 
dual vector bosons, acquired by them due to the Higgs mechanism, but rather by their full mass, 
which accounts also for screening in the gas of topological defects. Besides that, it turned out 
that the effect of screening leads also to the appearance of a power-like term in the large-distance 
asymptotics of one of the two functions, which parametrize the bifocal cumulant. 

To summarize, the correspondence between various field-theoretical models, considered in the 
present review, as well as the approaches to their string representation may be established by 
the Table 1, presented at the end of this Section. In conclusion, the proposed nonperturbative 
techniques provide us with some new information on the mechanisms of confinement in QCD 
and other gauge theories and shed some light on the structure of the vacua of these theories. 
They also demonstrate the relevance of the Method of Field Correlators to the description of this 
phenomenon and yield several prescriptions for the construction of the adequate string theories 
from the corresponding gauge ones. 

Further investigations of the problems addressed in the present review are planned to follow in 
at least two directions. First of them is a derivation of the string representation for the last up to 
now known model allowing for an analytical description of confinement, which is the low-energy 
SU{2) Af = 2 supersymmetric Yang-Mills theory (the so-called Seiberg-Witten theory) |3^. The 
second line of investigations is devoted to a better understanding of interrelation between confine- 
ment and chiral symmetry breaking. In this way, it is planned to develop further the approach 
proposed in Ref. |^ by virtue of bosonization of the equations obtained there. This turned out 
to be a rather hard topic due to the non-translation-invariant character of the corresponding in- 
teraction kernel. However, it is this kernel, which is responsible for the confining effects in the 
system, and therefore it should not probably be reduced to some more simple translation-invariant 
one. Investigations of this problem together with an application of other well elaborated methods 
known from the theory of NJL models might finally enable one to understand completely the 
interrelations between the phenomena of confinement and chiral symmetry breaking in QCD. 
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Table 1. Correspondence between various field-theoretical models and approaches to 
their string representation 



Model 


QCD within MFC 


Abelian-pro j ected 
theories 


Compact QED 


Mechanism of the 
string representation 


No integral over 
string world-sheets. 
String effective ac- 
tion is defined only 

W.r.t. Smin. 




S- independence of 
{W{C)) is realized 
by the summation 
over branches in V[p] 
in 3D or V\j^,] in 4D 


Mechanism of the 

mass goiioratioii 


Due to stochastic 

l)ackgrouii(l fields 


Higgs mechanism 


Debye screening in 

the iii()ii()pok> gas 


Type of propagator 
between the elements 
of the world-sheet (s) 


Nonperturbative 
gluon propagator 
(D-function) or 
propagator of per- 
turbative gluon in 
the nonperturbative 
background 


Propagators of the 
Kalb-Ramond fields 


Propagator of the 
Kalb-Ramond field 


Parameter of the ex- 
pansion of the result- 
ing nonlocal interac- 
tion between the el- 
ements of the world- 
sheet (s) 


Correlation length of 
the QCD vacuum, Tg 


Inverse mass of the 
dual gauge bosons, 
which they acquire 
both due the Higgs 
mechanism and due 
to the Debye screen- 
ing in the ensemble of 
topological defects 


Inverse Dcbyc mass 
of the dual gauge bo- 
son 
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7 Appendices 



A. Details of a Derivation of the Hamiltonian of the Straight-Line QCD 
String with Quarks 

Assuming that a meson as a whole moves with a constant speed (which is true for a free meson), 
i.e., R = 0, and bringing together quark kinetic terms ( |62D and the pure string action (0), we 
arrive at the following action of the QCD string with quarks 



i \2/ii 2/i2 2 2 2 



+ /i2 + y I 





+ I /ii(l - Ci) - /i2Ci + / d(3{(3 - Ci)u] (Rr) - J df3ur] (i?r) + J d(3{Ci - (3)w{rr) + 

/ 



1 

+ — 

"0 



C2(/ii(Ci - 1) + /i2Ci)r' - C2(/Wi + f^2) {Rr) + I dPuUiCi - Py - C2 (Rr) + C2Virr) + 



+ l{P-Cir [t,r]Vii?V-i (i?r)' + (/?-Ci) ((ri?) (rr) - (ff) (i?r)) + 



7]'^ ] [r, r]^ + 



+^( (/^ ~ Ci) ((rr)(fr) — (fr)r^^ + (^i?r^ — (^Rrj (rr) 



(Al) 



where we have performed a rescaling — > ay -^z^, — > ay -^z^, v — > ^^z^. Integrating then 
over 7], one gets in the zeroth order in — 



_ (rr) / /ii 

'/extr. — — \ P 



l^l + At2 



which together with the condition Rf = yields 



cr 



extr. 



/il + / d/S/Sl' 




r 

1 ' S2 

A*! + At2 + / d/3iy 





extr. 



(rr) A'l+A'a 



i^Jd(3u-} d(3f3u 





1 



/Ul + Ai2 + / 




Finally, in order to obtain the desired Hamiltonian, we shall perform the usual canonical 
transformation from R to the total momentum P, which in the Minkowski space-time reads 

J DRexp ij L (r, . . .) rfr = J DROP exp i J (PR -H{P,.. .)) dr 
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with if (P, . . .) = PR — L (ii, • • •) , and choose the meson rest frame as 

dL f R, . . .) 

P = ^ = 0. 

dR 

After performing the transformation from r to p we arrive at the Hamiltonian (35) of the main 
text with the coefficient functions a^, which read as follows 

CT^ }d(3 Ji ■ 

fli = — / — ,0,2 = 3-B - B, 
2 J u a 
^ 



as 



1 



2(/ii + /i2)(/Ul + 1^2 + J^o) 



(/il + /U2)(/il + /i2 + l^o) 



z/i (/ii + /ia) (z^i - 2/i2) - AiiZ/o(/ii + 2/i2) 



(34 = -^-B — B, 



a5 = 1^2 + 



+ 2^11^0 - 2fi2^i 1 



/Xl + /i2 + Uo 



+ 



/il + 1^2 



1 / il^l^O - MlJ'l + /i2))^(3z/o + 2(/ii + /i2)) 



/il + 1^2 



(/il + /i2 + t'o)^ 



+/il(/ilZ/o - 2z/i(/ii + /i2)) - 



/i?t'0 



/il + /i2 + 1^0 



and 



5 = 



Z/i(/ii + /i2)(z^l - 2/i2) + filMfJ'l + 2/^2) 



(/il +/i2)(/il +/i2 + t'o) 

Notice that during the derivation of H^^\ we have chosen the origin at the centre of masses of 

T . . . 2 

the initial state, so that Rr ^ 1, and the term / dru^ (Ri") on the R.H.S. of Eq. (A.l) has 
been disregarded. 



B. Path-Integral Duality Transformation 

In this Appendix, we shall outline some details of a derivation of Eq. (Q) from the main text. 
Firstly, one can linearize the term ^ {9fiO — 2gmB^)'^ in the exponent on the R.H.S. of Eq. (|82D 
and carry out the integral over as follows 



j DO^'^^- exp / d^x {d^e - 2(/„i?^)'| 



= J DC^De'^^- exp I J d*x 
J DC ^5 {d^C^)e^^{j d^x 



2f] 



^Cl + iC^{d^d-2g^B, 



27] 



]-^Cl + zC,{d,e-^--2g^B^ 



(5.1) 
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The constraint (9^C^ = can be uniquely resolved by representing in the form = duhf^^y, 
where h^j,^, stands for an antisymmetric tensor field. Notice, that the number of degrees of freedom 
during such a substitution is conserved, since both of the fields and /i^i, have three independent 
components. 

Then, taking into account the relation (|83|) between ^'^"^s- ^^d S^^,, we get from Eq. (B.l) 

j D6^^^-De'--^- exp j d^x {d,e - 2g„,B^f 




j Dx^,{i)Dh^y exp < j d'^x 



127] 



{B.2) 



In a derivation of Eq. (B.2), we have replaced DO^^'^^' by Dx^(^) (since the surface S, parametrized 
by x^(^), is just the surface, at which the field 9 is singular) and, for simplicity, discarded the 
Jacobian arising during such a change of the integration variable |^. In the literature, the above 
described sequence of transformations of integration variables is usually referred to as a "path- 
integral duality transformation" . In particular, it has been applied in Ref. [ p,l(j|| to the model with 
a global f/(l)-symmetry. 

It is further convenient to rewrite exp 

DGpu exp I / d'^x 
after which the 5^-integration yields 



as 



j DG^^ exp J d'^x [-Gl^ + iG^^^F^^^ 5 {e^yXpO^ {Gxp - 9mhxp)) 



DAp exp\ d X - {gmhpu + dpA^ - d^Ap) + ig^hp^F - 2iApj 



(5.3) 



Here Ap is the electric field, dual to the dual gauge field Bp. The dependence on this electric 
field actually drops out upon the hypergauge transformation (see e.g. Ref. [^]) hp^, — > hp^, — 
J- {dpA^ — d,,Ap). In particular, the interaction of the A^-field with the electric current cancels 
out during such a transformation. The outcome of this transformation together with Eq. (B.2) 
yield Eq. (|8^) of the main text. 



C. Integration over the Kalb-Ramond Field 

Let us carry out the following integration over the Kalb-Ramond field 

1 



Z 



Dhpy exp 



d x\ 



127] 



f72 



'i'^hpi^Yipi, 



(G.l) 



To this end, it is necessary to substitute the saddle-point value of such a functional integral back 
into the integrand. The saddle-point equation in the momentum representation reads 



2r/2 



30 



For the case when the surface S has a spherical topology, this Jacobian has been evaluated in Ref. |115] 
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This equation can be most easily solved by rewriting it in the following way 



iC.2) 



where we have introduced the projection operators [TT 

P^J.u,Xp = - (V^xVup - V^pVux) and l^u^\p = - {Sf^xSup - 5^p5ux) 
with V^u = 5pu — These projection operators obey the following relations 



^vfi^Xp 



'^fj,u,pX ^Xpjpu^ ^ fjLP,Xp^Xp,af3 ^pu,af) 



(the same relations hold for Ppu,Xp), and 



P, 



pv.Xp 



l-P 



Xp,af3 



0. 



(C.3) 



{CA) 



By virtue of the properties (C.3) and (C.4), the solution to the saddle-point equation (C.2) reads 



p2 _|_ ^2 



m? 



(i-P) 



^a(3ip), 



Xv,al3 



which, once being substituted back into the partition function (C.l), yields for it the following 
expression 



Z = exp< —n^T]^ 



d^p 



i + i^(i-p) 



(C.S) 



pi/,al3 



Let us now prove that the term proportional to the projection operator yl — Pj on the R.H.S. 
of Eq. (C.5) yields in the coordinate representation the boundary term. One has 



- P)pu,aP = T:{SupPpPa + Sf,aPuPl3 " ^uaPpPp - ^pPPuPa)- 



iC.6) 



Then, by making use of Eq. (C.6), the term 



2 2 
— TT 7] 



d!^p 



p^ 



(27r)^ p"^ + m? m? 
under study, after carrying out the integration over p, reads 

This is just the argument of the first exponential factor standing on the R.H.S. of Eq. (|85D , i.e. 
the desired boundary term. 
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D. Summation over the Grand Canonical Ensemble of Vortex Loops in 
the Dual Abelian Higgs Model 

In the present Appendix, we shall outline some steps of a derivation of Eq. ( |186| ) of the main text. 
Let us first consider the infinitesimal world-sheet element of the a-th vortex loop, which has the 
form 

Here, in order to distinguish from the index enumerating the vortex loop, we have denoted the 
indices referring to the world-sheet coordinate ^ by = 1,2. Analogously to the 3D case, 
it is reasonable to introduce the centre-of-mass coordinate (position) of the world-sheet = 
J d'^C,Xp{^). The full world-sheet coordinate can be respectively decomposed as x^^^^) = 1/^ + z'^{^) 
with the vector Zp{C,) describing the shape of the a-th vortex loop world-sheet. Then, substituting 
into Eq. ( p.72| ) instead of S^j, the vorticity tensor current of the vortex loop gas ( |185|) , we can 



perform the summation over the grand canonical ensemble of vortex loops as follows. 



7V=l^^- \i=l'' " / na=±l ^ a=l 



Here, is a certain rotation- and translation invariant measure of integration over shapes of the 
world-sheets of the vortex loops. Employing now the dilute gas approximation, we can expand h^j^ 
up to the first order in a/L, where a stands for the typical value of |2"(^)|'s (sizes of vortex loops), 
which are much smaller than the typical value L of |?/"|'s (distances between vortex loops) |^. This 
yields 



^^^(O/iN cos (tt / da^,y{z{^))h^^{x{C})] ~ / Dzpii[z\ cos {-KV^u^x[z\dxh^y{y)) 



E T^^^^^'^'^^i^^'i'^i'^^) ■ ■■^^2r.h^,2„u2Ay) J Dzpfi[z]Vp,^,,xAz] ■■■'Pp., ], [D.l] 



n=0 



where Vpu^xlz] = J dapu{z{^))zx{C)- Due to the rotation- and translation invariance of the measure 
the last average has the form 



/ Dzpfi[z]Vp,^,,xA^]---Vp, 

nl'2n,\2n [^] 



/ '!\2n 



l/il!yi,^2l'2^AlA2 ■ ■ ■ l/i2n-l!^2n-l,A'2n!^2n^A2„-lA2n + pCrmUtatiOUS 



(2n- 1)!! 

Substituting this expression into Eq. (D.l) and estimating the derivative w.r.t. ?/ as 1/L, we finally 
obtain 



31 



Similarly to the 3D-case, the approximation a L means that vortex loops are short living objects. 
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J Dzp{^)fi[z] cos i^n J - E -(2^ \ T) ' = cos y 

where we have introduced a new UV momentum cutoff A = which is much larger than 1/a. 

This yields the desired Eq. ( |186| ). 

E. Calculation of the Integral ([TF^ ) 

In this Appendix, we shall present some details of calculation of the integral (|193|) . Firstly, owing 
to the definition of the functions and , we have: 

J ^'^^^ J (27r)4 J (27r)4 J + + ~ 



flU 



(27r)4(p2 + m2)(p2 + JVf2)' 

Next, this expression can be rewitten as 



It is further convenient to introduce new integration variables a G [0, +00) and t G [0, 1] according 
to the formulae a = at and f3 = a{l — t) . Then, the integration over t yields for Eq. (E.l) the 
following expression: 

(47rm£,)2 y V / 

Such an integral can be carried out by virtue of the formula 

+00 / /? \ f 

J x^-^e-^-^^dx = 2 i^^j (2^^ , 3^/3 > 0, 3^7 > 0, 
and the result has the form of Eq. (|194]) from the main text. 
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